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eo Instationidre Warmespannungen in Hohlzylindern 
mit Kreisringquerschnitt * 


Von R. Trostel 


I. Einleitung. Die vorliegende Arbeit untersucht den Spannungszustand infolge instationarer 
rotationssymmetrischer Temperaturfelder in Hohlzylindern, deren Mantelflachen unter dem Ein- 
flu zeitlich und bis auf die Achsensymmetrie auch értlich beliebiger thermischer Randwirkungen 
| stehen. Die Ermittlung eines Warmespannungszustandes bedingt zunichst die Bestimmung des 
zugehérigen Temperaturfeldes, und daher wurde in Abschnitt II zunachst einiges iiber das achsen- 
symmetrische Temperaturfeld vorangestellt 1, *, . Hieran schlieBt sich in Abschnitt III die Span- 
nungsermittlung in einem endlich langen dicken Kreisrohr an. Das Ziel war die Ermittlung eines 
oberflachenspannungsfreien Warmespannungszustandes, der schlieBlich durch Uberlagerung der 
aus dem sog. thermisch-elastischen Verschiebungspotential 4 hervorgehenden Warmespannungen 
mit gewissen, durch Lovesche Funktionen *, > charakterisierten Randspannungszustanden gefunden 
werden konnte. In Abschnitt IV wird eine Lésung fiir den mittleren Bereich eines sehr langen 
Rohres mit langs der Zylinderachse unabhangigen thermischen Randwirkungen angegeben und 
hieraus der Fall freier Zylinderenden und des ebenen Verzerrungszustandes hergeleitet, sowie auf 
den Fall des ebenen Spannungszustandes hingewiesen. In Abschnitt V finden sich einige Naherungs- 
formeln fiir den Fall eines ebenen Temperaturfeldes, das durch Anheizen des Innenmantels 
entsteht, wobei fiir die Erwarmung des Innenmantels in erster Naherung des Anheizgesetz 
0;(t) = 0;, (1—e—’i') als maBgebend angesehen wurde. Einem Beispiel folgen schlieBlich in 
Abschnitt VI einige grundsatzliche Bemerkungen zum achsensymmetrischen Problem an diimnen 
Kreisrohren (Kreiszylinderschalen), wo auf die Angabe entsprechender Formeln jedoch aus Griinden 
der Platzersparnis verzichtet wurde. 


II. Das achsensymmetrische Temperaturfeld infolge zeitlich und értlich beliebig vorgegebener 
Randwirkungen an den Mantelflichen. 1. Grundlagen der Warmeleitung; die Warme- 
leitungsgleichung. : z 

Die Differentialgleichung fiir rotationssymmetrische quellenfreie Warmeleitungsprobleme lautet 
Bak Bio Ye 


40 =(53 7 ae) 7 ee HS . 


Hierbei sind c (cal/gr°C) die spezifische Warmemenge pro Gewichtseinheit, y (gr/em*) das spezi- 
fische Gewicht und A (cal/cmsec°C) die Warmeleitzahl. 


Gleichung (1) driickt das Gesetz der Erhaltung der Warmeenergie aus und zwar mit Verwendung 
des Fourierschen Warmeleitungstheorems. Hiernach ist die in der Zeiteinheit durch die Flachen- 
einheit hindurchtretende Warmemenge q (cal/em?sec) proportional dem Temperaturgradienten: 

ao oo 
a 203 a} =—Agrad d= —AI7"s 05 at (2) 
Fix die Formulierung der Randbedingungen hinsichtlich eines Warmeiiberganges in ein anderes 
Medium benétigt man noch eine zusatzliche Beziehung, die das Newtonsche Abkiihlungstheorem 
liefert. Entsprechend dem Fourierschen Theorem wird hier die aus dem erwarmten K6rper in ein 


* Gekiirzte Fassung einer von der Technischen Universitat Berlin-Charlottenburg genehmigten Dissertation. 
Der Verfasser ist den beiden Berichtern, Herrn Prof. Dr.-Ing. I. Szabé und Herrn Prof. Dr.-Ing. A. Teichmann 
sehr zu Dank verpflichtet. 

1 B. Riemann, Partielle Differentialgleichungen, Braunschweig 1882. 

2 H. S. Carslaw u. J. C. Jaeger, Conductions of heat in solids, Oxford 1947. 

3 R. Weyrich, Zylinderfunktionen, Leipzig und Berlin 1937. 

4 EF, Melan u. H. Parkus, Warmespannungen stationarer Temperaturfelder, Wien 1953. 
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umgebendes Medium der Temperatur Jy ttbergehende Warmemenge je Flachen- und Zeiteinheit 
TRand = Vand Neier (om— 9) “4 (3) 


also proportional zur Differenz der Temperaturen des Mediums (Pm) und der Korperoberflache (Ix) 
gesetzt. Hierbei ist 1 der Vektor der duBeren Normalen der Kérperoberflache und « (cal/em*sec°C) 
die Warmeiibergangszahl. . 

2. Die Rand- und Anfangsbedingungen. a) Die Randbedingungen an den 
Mantelflachen. Die in der aufzufindenden Lésung enthaltenen Randwirkungen an den Mantel- 
{lichen rr; und r=r, sollen mathematisch formuliert sein kénnen 

1) durch zeitlich und értlich beliebig veranderliche vorgegebene Randtemperaturen O;(z, t) und 
0,(z, t) des Zylinders, 

2) durch zeitlich und értlch beliebig veranderliche vorgegebene Temperaturen (y;(z, t) und 
Onma(Z, t), der im Innern bzw. am Aufenrande befindlichen Medien, . ; 

3) durch zeitlich und ortlich beliebig veranderliche vorgegebene, durch die Mantelflachen des 
Zylinders hindurchtretende Warmemengen q;(z, t) und q,(z, t) sowie 

; 4) durch Kombinationen zwischen 1), 2) und 3). 
Im Falle 1) lauten die Randbedingungen an den Mantelflachen 


Wr: z, t) = D (2, t) ’ Ora, By t) aa Ole, t) (4a) 
und im Falle 3) mit (2) 
, 00 . , oo 
AE ae MDs AFL tld. cH) 


Zur Darstellung der Randbedingungen im Falle 2) kann man zundchst nach (3) fiir die Mantel- 
flache r=r;, indem man (rand = Q(ri» 2, t) = fa(ri » igh) se Ose liek t)t und tt = {m5 Nos n,} 


Sey: 0} setzt QRand tt =— g(r;, 2, t) =— a (Om i(z. t) — O(r;, ee t)) und entsprechend fir 
v od sees i 
die Mantelflache r= rq schlieBlich q,(ra, 2,1) = — 4% (Onra(z, t) — O(ra, 2, 1) folgern, und damit 
erhalt man mit Beachtung von (2) 
, C ~ oo 
aris % 1) = 9 (Ones )— Or ))=— AZ 
oo 
Gel tae Zs t) Se es (Oualz, 1) cory B(ras Oe t)) =A Or + Zot ‘ 
also 
2 od A od 
Presta a = Omelet), — Oo F|, = Omelet). (Ae) 


Wir erfassen nun simtliche durch (4a) bis (4c) beschriebenen Randwertprobleme an den Mantel- 
flachen und damit auch den Fall 4) durch folgende allgemeinere Randwertaufgabe, in der jene als 
Spezialfalle enthalten sind: 


Oo oo 
Gy O (Taste te Cp aren == (25 t) d, O(r;, 2, t) + dye 


rei PHO t) - (5) 


Hier sind c,,c,,d, und d, zunachst willkiirliche Konstanten, und indem man beispielsweise 
¢ = d=1,c,=—d,=0, galz, t) = 0,(z, t) und ¢,(z, t) = 8; (z, t) setzt, folgt sofort der Fall 1) usw. 

Mit den Randbedingungen (5), von denen im folgenden ausgegangen wird, sind demnach beliebige 
Randwirkungen an den Mantelflachen erfaBt. 


b) Die Randbedingungen an den Stirnflachen z= 0 und z =I. Im Gegensatz zu den 
Mantelflachen soll die Lésung fiir das Temperaturfeld hinsichtlich der Stirnf lachenrandbedingungen 
nur folgende Spezialfalle enthalten: 

1) vollige Isolierung, 

2) Warmetibergang in ein Medium der Temperatur Null und 

3) konstant gehaliene Randtemperatur Null. 

Diese Falle werden durch die allgemeinere Randwertaufgabe 


a, Hr) + aS) = 0, — b, B(r, 0,1) + by 


erfaBt, wobei die Konstanten ay, ag, b,, und b, wieder entsprechend zu wahlen sind. 
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¢) Die Anfangsbedingung. Die Lésung soll einen zur Zeit t = 0 gegeniiber seiner Umgebung 
nicht erwarmten Kérper zum Gegenstand haben. Bezeichnet man die anfangliche iiberall als 
gleich vorausgesetzte Temperatur der Umgebung mit Null, so lautet die Anfangshedingung 


O(r, 2,0) = 0. (7) 


3. Allgemeine Bemerkungen zur Ermittlung des Temperaturfeldes. Eine un- 
mittelbare Entwicklung der den Rand- und Anfangsbedingungen (5) bis (7) angepassten Lisung 
fiir das Temperaturfeld ist nicht méglich. Hierfiir mu8 vorerst vorbereitend die Aufgabe gelést 
werden, das Temperaturfeld fiir den Fall zeitlich von t = 0 an konstanter Randwirkungen an 
den Mantelflichen zu bestimmen. Die Anfangs- und Randbedingungen lauten fiir diesen Fall in 
Spezialisierung von (5) bis (7): 


oo ry 
a Orb) +a] =0, %0(,0)+%5) 0, (8) 
obs iT. 0 
ao ad 
Cy O(Tas 2, t) + Cy caer == igl2) d, O(r;, 2, t) + d, Fe == fiz fun te 109) (9) 
(Fs 8,0) == O-. (10) 


Diese Aufgabe wird dadurch bewaltigt, daB man zunachst das Randwertproblem (8) und (9) durch 
die sog. stationire Lésung %,(r, 2) befriedigt, und zur stationaren Lésung dann noch eine Zusatz- 
lésung #,(r, z, t) instationaren Charakters hinzugefiigt, die die bereits durch die stationére Lésung 
erfiillten Randbedingungen nicht mehr verandert, jedoch nun auch noch zusammen mit %,(r, 2) 
der Anfangsbedingung (10) Geniige leistet. In Ziff.4 wird die stationdre Loésung %,(r, 2) 
ermittelt, wahrend Ziff. 5 der instationaéren Zusatzlésung #,(r, z,t) gewidmet ist. In Ziff. 6 wird 
dann durch Uberlagerung von %,(r, z) mit @,(r, 2, t) die Lésung fiir die Rand- und Anfangshedin- 
gungen (8) bis (10) gewonnen, und schlieBlich gelingt in Ziff.7 die allgemeine Lésung fiir zeitlich und 
értlich beliebige Randwirkungen an den Mantelflachen nach MaBgabe der Rand- und Anfangs- 
bedingungen (5) bis (7). 


4. Das stationdre Temperaturfeld. a) Allgemeines. Produktlésung der Warme- 
leitungsgleichung. Im stationaren Falle folgt aus (1) wegen 0%)/at = 0 die Differential- 
gleichung 


Fur ihre in Produktform darstellbaren Lésungen folgen nach EKinfiihrung des Produktansatzes 
O,(r, z) = R,(r) Z(z) gemaB der Daniell Bernoullischen SchluBweise mit der willkiirlichen Kon- 
stanten w, die beiden gewéhnlichen Differentialgleichungen 
dZ, 
dz 


+ w?Z(z) =0, — + = 20, g(r) = 0, (11) 


deren Lésungen trigonometrische Funktionen bzw. modifizierte Besselsche Funktionen Nullter 
Ordnung sind. Sie lauten 


Z,(z) = A, cos Wy 2 + B, sin wz , Ro(r) = Ag Ip(@o 7) + By Ko(mo 7) ; (12) 


wobei die Konstanten A, und B, (k = 1,2) sowie m) den Randbedingungen des vorliegenden 
Problems anzupassen sind. 

b) Lésung der stationdren Randwertaufgabe. Von der soeben gefundenen Produkt- 
losung 


Oy(r, z) = Ro(r) Z)(z) = (A, cos woz + B, sin dp 2) (A, I,(@9 7) + Bo Ko(m r)) (13) 


soll nachfolgend entsprechend den Ausfiihrungen von Ziff. 3 die Erfillung folgender Randbe- 


dingungen gefordert werden: 
oo 09| 
a, Bo(r, 1) + Oy = 5 


=0, 6, A(r,0)+%—= =0, (14) 
1 


7,0 
09| 00,,| 
oe = fa(2) » d, Do(ri» 2) = d,—* 


Ta.® 


= file) - (15) 


Ti 2 


Cy Oo(Ta, 2) + Cy 


Zunichst kénnen die beiden ersten Bedingungen (14) offenbar dann hefriedigt werden, wenn man 
1* 
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das Bestehen folgender Beziehungen fordert: | 
a Z(t as = 0,  ,Z(0) +h Ge] =o. (16) 
oder ausgeschrieben 
A, (a, cos wy 1— dy 9 Sin Wo 1) + By (ay Sin Wy I + a Wy 008 Wy 1) = 0, A, b, + By b,m = 0, (17) 
und aus diesem, in A, und B, linearen homogenen Gleichungssystemen kénnen fiir diese nur dann 
von Null oaeneiece Werte erwartet werden, sofern die Nennerdeterminante verschwindet, 


d.h. die Beziehung 
a, b, — a, b, (18) 
dy by + a, by 
besteht. Gleichung (18) ist eine transzendente Gleichung fur die HKigenwerte won (c= 1,25 Sas) 
die fiir das folgende als bekannt vorausgesetzt werden. 
Aus dem Gleichungssystem (17) folgt iiberdies eine Abhangigkeit der Konstanten A, und B, 


in der Form 


tg Wy l= Wo 


y COS Wy | — ay Wy Sin Wy | 
a, Sin yl + a, Wo Cos Wp |” 


ieee Zh 


so da fiir Z,(z) mit der neuen Konstanten 


Cau : 


1 a, sin yl + a, 1 cos w | 


geschrieben werden kann 
Z,(z) = C,[(a, sin wo 1 + ay Wo COS Wp L) cos Wy 2 — (a, COS Wy 1— ay Wo Sin Wy 1) sin wy z] = C,'Y (wo 2) « 
(19) 

Die den Randbedingungen (14) geniigende Produktlésung nimmt daher nach (13) mit (19) sowie 
mit den neuen Konstanten A = C, A, und B = C, B, die Gestalt 

B,(r, 2) = (A Iy(p 1) + B Ky(op 7)) F(a 2) (20) 
an, und da Y(qw, 2) und damit auch die durch (20) gegebene Produktlésung fiir alle nach (18) 
berechenbaren Eigenwerte w,,, die Randbedingungen (14) befriedigt, kann (20) in der Form 


SAG euea Ve = (An Ip(@on 7) + Bn Ky (Won 2)) (oon 2) = 2 (An Ion(r) + Bn Kon(r)) Yn(2) (21) 
verallgemeinert werden, wobei die bisher noch willkurlichen Konstanten 4, und B, aus den noch 
nicht beriicksichtigten Randbedingungen (15) berechnet werden kénnen. Indem man die Lésung 
(21) in (15) einsetzt, erhalt man namlich 


[An (6 Lon(Ta) + ¢20n Ton(ta)) + Bn (cy Kon(ta) + 2 @0n Kén(ta))] Pal2) = falz) » | 


= Pee (d, Ion(T:) + dg @on Ipn(r;)) Ae B,(d, Kon(ri) + dz @on Kon(ri)) | Pa(z) = fil) » | 


one: Striche die Ableitung nach dem Gesamtargument wo, r bedeuten sollen. 
Um hieraus die A, und B, berechnen zu kénnen, ist es erforderlich, die vorgegebenen Rand- 
funktionen f,(z) und f;(z) zunachst nach den Eigenfunktionen Y,(2) zu entwickeln: 


folt) = SOV), file) = Se Vale) (23) 


n= n=l 


a LMe 


(22) 


Die Fourier-Koeffizienten folgen mit Riicksicht auf die Mie a der Eigenfunktionen zu 


JS fole)Pn(z) dz f file) Pole) de 
0) et eat aos, (24) 
a YW (z) dz f YP? (z) dz 


womit die Entwicklungen (23) als bekannt angesehen werden kénnen. Fiihrt man nun (23) in die 
Gleichungen (22) ein, so erhalt man 


2 [An (c, Ip n(Ta) a Cy OOn In(ra)) -- ‘Be (c, Ko A | + Co Won Kon(Ta)) — c()] (2) = = 9 
B 


py [A, (d, Ipn(r;) + dy Won I9,(r;)) + By (dy Kon(r;) + dy Won Ko, (r;)) — ce] Y(z) = 0, 


Ie 


XXIV. Band 1956 Trostel: Instationare Warmespannungen in Hohlzylindern mit Kreisringquerschnitt 5 


und diese Beziehungen sind offenbar dann erfiillt, wenn sie fiir jedes einzelne Reihenglied n be- 
stehen, d. h. wenn die Konstanten 4, und B, folgendem linearen Gleichungssystem geniigen: 


A,(c, Ip n(ra) “= Cy W0n Ion(Ta)) =| B,(¢, Ko1(ta) + Cy W0n K9,(ra)) = c() ? | 


D j , 25) 
A,(d, Ip ,(7;) Se d, W0n Ipn(r:)) ae B,(d, Ko,(7;) ie d, M0 n Ko n(Ti)) = cl) * | ( 
Hieraus kénnen die 4, und B, mit den Abkiirzungen 
i t ’ il / 
at’) — AC (c, Ko ara aa Co Won Ko ara)) 5 x(a) => An (d, Ko,(r;) + d, Won Ko n(7;)) > 
(26) 
4 1 ; / 1 , 
BO =A. (c Ip n(T.) fe C2 Won Tn(ra)) ’ Be) Fahy ae (d, Io n(ri) =I don I.n(ri)) ’ 
wobei 
A.= (c, Io 2 Ee = Cy M0n I) n(Ta) (d, Ko n(Ti) ++ d, Mon K9,(r;)) 
a (e, Kon(Ta) == Co W0n K9,(ra)) (d, Io n(ri) a d, Mon I},(r;)) 
ist in der Form 
4, = 00) eO-+ af ol), By = Ae) + A eo 1) 


gewonnen werden, und damit ist die durch (21) gegebene den Randbedingungen (14) und (15) 
geniigende stationare Temperaturverteilung gefunden: 


Dor, 2) = 3? [(al e + a 6) Ton(r) + (BO ef + BO ) Kon(r)] Pa(2) » (28) 
n=1 
oder einfacher mit der Abkiirzung 
yr) = (a cO + a4) c@) Io,(r) + (BO c(i) + Bl) c@)) Ko,(r): (28a) 
Bo(r, 2) = SY nalr) Pq(2) - (28b) 
n=l 


5. Die instationidre Zusatzlésung #,(r,z,t). a) Allgemeines. Produktlésung der 
Wiarmeleitungsgleichung. Die maBgebende Gleichung fiir das instationdére Temperaturfeld 
0,(r, z, t) ist die Differentialgleichung (1) und wir gewinnen hieraus fiir die in Produktform darstell- 
baren Lésungen mit dem Ansatz @,(r, z, t) = R(r) Z(z) T(t) die Beziehung 


RZ re Ye f Re 
= ae fe 
i Re GF pe 


(29) 


aus der nach Differentiation hinsichtlich t, z bzw. r mit den zunachst willkiirlichen Konstanten A, 
und A, die gewohnlichen Differentialgleichungen 


i iG oy eke Rig Ri ee 
Faas) Gn RGR ie ee ey 
hervorgehen, wobei mit Riicksicht auf (29) 
w® = Ai + 23 (31) 


sein muB. Die Lésungen dieser Differentialgleichungen lauten 
ae 


2 
T(t) = Ay els) ', Zz) = A, cosAyz+ B,sind,z, R(r) = Ay Jo(Azt) + By Noor), (32) 
wobei Jy baw. Ny Besselsche bzw. Neumannsche Funktionen Nullter Ordnung bedeuten. 


b) Das Randwertproblem des Temperaturfeldes %,(r, z,t). Gemai® den Ausfiihrungen 
in Ziff. 3 wird von der durch (32) gegebenen instationdren Zusatzlésung (sofern man noch Ay = | 
setzt, ohne hiermit die Allgemeinheit einzuschranken) 


o 


2 
B,(r, 2,t) = R(r) Z(z) T(t) = eo z) [A, cos A, 2 + B, sin A, z] [Ay Jo(Ay 7) + Be No (Ae 7)] (33) 


verlangt, daB sie bei der nachfolgenden Uberlagerung mit der stationairen Lésung #,(r, 2) die be- 
reits durch die stationare Lésung erfiillten Randbedingungen (14) und (15) nicht mehr verindert. 
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In diesem Sinne hat man also vo #,(r, z, t) die Erfiillung folgender Randbedingungen zu fordern: 


a0, ad 34 
a (rt) +a,54 = 0, by 8,(r, 0,) + baz], = 9; (34) 
ae Lae) d, 8 th +d OD ses. 0 (35) 

c, O,(Tay Zt) + Cy aril ‘ LUT Ee 25, a : 


Diesen kann mit der Lésung (33) geniigt werden, wenn folgende Bezichungen bestehen: 


dZ ae 
ay ZI) + ay) = 0, — , 2(0) + be al epee (36) 
dR dR) 
c, R(ra) + al 0, 4d, R(r,)+¢4 Glan, == @), (37) 


oder, ausgeschrieben (wobei Striche wieder eine Ableitung nach dem Gesamtargument A,r bedeuten) 


A, (a, cos A, 1— a, A, sin A, 1) + B, (a, sind, 1+ a,d, cosa) =0, (38) 
A; b, a B, by A; a ? 
A, (¢, Ty(Ag Ta) + Cy Ay Jo(Ae Ta)) + B, (¢ N(Ag Ta) + €2 Ay No(s Ta)) =0, (39) 
A, (d, Jo(dg ri) + dy Ag J0(Ae r:)) + B, (d, No(Ao 7:1) + do Ap No(s ri)) — 


Damit (38) bzw. (39) fiir von Null verschiedene 4, und B;,(k = 1, 2) erfullt sind, ist es notwendig, 
daB die Nennerdeterminanten der Gleichungssysteme (38) und (39) verschwinden, d. h. daB 


A,(A1) = be A; (a, cos Ay 1— ay A, sin A, 1) — b, (a, sin A, 1 4- a, A, cos A, 1) = 0 


bzw. 
a, b, — a, b 
Bt oo Fa bel oe 
sowie 
A,(Ay) = (Ch Jo(Ag Ta) + Cg Ay Jog =) (d, No(Ae Ti) + dy Ag No(Az r:)) | (41) 
ae (c, No(Ag Ta) + €2 Ap No(A Ta)) (d, To(Ag 71) + dy hy Jo(Ae r:)) =0 


ist. Hieraus sind nunmehr die /,, (s = 1, 2,...) und A,, (k = 1, 2,...) ermittelbar. Man erkennt 
iibrigens aus der Gleichheit der Eigenwertgleichungen (18) und (40), daB 

Ave = W0s (42) 
ist. Aus den Gleichungssystemen (38) und (39) folgt eine Abhangigkeit der Konstanten A, und B, 
in der Form 


dy Jo(Ag ri) + dg hg Jo(A27:) 
?d, No(Agr;) + dg Ag Noda 7;) ” 


ae a, cos A, 1 —a, A, sin A, 1 
ie 1a, sin A, + a,/, cos A,1 ’ 


p24 


so daB mit den neuen Konstanten 


C= Ais Cy A, . 
tO aysindyt “2 a./, cosy.” “2 d, Ni(er;) + dad, Nizr;) .” 


fiir Z(z) und R(r) 
Z(z) = C, [(a, sin A, 1 + ay A, cos A, 1) cos A, z 
— (a, cos A, I— ag A, sin A, 1) sin A, 2] = C, P(A, z), 
R(r) = Cy [(dy No(Ag ti) + dy Ay Noe ri)) Jo(Aar) 
—  J60ar)-+ 44, Jada) Mller] = Gordian), | 
und damit fiir das Produkt R(r) Z(z) mit der neuen Konstanten C = C, C, 
R(r) Z(z) = C By (A, r) YA, 2) 


geschrieben werden kann. Nun erfiillt R(r) Z(z) die Randbedingungen (36) und (37) fiir simtliche 
nach (40) bzw. (41) ermittelbaren Eigenwerte, so da® wir den Produktansatz wieder verallge- 


| (3) 


hh 
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meinern koénnen 


ie) io) == (Oks he co es) — Oks 2 
Beg = > SG, Gun Va.)e 0) = SF Sodan Vee us) 
s=l1 k=1 enh joesil 
wobei nach (31) 
Wks = Ais + An (46) 


ist. Hier liegt also, ebenfalls wie in Ziff. 4 eine Entwicklung des Temperaturfeldes nach den Kigen- 
funktionen vor, und die in (45) verbliebenen willkiirlichen Konstanten C,, sind der bisher noch 
unberiicksichtigt gebliebenen Anfangsbedingung des Problems vorbehalten. 

c) Das Anfangswertproblem des Temperaturfeldes 3,(r,z,t). Wie bereits in Ziff. 3 
erwahnt, erfiillt die stationare Lésung (28) zwar die Randbedingungen (8) und (9), jedoch nicht dic 
Anfangsbedingung (10), und es mu8 daher von der instationaren Zusatzlésung ,(r, 2, 1) verlangt 
werden, daB sie der Anfangsbedingung 


O,(r, z, 0) =— A(r, 2) (47) 
genigt, d. h. zur Zeit t= 0 gleich der negativen stationaren Lésung wird, so da sich die Wirkungen 


von §,(r, ) und #,(r, z, t) zur Zeit t = 0 in der Tat aufheben. Indem man in (47) die Beziehungen 
(28b) und (45) einsetzt, erhalt man zunichst 


3 e Ge Poult) P(2) = 3 ne(r) Wiz) ; 


woraus fiir jedes Reihenglied s 


5 Cys Box(r) = — 7,(r) (48) 


gefolgert werden kann. Beachtet man nun die Orthogonalitatsrelationen der Kigenfunktionen 


Pox(r) 


Ta Tq ’ 
fr Doi(r) o(r) dr =0 fiir i +k, fr i.e) dr =5 Oe + BODE, ) 


so erhalt man schlieBlich fiir die in (48) verbliebenen Konstanten C,,, indem man Gleichung (48) 
auf beiden Seiten mit r ®o;(r) multipliziert und hernach in den Grenzen r; < r < rq integriert 
Ta 
2 


Ce= — rs(r) Pox(r) dr (50) 
: na, win | ea 


LG 


und damit letztlich fiir die instationare Zusatzlésung #,(r, z, t) aus (45) mit (46) und (50) 


UGy 
e 


ae pe (si) 
e 2! Bir) Be) | . “ 
B,(r, 2, t) =— 2 > > (2h, + HOM rys(r) Dox(r) dr (51) 


Gale ise! , 
6. Das Temperaturfeld fiir die Anfangs- und Randbedingungen (8) bis (10). Dieses 
wird gema8 der vorangegangenen Ausfiihrungen durch Uberlagerung der stationaren Lésung 
9,(r, 2) mit der instationaren Zusatzlésung #,(r, z,t) gewonnen. Dementsprechend erhalt man 


mit (28) und (51) : : 
ie Los \ 
co | °° —(") ‘© ) | 

Or, z, t) ad Bol; z) ae O,(r, 2, t) ea oe ns(r) ae 2 = (22, + pans I ns(r) Do 4(7) dr Ps) y (52) 
sk | k=1 kJ Ir 5 
L e J 


7. Die allgemeine Lésung fiir das achsensymmetrische Temperaturfeld nach Mab- 
gabe der Rand- und Anfangsbedingungen (5) bis (7). Fiir die Bestimmung dieses Tempe- 
raturfeldes wird von der Lésung (52) ausgegangen, die man mit 7(r) = 0Iy.(r) + BO Ko.(r) 
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und 7(r) = a Ips(r) + BO Ko,(r) in der Form 
Ta | 


| _ (@ks\? | 
Oz 1) — > | nO(r) — 2 > : ( J Poul?) ry (r) Dox(r) dr | c) Y,(z) 
c= | fn [r?(@on+ Pople 
i . 


os J 
[ie Te | 
Deedes. 
sap Pi Gel ee, | 
a ()(r) —2 oe fr n(r) Dor(r) dr cl (2) (53) 
a ‘i 2 d awa sion | | | 
J | 


schreiben kann. Sie charakterisiert ein Temperaturfeld zur Zeit t infolge der zur Zeit t = 0 ein- 
setzenden und zeitlich konstanten mittels der Funktionen f,(z) bzw. f;(z) gekennzeichneten Rand- 
wirkungen. Im Falle zeitlich veranderlicher Randwirkungen Pa(z, t) baw. p;(z, t) sind nun die 
Fourier-Koeffizienten c bzw. c® nicht mehr konstant, sondern Funktionen der Zeit, was sofort 
aus den Entwicklungen von ¢,(2, t) bzw. q,(z, t) nach den Eigenfunktionen Y7,(z) hervorgeht 


COW (z), (54) 


rel Pal t) aa 2 fa) ¥ (2) : pilé t) ms ah 


seed s 


1 
gale! deg, J Gales t) Vale) de 


1 
ES t + uy Pilz, t) Ps(z) dz 
) (2) AG) oe === ¢) 
to dog dt a is } W%(2) de cl) (t), 
0 s 


iy elementes dr die Zusatzwirkungen 
S OPa i 
Se 0; Pa(%,T) = a dr bzw. Ot Pilz, T) = we dt 
| - d. h. die Zuwachse 
t = (2) de) 
ABD dc = Sire dt bzw. dc = pet dt 
5 dt 3 dt 


Betrachten wir zunachst vereinfachend allein Anderungen am auBeren Rande, und nehmen wir 


an, daf} die zur Zeit T entstandene Wirkung dc(r) von nun an zeitlich konstant weiterwirken | 
wiirde (siehe die obere Figur der Abb. 1), so wiirde das hierdurch zur Zeit t (nachdem also t— 1 | 


Sekunden nach Einsetzen der Randwirkung dc) vergangen sind) vorhandene Temperaturfeld 


mittels der Lésung (53) dadurch erfaBt werden kénnen, daB man in ihr c) durch de® sowie t | 


durch t— 7 ersetzt, und c = 0 setzt: 


Wks \2 ke 
— oo. (" (t—t) (a) 
dj(r, z, = WO (r) 2 g DoyAr) (a) dec\ 
2 L 2 [r? (D5, ap Do; Nee Ms (") ‘i Pox(r) dr P.(2) dt dr ; 
Ti al 


Die Gesamtwirkung zur Zeit t kann man sich nun entstanden denken durch die Summation samt- 
licher wahrend der einzelnen Zeitintervalle dr entstandenen Elementarwirkungen dc zuziiglich 
; | 


der Anfangswirkung c((0) (Abb. 1), womit man 
iT 


4 Ta | 
Ee as | 


Ha) sh = = (a) _— ~ AL @. Dy u(r) a | a 
(r, 2, t) > (r) Ps, (2 @a,-4 oy 0)(r) r Dox(r) dr | ,(z) (0) 
L Ti 4 
t [ Tq | 
i a = (“) 2 (=) | 
: , 6 a D (2) 
4 PS f(r) 2 >” Sout) f(r) + Box(r) dr | ¥(2) Se 


GS + opr 
0 | ev | 
= Ly i} 


und es entstehen nun wabrend eines beliebigen Zeit- _ 


at 
t of) = 1 = c(t), 
JS ¥2(z) dz 
Le ig t-t = : (55) }} 


| 
| 


| 
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also 


oC 


IMT, 2,0) = S| nr) (t) 


a 


. (r) 5 zo( a)? i * es) dit : 
=a DN [r2 (@ a @52)) c((0) e ( a val e ( a a at | W(r) r Doi (r) dr | P,(z) 


0 7 


erhalt. Beachtet man nun noch die Relation 


t 
(*#)* (= ke) 
a —T) (a) Oks (t—7) 
Oe fot se dz = c(t) =a | cW(t) e aaa dr, 
0 


die sich mittels der partiellen Integrationsregel verifizieren laBt, so folgt 


G(r, x, t) =D) fw —2 See We J r ni(r) Bys(r) ar c(t) ys(2) 
r 0k "6 


STN) 2 Ghee Pork) Pole) ie (a) pies Bigat! 
Papas [r? (D2, + Di Die {| Cs (7) UR (r) r Do x(r) € dr dr . 
0 7; 


Nun verschwindet jedoch das erste Glied dieser Formel, da der in diesem auftretende Summen- 


ausdruck 
Pr fe 
2», [r? (3, + Oi on |”! ea! 


gerade die Entwicklung der Funktion 7%(r) [vgl. (48) und (50)] nach den Zylinderfunktionen Dy ,,(r) 
darstellt, und damit verbleibt 


o 2 
ks 
fo, ] —— — 
ES) (=) 


G(r, z, »=>> 2 Wis Pox(r) Vole) | | (zx) (r) r Dox(r) | 


AEP + ODI | 


dr dt. 


Ersetzen wir in dieser Formel c{(rt) durch c)(r) sowie n@(r) durch 7(r), so haben wir die Lésung fiir 
das Temperaturfeld bei alleiniger Innenrandwirkung, und schlieBlich gewinnen wir fiir das Tempe- 
raturfeld infolge des Einflusses zeitlich und értlich beliebig verdnderlicher, achsensymmetrischer 
Randwirkungen an beiden Mantelflachen 


t Ta ks ‘ie 
2 wt, Dy (r) Pelz) es are ie lira) 
Br, 2,1) = >, » Se roars | | (ote) ni) + ee) nf] + eal dr dr 


= EO 7; (56) 
bzw. mit den Abkiirzungen 
ie ; 2 wi, 
#1 a + DME 
— (Hs) the (57) 
Ty.(t) = if [e(9(x) m.(r) + er) 1f(r)] r Doulr)e dr de | 
einfach 
iis, > Dus Poy (r) 2 .(z) Test) (58) 


enh eh 


III. Der achsensymmetrische Warmespannungszustand am dicken Kreisrohr. 1. Allgemeines; 
das thermisch-elastische Verschiebungspotential. Das im Abschnitt II behandelte 
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achsensymmetrische Temperaturfeld verursacht einen achsensymmetrischen Warmespannungs- 
zustand. Simtliche Spannungs- und Deformationsgréfen sind somit unabhangig vom Polarwinkel g, 
so da® sich die allgemeinen Spannungs- und Deformationsgleichungen bereits erheblich verein- | 
fachen, zumal die Rotationssymmetrie das Verschwinden der Tangentialverschiebungen wu, sowie | 
der Schubspannungen 7,,, und t,, bedingt (Abb. 2). 

Um zu den fiir die Berechnung des Problems maigebenden Beziehungen zu gelangen, werden die | 
Hookeschen Gesetze sowie die Gleichgewichtsbedingungen am Element herangezogen. Die Massen-— 
beschleunigungen infolge des zeitlich verander- 
lichen Temperatur- und damit Spannungsfeldes 
bleiben unberiicksichtigt, d.h. es werden zu 
jedem Zeitaugenblick statische Verhaltnisse 
vorausgesetzt. 

Es verbleiben im rotationssymmetrischen 


Falle die Gleichgewichtsbedingungen 


2 no) OF | 
(59) 

0 0 

5, (0 Oz) ta ( tra) = 0 | 


und die thermisch-elastischen Bedingungen 


dv 1 
& =a = Flo, — 7 (Op +0,)] tu, 
v iL 
Ey = 7 = Te [o, y (o, o:)] 1 Of D d 
60). 
ow 1 ae 
eee al 
Vo» dw | ov = zt on ae 
Abb. 2. cee or , On : rs 


wobei a, den Warmeausdehnungskoeffizienten bedeutet und mit v bzw. w die Verschiebungen in 
radialer bzw. axialer Richtung bezeichnet werden. Aus diesen 6 Gleichungen kénnen die 6 Un- 
bekannten o,, 0,, 0,, Tz > v und w berechnet werden, wobei hier, entsprechend der geringeren Anzah| 
unbekannter VerschiebungsgréBen der Weg der Elimination der Spannungen beschritten werden 
soll. Demgema8 werden mittels der Deformationsgleichungen (60) zunachst die Spannungen in 
Abhangigkeit. von den Verschiebungen dargestellt : 


i ey v Ou v Ow 1l+y 

i | ie ee | Pele 

0; arene sor; (5; al 4] Teor, 00). 
ee st ek PL) _ (ev , vv, Ow i) Say 8 

P ee ape r Toe Ap. le aye As =a?) | 
rc E Ow ee CS ow i is (61). 
LER | cae Teale arr ae pany] 
aes E Ov | Ow 

2) oe | 


. . . . . . | 
und nunmehr in die beiden Gleichgewichtsbedingungen (59) eingesetzt, womit zwei gekoppelte 
partielle Differentialgleichungen fiir die Radial- und Axialverschiebungen v und w entstehen: 


hoe <2) a DP oe 02 1 oO fow Ov 00. 
a ee ( ant 0) ay 3a Tay oe cE . 3 Us ores e (028 
ae i fll 0 \ {ov Ow Ob 
(l v) Aw | 5) E | x) bs a (d | v) Xs Eye =" (), (63) | 


Zur Lésung dieses Gleichungssystems wird das _ thermisch-elaStische Verschicbungspotential 
F(r, 2, t) eingefiihrt, das in folgender Form durch seine partiellen Abteilungen definiert wird: 
Seen peels 


v=, p= ; 
i i (64)) 
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Mit diesem Ansatz gehen die beiden Verschiebungsgleichungen in 


0 1ty a ly 
or AF, — 758] =0, [4h — jt iae — 0 
uber, woraus fiir F\(r, 2, t) die Poissonsche Differentialgleichung 
a ee en 
AF, = jae ey) OD (65) 


gefolgert werden kann. 


Ist zunadchst ohne die Beriicksichtigung der vorliegenden Randbedingungen ein partikulares 
Integral der Gleichung (65) gefunden, so ergeben sich nach (61) und mit Beriicksichtigung von (65) 
folgende, i. allg. den Randbedingungen nicht geniigende Warmespannungen: 


Ea: : E @F. Eins « E 1 OF 
(a9) 1 (i) ee ere. eek oe Sepa 
eo ] -- ti Pit, 7 | oe 1—y l+yvr or’ 
a 2 (66) 
(1) Ear 9 of OF, sO. pohie OLR y 
a =) l+y dz? re 1 + » Ordz” 


2. Der Randspannungszustand; die Lovesche Funktion. Da der nach (66) berechen- 
bare Warmespannungszustand J i. allg. die Randbedingungen nicht erfiillt, ist diesem ein zusatz- 
licher Spannungszustand zu iiberlagern, der die infolge des Zustandes 1 hervorgerufenen Rand- 
spannungen an die tatsachlichen Gegebenheiten anpafit. Dieser Zusatzspannungszustand rihrt also 
allein von Oberflachenspannungen her und seine weitere Verfolgung ist somit ein reines elastizitats- 
theoretisches Problem. Die maBgebenden Gleichungen sind demnach durch (59) und (60) bzw. durch 
(62) und (63) gegeben, sofern man # = 0 setzt. In diesem Falle fiihrt man zur Lésung der Verschie- 


' | bungsgleichungen (62) und (63) die Lovesche Verschiebungsfunktion L(r, z, t) ein. Man definiert sie 


in folgender Form durch ihre partiellen Ableitungen: 


Tee é VA aL 
ei a ae so pl) SE ale 


(67) 
Hiermit kann man die erste Verschiebungsgleichung (62) identisch befriedigen, wahrend aus der 
zweiten Gleichung (63) 

ALL (68) 


hervorgeht. Sind entsprechende, den Randbedingungen des Problems geniigende Lésungen Ey (7, 2; 0) 
der Gleichung (68) gefunden, so folgt mit (67) aus (61), worin ebenfalls # = 0 zu setzen ist, fiir die 


Spannungen 
Os a+ 0 ») (? AL. — 2), 
= (lee xa —2) 5 (e z) AL. — 3); | 
ay = a+») “i —2) 5 (a— y) Aly, — oo 


3. Das partikulare Integral der Poissonschen Differentialgleichung. Die Auf- 
findung der partikularen Lésung F,(r, z, t) wird durch die Struktur der Lisung (58) fiir das Tem- 
peraturfeld auferordentlich vereinfacht. Zerlegt man die Funktion F, in die Teilfunktionen F';, 
in der Form 


R= 3 2 Fu, (70) 
so kann (65) mit (58) auch in der Gestalt 


ee a, Dye P or) VG) Ts, (2) (71) 


sare ica a 5 i 
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geschrieben werden, und dieses ist immer dann erfillt, wenn fiir jedes Glied der Reihe | 


AF ue = 52 0% Dus Pox(r) Pz) Tes(t) (72) 


besteht. Zur Lésung von (72) wird mit der zunachst unbestimmten Konstanten y; 5 
Firs = Yrs Dor(r) V(z) Th s(t) (73) 


gvesetzt, und da fiir die Zylinderfunktionen ®o;(r) = Do(Ax, 7) und die trigonometrischen Funktionen 
W(z) = (Ai, z) nach (30) die Differentialgleichungen 


LD 1 d® F d? 
sah 4 1 eh — 78, us, 


: Soe iis Y, (74) 


bestehen, ergibt sich, sofern der Ansatz (73) in (72) eingefiihrt wird, mit Beachtung von (74) schliefi- 
lich fiir die Konstante 
oo 1 +- v “ Di 1 -- v Dx s 
Uke oN gt 82 ee eae Dey ane 


und damit fiir das partikulare Integral der Poissonschen Gleichung 
SEN L+ 9 O&% QTY Q7 Gon) Me) Tes | 
= F, eS 2 ny r, (75), 
spo a eee | 
Fur die Spannungen erhalt man also gemaf (66) mit (58) und (75) 
2 Eo, ae s (a + Fx) Port) + Be Ohl) ¥@) Te s(0) 
s=1 


GOT, Z, 8) = i 


[v2 (3, + BEN | 


2E 4, >; [Ais + 28x) Pole) + Ca!) Pi eM Vo@) Tes) 
— @ (ly) A [r? (Pi, + Doe 


. 2 Eo, Ax Por) Vs(z) Ty 
o(r, ad) t) z= a eee (1 se ”) s = [r? (P? k ate Pi ADK ° 


2 Boy SS) Ays Anz Pox(r) P3(z) Ty s(t) 
a (L—v) oy 2 (DR, + DERI 


a? (1 —») 


or, Zt) 


Mr, 2, ) = 


und insbesondere fiir die Oberflachenspannungen 


a) an der AuSenmantelflache r — r, 


25 ey > » (Ar «+ Abe) Por (ra) st Ao k Py Ka) | ¥(z) Ty s(t) 


o)(rq, 2, t) = — 


a (1 — ») [r? (Gx + Pople 2 aw 
7 
TO(Ta, z,t) = ites >> S71 An sha t PD Be) =) Tus) | 
s=l k=1 [r? (®o, + Poi)1," 
b) an der Innenmantelflache r = r; 
AM(r;, z, 1) =— coe > Mais + Bi) Poulri) + Man Ponls)] ze) Tet) | 
OE) eat eie , FeRt ae (78) 
2Eo Das Aan Phy ri) Vi(z) Ty 9(t) | 
(r;, 2, t a ae 1s2k Pok ks(l) 
ral Ye Ya? d=0) > Da [r2 (2, + ay 3 
c) an der Stirnflache z = | 
2B a 23, Dyi(r) W(L) Ty (t) 
o(r, i f= e ok s s | 
é ~ a@(l—y) > a [r? (@2,, + @ ale ‘ : 
| 


(1) Re Pie Ar shar Pordr) Pl) Tr, s(t 
WNC cera e eh Feibie 


9 
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d) an der Stirnflache z = 0 


o)(r, 0, t) =— 2 Eo, SS STARE Por(r) Y(0) Tr s(¢) 

Oe ec ads Oe a (80) 
1(r, 0, t) ———2 S ST Aas Mak Ph i(r) ¥4(0) Ti s(t) 
i pS) pe! (03, Ochs 


4, Die Spannungszustinde zur Erreichung randspannungsfreier Oberflachen. 
In den nachfolgenden Untersuchungen wird es sich darum handeln, mit Hilfe bestimmter zu iiber- 
_lagernder Randspannungszustiinde resultierend einen randspannungsfreien Warmespannungs- 
zustand zu erzeugen, d.h. die durch (77) bis (80) gegebenen Oberflaichenspannungen zum Ver- 
if schwinden zu bringen. Zunachst wird dieses fiir die Mantelflachen r = r; bzw. r = Ta durchgefiihrt. 

a) Der Randspannungszustand 2 zur Erreichung randspannungsfreier Mantel- 
flachenr—=r; baw.r=r,. Die maBgebende Beziehung ist die Differentialgleichung (68) fiir die 
Lovesche Funktion L,(r, z, t). Zur Lésung von (68) wird zunachst 


AL, = L,(r, z, t) (81) 
gesetzt, womit Gleichung (68) die Form 
t AAL, = A(AL,) = AL, = 0 (82) 
/ annimmt. Fir L,(r, z, t) wird nun ein Produktansatz in der Gestalt 
L,(r, 2, t) = SS [Pins Tos(r) + pres Kos(r)] We(z) Tea(t) (83) 


sai k=l 


benutzt, wobei pi;,; bzw. p2;, willkiirliche Konstanten und Y{(z) die Ableitungen der durch (43) 
definierten trigonometrischen Funktionen nach ihren Gesamtargumenten /;, 2 bedeuten. Dieser 
Ansatz befriedigt wegen 

he d? 


| 5 i (2)] =— A, Pe) 


a? ld 
(is i ra, AG eee 


: die Differentialgleichung (82) identisch, was durch Einsetzen verifiziert werden kann, und man hat 
demnach fiir die Lovesche Funktion nach (81) die Differentialgleichung 


| AL, = = = [Pixs Ip s(r) + prs Kos(r)] Ps(2) Tr s(t) - (84) 

Zu ihrer Losung wird mit der noch unbestimmten Funktion R,,,(r) der Ansatz 

! YEAR ey = Tint aA Tacs 4) — = = Reo) 97, (2) Ts) (85) 
s=1 k=1 s= = 


versucht, und indem man diesen in (84) einsetzt, kann man schlieBlich fiir jedes einzelne Reihenglied 


| (k, s) 
| ALons = A [Res(r) Vi(2) Tes(t)] = [Pine Ios(r) + pers Kos(r)] Y5(2) Trs(t) (86) 


i folgern. Hieraus erhalt man schlieBlich unter Beachtung von 


© Wi (2)] = — Bs) 


! folgende gewohnliche Differentialgleichung fiir R,,(r): 


Apr p dR ohh 
| et ta Ais Ra slr) = Pats Los(r) + Pare Kos(y) - (87) 
) 
} 


Thre Lésung setzt sich aus dem aus modifizierten Besselschen Funktionen bestehenden Fundamental- 
system der homogenen Gleichung sowie einer partikularen Lésung der inhomogenen Differential- 


| * Diese Differentialgleichung wird natiirlich auch von K,(r) erfiillt. 
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gleichung zusammen. Man erhalt schlieBlich das einfache Resultat 
Res(r) = 54 (Pres Los(r) + Paes Kos(r))-+ pans Tos(r) + pars Kos(r) » (88) 
1s 
so daB fiir die Lovesche Funktion nach (85) 


L251) = > > Ore [Pixs Zos(r) + poxs Kos(r)] + pans Los(r) + pars Kos(r)} Us (2) Txs(t) (89) 


sl h— 


gefolgert werden kann. Damit ergeben sich fiir die Radial- und Schubspannungen nach (69) 


E SNe D 
Ta FHC =H 2 2 Paks [Ars t Hoe + (129) Toa] 
=f P2ks [Ais Lf Kos ae (l— 2 v) Kos| aie 2 P3ks his Ios a 2 Paks Ais Kos} As P,(2) Ty, (t) M4 


E [2,2 co ; 
w= 3 a+» Sys {Piks [Ast Jos + 2 (1—1) Ios] 
+ P2ks [Ais r Kos an 2 (1— v) Kos] ai 2 P3ks Ais Ips = 2 Paks Nis Ko.}Ais VG) T;, s(t) 


(90) 


und speziell fiir die Randspannungen an den Mantelflachen 


E ei j 
Oro Fat) = ee Coa pa > {Dies is Teledt) il ae ye 
+ Prks [Ais T; K6 (ri) ele (l— 2 v) Ky ;(r:)] a 2 P3ks i I(r) 
ae 2 Paks ic K@.(r;)} Ais T(z) T;, <(t) ? 


(91) 
Wr) = re Oe) Dy (Pres ners odes) +209) Toate) | 
| 


a P2ks [Ars ue Ko s(r;) a Z (1 = v) Ko(r;)} ++ 2 P3ks his Io 5(r:) 
Sle 2 Paks Ais Ko,(r;)} Ais YW, (2) Ty, s(t) ’ 


E wh Se / 
GP)(rq, 2) t) = ome ie 1) al any 2 Da ay {Piks [Ais Tq Iq (Ta) i (ee 2 v) I) (Ta) | 


+T P2ks [Ais Ta Ko (Ta) se (1 —2 v) Ko .(r.)] =— 2 P3ks His Ip 5(ra) 
Spa Paks Ais 05(Ta)} As a) T;, s(t) $ 


ane y (92) 
E : : , 
Nie % ) = seas spy {Pies [Ars ra Too(Ta) + 2 (1—») Ioa(ra)] 


T P2ks [Ais Tq Ko s(ra) =| 2 (1 ra v) Ko (ra) | =e 2 P3ks Ais Aes) 
ai 2 Piks His K9,(r.)} Ais 17. (z) Ty s(t) . 


Nunmehr kénnen die bisher willkiirlichen Konstanten Puks (4 = 1, 2,3, 4) aus der Bedingung 
der Randspannungsfreiheit der Mantelflachen berechnet werden. Die entsprechenden Randspan- 
nungen des Zustandes (2) miissen demnach gleich den negativen Randspannungen des Zustandes (1) 
sein, also | 


Mr 1) =I Mr at), Or) =—— Ar 2,1), | 


TZ 


0) (ra, z, t) = — 0) (ra, Z, t) p tO. 5 be oO) (ra, Z, t) ‘ | (93) 


Mit (91) und (92) sowie (77) und (78) entspricht das Gleichungssystem (93), sofern man seine Giiltig- 
keit fiir jedes einzelne Reihenglied (k, s) fordert, einem linearen Gleichungssystem fiir die Kon-. 
stanten p,,x5 (44 = 1, 2,3, 4), aus dem diese ermittelbar sind. Die Lovesche Funktion (89) ist damit 
vollig bestimmt, und mithin sind auch samtliche Spannungsgrofen des Zustandes 2 gegeben. | 
Mit den Spannungszustanden J und 2 ist allerdings noch nicht die Spannungsfreiheit der Stirn- 
flachen z = 1 bzw. z = 0 erreicht. Hier verbleiben Schub- und Normalspannungen, und von den 
Schubspannungen kann man, da sie ja an den Mantelflachen bereits zum Verschwinden gebracht 
worden sind, von vornherein sagen, daB sie fiir r — r; bzw. r=rg Null sein miissen (Abb. 3), so 
daf sie im Zwischenbereich sicherlich keine grofen Werte annehmen kénnen. Da sie dariiber hinaus 
wegen der Rotationssymmetrie des Problems an jedem Rande ein Gleichgewichtssystem bilden, 


| 


4 


| der Zylinderenden die tatsachlichen Span- 
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kénnen sie nach dem De Saint Venantschen Prinzip in einiger Entfernung von den Randern 
kaum noch einen Einflu® auf den Spannungszustand des Rohres haben. Dasselbe gilt fiir die 
Axialspannungen. Da namlich in den bisherigen Lésungen die Méglichkeit verschiedener wirme- 
_ technischer Gegebenheiten an den beiden Stirnflaichen enthalten ist, so ist der Verlauf und die 
GréBe der von den Zustanden ] und 2 herrithrenden Axialspannungen an den beiden Stirnflaichen 
i. allg. sicherlich verschieden, so daB dementsprechend auch zwei voneinander (betrags- und vor- 
_ zeichenmafig) verschiedene Axialspannungsresultierenden P,) bzw. P,, erwartet werden kénnten 
_ (Abb.4). Nun muB sich jedoch andererseits, da bei den vorangegangenen Betrachtungen in jedem 
_ Falle vom Gleichgewichtsprinzip ausgegan- 
gen wurde, der Zylinder unter dem EinfluB 
der Spannungszustinde J und 2 im Gleich- 
- gewicht befinden, und insbesondere fordert 
die Gleichgewichtsbedingung der axialen 
Krafte mit Beriicksichtigung der Tatsache 
randspannungsfreier Mantelflachen 


Peep awl) < (94) 


Dies kann nach dem vorangegangenen nur 
dann erfillt sein, sofern P.)» =0 und 
P.,=0, d.h. sofern die Spannungsresul- 
tierenden jeweils selbst verschwinden. 
Der durch die Spannungszustande I und 
2 beschriebene Wdarmespannungszustand 
erfaBt dementsprechend bis in die Nahe 


nungsverhaltnisse nahezu exakt, so daB die Abb. 3. Abb. 4. 


| bisherigen Uberlegungen fur den Fall nicht 


zu kurzer Rohre praktisch hinreichend sind. Um jedoch auch noch den Fall sehr kurzer Rohre zu 
erfassen, sei nachfolgend versucht, den Weg fiir eine Lésung des Stirnflachen-Randwertproblems 
aufzuzeigen. 


b) Die Randspannungszustande zur Erreichung randspannungsfreier Stirn- 
|flachen. « Allgemeines. Wir sehen fiir das folgende die Spannungen der Zustande J und 2 
| als bekannt an. Sie sind Funktionen des Ortes und der Zeit, wobei der zeitabhaingige Teil Ty, ,(t) 
} mit den ortsabhangigen Groen jeweils multiplikativ verbunden ist. Die durch die Zustande / und 2 


| charakterisierten Spannungen kénnen also in folgender Form dargestellt werden 1: 
| 


Gor, (2. be > Dy GUA) (r, z) Ty.(t) , ol t)(r, z,t) = oy Sy GUL+2(r, z) Tys(t) 5 
. are s=1 k=1 (95) 
| alt)(r, 2,1) = SS Mr, 2) That), Nr, 2, 1) = SA Mr, 2) Trolt) - 
if s=lk=1 — s=1 k=1 
'| Speziell erhalt man hieraus fiir die Spannungen an den Stirnflachen 
| = <7 oo os) 3) 
ort) = SS or, I) T,.(2), 02, 1,1) = SD r(x, 2) T,.(2) , 
| , lbs. aL eat (96) 
| or, 0,2) = 3 SGlsMr, 0) Tr), Nr, 0.) = SY LI, 0) Teall), 
| ; s=1 k=1 R= il peal 


und es handelt sich nun um die Aufgabe, einen Randspannungszustand, (d. h. eine Lovesche Funk- 
tion) aufzufinden, der an den Stirnflachen den Randbedingungen 


o,(r, 1, t) = ss ; G,49(7 1) Ty. s(t) =— of V(r, I, t) = — Piee Ut (r, 1) Trs(t) » 
s=l k= a yo 
OAT, 0, t) = > y Oe Al, 0) Tie) a or, 0, t) =a > = o+2)(r, 0) T;. s(t) r) 
goal b= =i 


1 Zur Vermeidung unnotig vieler Indizes wird 1 anstelle von T,, geschrieben. 
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j (97) 


u(r, l, t) = Ss Ty o(Ts l) Ty, s(t) =— rl F2)(r, l, t) a P54, Pa TUT) (r, yD) Ty, s(t) ’ 


or, 0,1) = 3) 3 G.(r, 0) Tealt) =— e+ %(r, 0,1) =— YS AH, 0) Test 


geniigt, ohne hierbei jedoch die bereits an den Mantelflachen erfiillten Randbedingungen hinsicht- 
lich der Radial- und Schubspannungen zu verandern. Neben (97) miissen also noch zusatzlich 
folgende Randbedingungen gefordert werden: 


OATas 2, t) =Q) ? Or(Tis = t) = 0 o) Tl Fas z, t) = () ° t(Tis 2) t) = Oe (98) 


Die Bedingungen (97) und (98) sind erfiillt, wenn fiir jedes einzelne Reihenglied (k, s) folgende Be- 
ziehungen bestehen: 


6,45(T, t) =— ot 2)(r, LG. O43 F Oa ol t2)(r, Oy; | 

= bs & Po 9 

Tadlts 1) = — WENT). T(r, 0) =~ HLH, 0), | ve 
OeAtas z) = 0 ’ Onna lis 2) =a) o) TAs 2) =U) ’ Cis Tes z) =0i (100) 


Wie bereits aus (99) und (100) ersichtlich, separieren wir fiir das folgende die Zeitfunktionen Ty, ,(t)> 
womit die Aufgabe auf ein Randwertproblem fiir die quergestrichenen Groen 0 und T hinauslauft- 
Die Untersuchungen beschrianken sich auf die Betrachtung eines Reihengliedes (k, s). Die Ver- 
allgemeinerung auf siémtliche Reihenglieder bietet keine Schwierigkeit. 

Die Ermittlung des Stirnflachen-Randspannungszustandes erfolgt in zwei getrennten Arbeits- 
gingen, denen zwei voneinander unabhangige Lovesche Funktionen L,(r, z) und L,(r,z) zugrunde | 
liegen. In beiden Fallen werden zunachst die Randbedingungen (100) befriedigt und nach Uber- | 
lagerung beider Funktionen schlieBlich auch die Bedingungen (99). i 


6) Der Randspannungszustand 3. Die mafgebende Gleichung ist die Differential- | 
gleichung (68) und zu ihrer Lésung wird entsprechend Absatz a) 


| 

ALG) = Le) (101) | 

gesetzt, so daf} aus (68) | 
AAL() = A(AL®)) = ALG) = 0 (102) | 


hervorgeht. Die Lésung von (102) wird nun jedoch nicht wie in Absatz a) durch modifizierte || 


AX Risn(r) (An etm * + By etm *) = S [qin Jon(t) + gon Non(r)] (An e™* + B, e-¥m) (106) || 


ot n=1 


gefolgert werden kann, und dies ist erfiillt, sofern man das Bestehen dieser Differentialgleichung fiir \ 
jedes einzelne Reihenglied (n) fordert. Es ergibt sich schlieBlich folgende gewohnliche Differential- 
gleichung fiir Reeth): | 


| 
Besselsche — und trigonometrische Funktionen, sondern durch Besselsche — bzw. Neumannsche || 
sowie Exponentialfunktionen dargestellt. Es ist namlich mit den zunachst willkiirlichen Konstanten || 
ins Gan» An und B, sowie Un | 
! 
= oo 
LET, 2) = = [qin Jo(ln 7) + gan No(un 1)] (An en * + B, e~ Hn #) (103) | 
eine Loésung von (102), was durch entsprechendes Einsetzen leicht verifiziert werden kann, und man | 
erhalt damit nach (101) folgende Differentialgleichung fiir die Lovesche Funktion: | 
AL) _— = (1 n Jo(Un 1) + Gan No(Un r)] (A, e!n* + B, emsn “), (104) | 

Zu ihrer Lésung wird mit der zunichst willkirlichen Funktion R,, ,(r) | 
Tt = 2 Ul, = 2 Resn(r) (dn ot? + By e—Hn*) (105) | 
gesetzt, womit aus (104) mit Jo(un r) = Jon(r), No(un r) = Non(r) 


Rin 1 dRy sn | 
Wd 3 AE: 4 [i Risa = Gin Jon(r) = q2n No n(r) ’ (107) i 
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fiir deren allgemeine Lésung man letztlich nach einiger Rechnung mit den willkiirlichen Konstanten 


dun (uU = I, 2, 3, 4) 


LF / ’ 
Ry sn(r) os Jon =f an Non) = 930 Jon + Yan Non (108) 
erhalt. Fiir die Lovesche Funktion folgt dann gemaB (105) 


- r 1 1 
L®) aa SS ae (1 n Jon = q2n Non) + Wan Jie air Jan Non (An ebn # + Inte e7 En )) (109) 


n=1 
und fiir die zugehérigen Radial- bzw. Schubspannungen gemab (69) 
E = ' 
we = eo a as ye RET) ap {qin [Ln ie Fe ta (1— 2 v) Jon] a Gan [Un if Non 
a (1— 2 v) Non] ae 2 Qn We Jon + 2 Gan Non} Un (An Ce B,, @ Pn a) ” 
s E = , 
We Ey dn) dy {Me [tn Jon — 2 (1) Jn] + don [pln 7 Non 
rte 2 C v) Non] Sif 2 q3n be. Jon ie 2 G4n ie Non} Ln (An Gene ae B, en z) . 


Die Randbedingungen (100) werden von diesen Spannungen auf jeden Fall dann erfiillt, wenn ihnen 
jedes einzelne Glied der in (110) vorhandenen Reihenentwicklungen geniigt, und hieraus folgern 
wir folgendes lineares homogenes Gleichungssystem fiir die q,,,(u = 1, 2, 3, 4): 


{qin [Un T Jon + (L—2 v) Jon] + Qon [unt Non + (L—2) Non] 7 | 

Si 2 Qn it Jon aie 2 dan ING le = 0 ) 
{din [unt Jon + (L— 2%) Jon} + dan [ant Non + (1-2) Non] 

1 2 Gon, Jon + 2 Gin Le Noatra ee Boe 
{din [tn 7 Jon — 2 (L—v) Jon] + dan [Un 7 Non — 2 (L—¥) Non] | oe 


(110) 


= 2 g3n fe Jon 36 2 gan Nowe ie =0, 
{qin [Un r Jon— 2 (l1— vr) Jen T G2n Wis r Non— 2 (1— ») Non] 
4 2 dan be don + 2 On pe (Noni Ok, J 


dessen Nennerdeterminante Null sein mu}, damit wir von Null verschiedene Werte Gun erwarten 
kénnen. Also ist 


A(u) = 0 (112) 


die Eigenwertgleichung zur Berechnung der y, (n = 1, 2,3...), die wir fur das folgende als bekannt 
voraussetzen wollen. 

Allerdings hat das notwendige Verschwinden der Nennerdeterminante gema (112) eine lineare 
Abhangigkeit dreier Konstanten (Beispielsweise q2,, q3n und q,,) von der vierten, willkiirlich wahl- 
baren Konstanten (z. B. q;,,) in der Form 

qqkn = Yin Ove (113) 
zur Folge, wobei o{% eine aus den Koeffizienten des Gleichungssystems (111) bestehende Konstante 
bedeutet, so da anstelle der bisher zur Verfiigung stehenden vier Konstanten q,,, nunmehr nur noch 
eine willkiirlich wahlbare Konstante verbleibt. Da in dem von z abhangigen Teil der Loveschen 
Funktion noch zwei willkiirliche Konstanten (A, und B,) enthalten sind, diese jedoch letztlich 
multiplikativ mit der Konstanten q,, zusammenhangen, verbleiben effektiv in der Loveschen Funk- 
tion (109) je Reihenglied nur zwei willkiirliche Konstanten, und man ist somit nur in der Lage, zwei 
der insgesamt vier Spannungsrandbedingungen hinsichtlich der Stirnflachen zu befriedigen. Es 
mu8 daher nach einer weiteren Loveschen Funktion L(*)(r, z) gesucht werden, die nun auch noch die 
restlichen willkiirlichen Integrationskonstanten fee 

Bevor jedoch dazu iibergegangen wird, seien zunadchst noch die der Loveschen Funktion (109) 
gemiB (69) mit Beachtung von (113) entsprechenden Axial- und Schubspannungen niederge- 
schrieben. Es folgen schlieBlich mit den Abkiirzungen 


2( ri BO) v) {Un iJon “ie 2 Cane v) on Oy [Ln r Non ale 2 (2— 9) Non] | 
ee O13 we Jon 2 Ove Tis Non} = = fal? ) ° 


9 
a 
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E fn Dips ae ey / ee (n) N, ps. OY il mie No t (114) 
ae) (ee a) {unt Jon 2 (1 v) Jon i OD [Mn T On ( a On] 
== 2 dy3 Ww Jon ai 2 O12? we Non} Sept ay 
Jin Ap, — HG > Jin Be — Py it J 
die Spannungen 
AY, (r, 2) — 3 Jule) Othe etme — ab? emt), 
ae (115) 


AN 2) = SF galt) (ote etm + nem tn) 
Aus (115) ist die bereits erwahnte Schwierigkeit hinsichtlich der Konstanten sofort zu erkennen. 
Sie sind fiir die Befriedigung der vier Randbedingungen (99) nicht hinreichend. Nachfolgend wird 
daher eine weitere Lovesche Funktion ermittelt. 


vy) Der Randspannungszustand 4. Analog den Ausfithrungen in Absatz a) betrachten 
wir nun noch einmal eine aus modifizierten Besselschen Funktionen und trigonometrischen 
Funktionen bestehende Produktlésung. Entsprechend (89) ist also mit den willkiirlichen Kon- 
stanten Omn (m = 1, 2, 3, 4), A,, B, und (,, sofern man noch zur Abkiirzung I(@,,r) = Ip,(r) und 
Ko(U@nr) = Kon schreibt, 
Lr, z) = Dy be [Qn Lon + Q2n Kon] + @3n Lon + Q4n Kon} (An sim Un 2 + B, cos ft, 2) (116) 


An 


ebenfalls eine Produktlésung der Gleichung (68), und fiir die Radial- und Schubspannungen er- 
geben sich schlieBlich mit (116) gemaB (69) die Werte 


a E = j = ; 
a@), = TF HII {in [Hn t Lon + (1— 2) Ion] + v2n [Ht Kon + (L— 2 ¥) Kon] 
+2 03n M2 Tn + 2 O1n B? Kon} Hn (Ap 608 fin 2— B, sin fin 2) , 
(117) 
Bee E = <3 j 2 : 
7h) - 2 d. + ”) (Pe {Oin [Un r Lon aie 2 (1 ae v) Ton] = 02n [Un r Kon ++ 2 (1 = v) Kon] 


+ 2@3n (2 lon + 204m U4 Kony fn (A, sin i, 2 + B, cos Bee 


Die Befriedigung der Mantelflachen-Randbedingungen (100) ist mit (117) auf jeden Fall dann ge- 
wahrleistet, wenn fiir jedes Glied (n) der in (117) auftretenden Reihenentwicklungen folgende 
Beziehungen bestehen: 


{O1n [Un i Ton ale (1 2 v) Ton] a8 O2n [Mn r Kon SEF (1 — 2») Kon] ie 

His O3n Wed Ton ae 2 O4n pe on irers — 0 ; 
{O1n [Ln r Ton ai (1 ae v) Ion] == 02n [Mn hi Kee + (1 — 2) Kon] = 
at 2 Q3n Me Ton + 2 O4n Th Koatrers tl) 
{O1n [Mn Yr Ton + 2 (1 — v) Ibn] SP O2n (i; Lie Kon — 2 (1 —?) Kon] 
Se 2 O3n ee Ton ap O4n We Beyer — 0 ° 
{Qin [Hut Lon +2 (1—») Ton] + 020 [Hn ? Kon + 2(1—v) Kon} 

=. 2 O3n pee Le mitts 2s O4n pa Cee, == 0 3 
Aus diesem linearen homogenen Gleichungssystem hinsichtlich der Konstanten Oin (6 =, 2534) 


kénnen diese nur dann als verschieden von Null hervorgehen, sofern die Nennerdeterminante 
des Systems (118) verschwindet. Also ist 


(118) 


A(u) = 0 (119) 


die Eigenwertgleichung zur Bestimmung der ff, (n = 1, 2, 3,.. .), und das Verschwinden der 


Nennerdeterminante bedingt nun wieder eine lineare Abhangigkeit dreier Konstanten von der_ || 


vierten willkiirlich wahlbaren Konstanten beispielsweise in der Form 


s(n) 


Okn = Qin OV « (120) 
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Die der Loveschen Funktion (116) gemaB (69) entsprechenden Axial- und Schubspannungen 
erhalt man dann schlieBlich unter Beachtung von (120) mit den Abkiirzungen 


2(1 + v) Soo r Ion + 2 (2—») Ion + 8? [Unt Kon + 2 (2— ») Kon | 
= 2 Op ee Lon SF 2 OP a2 Kon = f,(r ge 


FOE) ay (lot Hon + 21>) Lon + HP int Kon +2 (1—>) Koo] (121) 
+ 26) G2 Ton + 2 SP W? Kon = Balr) 5 
Din dn =#HTO, On Be SHS! 
nach einiger Rechnung zu 
o@ (r, z) = PD (r) (#22 cos L, z— x$,” sin Up 2) 


(122) 


oo 


t(r, 2) = 28 n(T) (#92 sin UW, z + HP, cos Up 2) . 


6) Ermittlung der ere x aus den Randbedingungen 99). Uberlagern wir 
die beiden Loveschen Funktionen L(3) und L(4), so sind in den zugehérigen Beziehungen fiir die 
Axial- und Schubspannungen 


pie = BD, +O, = S fale) Oc? ote — BP etn 


Ae flr), ( #F COS Un 2 — *F 2 sin My 2) | (123) 


ft > > gn(r) (Lee? etn ® 1 BD etn *) 


i = 
+ gn(r) (#? sin Un 2 + 2B cos My 2) J 
geniigend willkiirliche Konstanten x%°, x, #39 und x$;,° enthalten, um die Spannungs- 
randbedingungen an den Stirnflachen zu Demidieen Man erhalt zunachst aus (123) fiir die Rand- 
spannungen 


a) an der Stirnflache z — 0 


Ginn (ey 0) — >; be es 1 fa(r) a xe, s) f,(r) ; 


roi (124) 
TAF, 0) = gy) [Gr,$ ae oR | gn(r ) - = xB, u &n(T) ’ 
n=1 
b) an der Stirnflache z = 1 
Crea l,l) = » [oh etn? — yh? een! f(r) + [KG cos Wn, 1— 2H, sin Lp T] fal) : 
oe (125) 


Teale) = SS [xP etm! + xP eH] gy(r) + [BP sin fin E+ 2? 008 Hin U] Bul) 


und damit ergibt sich, sofern man noch 


ch? = oe “es See HP? 3 ce s) = x}: s) : pee avi Sp ae x” ’ of: $8) _ = x} s) (126) 
setzt, aus den Randbedingungen (99) folgendes Gleichungssystem: 

& cn” fal) + Cfn” far) = — OO FP(r, 9) , 

a1 

>, cf, s) r) + ¢§ cy, DF 2nr(r ) ee sar, 0) A 

n—1 Bul ) i (127) 


[e$ Cof pin E+ oS,” Sit in UT] fa(r) + [6? cos fn 1— ef", sin i, T] f(r) =—54,,(r 1), 


[cf Sint pn E+ cf? Cof un UH] ga(r) + [ets sin i, 1+ Cf? cos Hal] Bn(r) =— TYS?)(r, I). 


9* 


“a 
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Man bildet zunachst aus den beiden ersten Gleichungen von (127) 


S eh? fale) =— [BQN 0) +S A? ale] =— RHC), 


n=] Hea) (128) 
s cf s) &n(r N= — [rl tX(r, 0) ae os Gee gar) | Lo eee RY: OG ) 
i 


und bestimmt nun die c(®,” und c$*, wie bei der Entwicklung einer Fourierreihe aus der Be- 
dingung, daB die Integrale der Fehlerquadrate minimal sein sollen. Man fordert somit 


Ta oe es ; : OF > 
Bib? = [URE %(r) + 3 oti? falr) Pdr > Min , tee leer 
‘ also (129) 
(k, s) " ¢k, 8) (k, s) 2 o une “ =0 | 
Fy ee ee | [Ry : (r) a5 2 C2; n &n(r )/? dr— Min 9 OcSk s) ache : i vp 
Lai n=1 “1 pees «2+ 00 


und dies entspricht Rare zweimal unendlich vielen Gleichungen: 


3 cf i? J fal r) dr + PRE) HO dr = 0|;21...005 
Abe (130) 


2+ 9 


Dy cf J'salr ) g(r) dr + eG g(r) dr == 


n=" ia Ts 


bzw., wenn man hierin noch die Funktionen R{**(r) und R¥(r) gemaB (128) durch die urspriing - 


lichen GréBen ausdrickt, 


[es i? J IalVS r) dr + 6 J Int) r) dr] + foun 0) f(r) dr = 0]: -1...0> | 


in 


=| 


on 


2 lee cal Sn(r) gi(r) dr + 0, iE u(r) gi(r) dr] + i AL F2(r, 0) gir) 
Aus diesen zweimal unendlich vielen feos Cine fiir die viermal unendlich vielen Unbe- 
kannten cf”, c,°, c,° und ¢cf,° (n = 1, 2...) kénnen die Unbekannten c{*,? und cf? (n = 1,2...) 
durch die Wabeereica Cf und 77 (in — 1 2...) ausgedriickt und dhe in den beiden letzten 
Gleichungen von (127) ersetzt werden, so daB in diesen nurnoch die ci,” und 672” (n = 1, 2754 
verbleiben. Auch hier ist man dann in entsprechender Weise mittels der Fehlerquadratmethode 
in der Lage, die ci’? und c$,° (n = 1, 2,...) zu ermitteln. Sind diese bekannt, so kénnen nun- 
mehr riickwartig die Keustanten ce und oh, © aus (131) berechnet werden, und mit diesen lassen 
sich gema{s (126) die urspriinglichen Konstanten x%,, xP, x9 und x$,° (n = 1,2...) be- 
stimmen. Die Loveschen Funktionen Li) und Li), d.h. der Rend gsrecee eed zur Er- 
reichung randspannungsfreier Stirnflachen sind damit vollstandig bestimmt, wenn wir noch 
samtliche Teilspannungszustande, d. h. samtliche Loveschen Funktionen L;,, aufsummieren: 


= 2, 2b En) (132) 


IV. Ebene Probleme. 1. Allgemeine Lésung fiir das ebene achsensymmetrische 
Temperaturfeld. Die entsprechenden Beziehungen des Abschnittes II vereinfachen sich er- 
heblich, wenn das Temperaturfeld unabhangig von z wird, was beispielsweise im mittleren Bereich 
eines sehr langen Rohres mit von z unabhingigen thermischen Randwirkungen ¢,(t) bzw. 9;(t) 
der Fallist. Wir erhalten hierfiir mit J,, = 0 (s = 2,3. ..) und dem Grenziibergang Ay; > 0 schlieB- 
lich mit den Abkiirzungen 
ni) (em +) — es (ain +B) 


r; 


Fier a 


(133) 


r, c,d c, dy 
an ol, Sas : 


YT; ify 


fiir das Temperaturfeld 


Ak a6 : 3 Ta “a NG ay: 
Hr, t) = 2) Eu | [reo rontre alee dee (134) 


a? [r2(@2,, + Of Die 
0 
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wobei die Werte A,, = @, durch (41) definiert sind. Analog (5) bis (7) geniigt diese Lésung fol- 
genden Rand- und Anfangsbedingungen: 


Cy Bra, t) + cy es iets Palt) . d, Hr;, t) + d, ee gilt) , Hr, 0) =0. (135) 
a,t i, t 

2, Wdarmespannungen im mittleren Bereich eines sehr langen Rohres infolge 
eines ebenen Temperaturfeldes. Samtliche SpannungsgréBen sowie die Radialverschiebungenv 
_ (Abb. 2) miissen in diesem Falle von z unabhangig sein. Hieraus folgert man schlieBlich aus den 
allgemeinen Beziehungen (59) bis (63) 

1) das Verschwinden der Schubspannungen 1,, , 

2) ein Ebenbleiben der Rohrquerschnitte z = konst. wahrend der Deformation, d.h. mit den 
zunachst willkiirlichen Zeitfunktionen A,(¢) und A,(t) fiir die Axialverschiebungen 


w(z, t) = A,(t) z+ A,(t) , (136) 


3) aus der Verschiebungsgleichung (62) folgende Differentialgleichung fiir die Radialverschie- 
bungen v(r, t) 


Onl ere 1+» of 
ro 3 a) 0 earhe y a) 
Thre Lésung lautet mit den zunachst willkiirlichen Zeitfunktionen C,(t) und C,(0) 
iep al Gt C,(t) 
vt) papas [Addr + Br 4 SO, (138) 
und damit ergeben sich nach (61) fiir die Spannungen, sofern man noch (136) beachtet, 
E Gi CG (1 + v) (1 —27) oa, ; 
et) ary a=al(s +14) 3 (L 29) hy Tae Dens 
E Cc C I+ 1—2 1 
B(r, t) = (Gee 9) ( = +y 4,) | 3 (1— 2 ») Ce ane : ee (6 = [ p ar) (139) 
E esr) Lire 2B) ig 
cA t) => al ae ”) al ero ») > C, (1 v) A, ae 1 er Xt al . 
In (139) werden zunachst die beiden willkiirlichen Funktionen C,(t) und 4,(¢) gemaB 
= sty reg (140) 


zu einer neuen willkirlichen Zeitfunktion B(t) zusammengefaBt, und daraufhin aus der ersten 
Gleichung von (139) mittels der fiir die Radialspannungen o, zu fordernden Randbedingungen 
o,(r;,t) = 9, o,(ra, t) = 0 die Funktionen B(t) und C,(t) berechnet. Es folgt schlieBlich fiir die 
Spannungen, sofern man beachtet, daB man in der Beziehung fiir o, die Zeitfunktion C,(t) gemah 
(140) durch B(t) und A4,(t) ausdriicken kann, 


. (141) 


r 
if a 


GF, t) = ole noe ((e) HE 7 [ oes L ae | aoe 


vi 7 


2 {hl = 9 
Oi) 5, a | ro dr—-—* A,() |, (142) 


wobei zur Abkiirzung FE a,/(1— v) =o gesetzt ist. Die bisher noch frei wahlbare Zeitfunktion 
A,(t) beriicksichtigt die Lagerungsverhaltnisse des Rohres an seinen Enden. Bemerkenswert ist, dah 
sich ihr EinfluB allein auf die Axialspannungen beschrankt. 

a) Der ebene Verzerrungszustand. Liegt das Rohr zwischen zwei in z-Richtung un- 
nachgiebigen Wanden, so ist ¢, = dw/dz = 0, wofiir A,(t) = 0 erforderlich ist. Man erhalt dann 
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aus (142) 
"a | 


2” 
0,(r, t) = — o)| 0-3 | r } dr 


(143) | 


d) Der Fall freier Zylinderenden. Ist die freie Dehnbarkeit des Rohres in Richtung der 
Rohrachse gewahrleistet, d.h. wirken keine duBeren Krafte in axialer Richtung auf die Rohr- 
enden, so erfordert das Gleichgewicht der zur z-Achse parallelen Krafte fiir ein an der beliebigen 
Stelle z abgeschnittene Zylinderstiick, daB die in der Schnittflache z wirkenden Axialspannungen 0, 


fir sich im Gleichgewicht sein miissen. Aus i 0; dF = 0 erhalt man schlieBlich mit (142) | 


die Zeitfunktion A,(t), womit letztlich fiir die Achsialspannungen 


i 
v 


o,(r, t) = — 0 — | rédr (144) 
hervorgeht. 

Es sei noch vermerkt, da die Formeln (141) fiir o, und o, mit o) = Ea, die entsprechenden 
Beziehungen fiir den Fall des ebenen Spannungszustandes darstellen. 


V. Niaherungslésungen fiir das ebene Problem bei diinnen Rohren und nicht zu schnellen Anheiz- 
vorgangen 0:;(t) = 3;, (1 — e— ’i') des Innenmantels. 1. Das ebene Temperaturfeld bei diinnen | 
Rohren (Abb.5). Von der Lésung (134) ausgehend erhalten wir fiir den Fall groBer Kriimmungs- | 
halbmesser mit Verwendung der asymptotischen Entwicklungen fiir die Zylinderfunktionen @o;(r) 


Abb. 5. Abb. 6. 


das ebene Temperaturfeld in diinnen Rohren, das mit der entsprechenden Lésung fiir eine Platte 
ibereinstimmt. Es ergibt sich schlieBlich mit den Abkiirzungen 


Ww 2 
Pily) = d, sin wy, y — dy ox COS Wh Y » D, = _ 2 O% 


. er Aus 
2|(q2 + d2 pi— Ay, 
( BE EOD EP dy |p (145) 
(y, t) = pilt) [ce (h — y) + 2] — Galt) [— diy + 4,] 
™Y> c, d,h + c,d, —c, d, 
fiir das Temperaturfeld 
ee Die eI a(k3 i —1) 
Blyst) = & Di Pay) | [aly t) Pay) € ) dy dr, (146) 
= 00 
wobei die Eigenwerte @, (k = 1, 2,...) aus der transzendenten Gleichung 
5 Be c,d, — cy dy 
tgoh=o aE ET (147) 


hervorgehen. Mit (146) werden analog (135) die Rand- und Anfangsbedingungen 


oo oo 
OO) ees = gal), 4 1) +S] = gle), —(y,0) =0 (148) 


‘ 
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befriedigt. Wir untersuchen nachfolgend Spezialfalle, die wir fiir einen Anheizvorgang 
Bi(t) = 0; 5(1 — e— "i*) (149) 

auswerten, Hierbei bedeuten y; ein willkiirlich wahlbarer, die Anheizgeschwindigkeit charakte- 
risicrender Koeffizient (Abb. 6) und @;, die stationire Endtemperatur des Innenmantels. 

a) Isolierter AuSenmantel. Die Isolation des AuBenmantels verhindert jeglichen Warme- 
austritt, und die Randbedingungen lauten daher 
a0) 
oy a t 
Hier kann also aus (146) die entsprechende Lésung dadurch gewonnen werden, da man c, = 0, 


¢, = d, = 1, dy = 0, pa(t) = 0 und ¢,(t) = 0,(t) setzt. Aus der Kigenwertgleichung (147) folgen in 
diesem Falle 


G(T t) = 0, dh. 0, 90,1) =8,(t), dy, 0)=0. (150) 


2k—1 2 
OF ——————— 


Te ae 


und man erhalt schlieBlich, nachdem man die Integration hinsichtlich y ausgefiihrt hat, 


Ss t Meal eNe 
9 — % 28 2k—l1 =a eer (t —z) 
Wo pip Servant stag fame Peal te 
e 6 
und speziell fiir den Temperaturverlauf am AuBenmantel 
oo t 2k—l «\* 
x : (¢—1) : 
Balt) = Hh, t) = (a ie (— 1)! (2k —1) 9,(z) e ( 2 a) dr. (152) 


() 


Wertet man in (152) das Zeitintegral fiir das durch (149) gegebene Anheizgesetz aus, so erhalt man 


schlieBlich mit Beachtung von 
2 y\k—1 . 
>§ a eke =F 


k=1 
sofern man noch zur Abkiirzung 
am \2 
(555) inl 

setzt, 

‘ S41) | er I)? gk Irs] 

0,(t) — Ne Ome ( ) B ( ) Vv; e 
a gol 2K —I1 [P(ZK—IP—% SO are 


und fiir die Differenz Ai(t) = 9;(t)— 3,(t), die sich in der nachfolgenden Naherungslésung fiir 
die extremalen Spannungen im wesentlichen als proportional zu den Spannungen erweisen wird, 
ee Sy ea 
== a a a vb te = a ~ = 
AB(t) = Bit) — Balt) = Bie mS > sr .—) (153) 


k=] (2 3 


(Abb. 7). Bricht man die Reihenentwicklung (153) nach dem ersten Gliede ab, so erhalt man 
fiir die Temperaturdifferenz eine Naherungslosung 


A eee gt 


AY(t) a Se — 
Ge 
Vi 


die, wie eine entsprechende Restgliedabschatzung ergab, fiir alle Zeiten ¢, 
fiir die , 


(154) 


Abb. 7. 


prep =08 (155) 


ist, bis auf einen Fehler von héchstens 7,5 °% die tatsachlichen Temperaturdifferenzen A(t) 
approximiert. Fiir Stahlrohre hat man /[sec™] = (7/2 ah)? = 0,304/h?, wobei h die in cm einzu- 
setzende Rohrwandstarke bedeutet, so das damit aus (155) ft’ = 0,304t'/h? = 0,8, d. h. t! = 2,63 h® 
folgt. Man ersieht hieraus, daB der vernachlassigte Rest um so schneller belanglos wird, je kleiner 
die Rohrwandstiarke h ist. Betrachtet man beispielsweise ein Stahlrohr mit einer Wandstarke von 
h — 3m, so wird das vernachlissigte Restglied der Reihenentwicklung bereits nach t’ = 9-2,63 
= 23,7 sec kleiner sein als 7,5°%. 
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Aus einer Extremalrechnung fiir (154) erhalten wir die groéBte Temperaturdifferenz 
1 


Dis _ A 


B 


gine 
Op; v7 


max Ad cS. 4 ieee (156) | 


zur Zeit 


und indem wir ft, > ft’ = 0,8 setzen, erhalten wir schlieBlich B/v; = 1 als Kriterium fiir diejenige 
Anheizgeschwindigkeit, bei der (156) mit einem Fehler von héchstens 7,5% gilt. 

Fiir die Bestimmung von »; wollen wir annehmen, daB der Anheizvorgang am Innenmantel dann 
als abgeschlossen anzusehen sein soll, wenn die stationire Endtemperatur 3;, zu 99%, erreicht 
worden ist. Demnach hat man, wenn man mit f, die bis zu diesem Zeitpunkt verstrichene (Anheiz-) 


In 100 4,6052 
Zeit bezeichnet, 9;(t,) = 0:5 (1— ei) = 0,99 8;, woraus »; =—— = —" 


; hervorgeht. 
Hiermit folgt aus (156) 


1 


Ad 4. [450052 \1—(4,6052/6 te) 
—% = || Bt, =f(6 t,) (157) 


(Abb. 8), und zwar mit einem Fehler von héchstens 7,5°% fiir alle Anheizzeiten t,, die der Beziehung 
Bt, = 4,6052 geniigen. Der Wert f(t.) ist im Bereich 0 < 6 t, < 4,6052 weitaus ungenauer, und 
nimmt fiir t,=0 den Wert 4/z an (gestrichelt), wahrend in Wirklichkeit in diesem Falle 
max A?/);,=1 sein mu. Die strichpunktierte Linie zeigt den genaueren Verlauf von 
f(B t.) in diesem Bereich. 

b) Warmeiibergang in ein am AuBenrande befindliches Medium der Tempera- 
tur Null. Entsprechend den Ausfiihrungen unter a) kommen wir auch hier zu einer Nahe- 
rungslésung, wobei nunmehr jedoch f =f(f t,, N) noch in komplizierter Form von der Nusselt- 
schen Zahl N =A//a h abhangigist. Es ergibt sich die in Abb. 9 qualitativ skizzierte Kurvenschar, 
deren Grenzkurve N = co wieder den Fall des isolierten AuBenmantels darstellt. 

2. Naherungslésung fiir die Spannungen. Nachdem wir in Ziff. 1 Naherungslésungen 
fiir A9(t) und damit auch fiir 9,(t) bei vorgegebener Innentemperaturfunktion #);,(1— e—i') be- 

stimmt haben, gehen wir mit diesen nun als vorgegeben anzusehenden Rand- 


ee =f (Pte) temperaturfunktionen auf die Lésung (134) zuriick, in der gemaf des Rand- 
x und Anfangswertproblems 
i Hr t) =O(t), Ar, ) = 9(t), —-P(r,0) =0 (158) 
i \ op de 2 Virey ad, = 0,rq, (3) == HO) sano ead 
a 8;(0) n ™ — 8,(0) In 
<< (r,t) = ae (159) 
05 i. 7 a 
reds 
eee a Ve ees ! | 1 fitz 
0 5 70 5 ros 
Gt 4 605 — 
Abb. 8. Abb. 9. 


zu setzen sind. Indem wir in (134) auf das Zeitintegral die partielle Integrationsregel a d d 
gemaf (48) und (50) beachten, daB 8 get anwenden un 


r 
a 


. 2 Dy KA) ; ’ 
JG ae mr, t) + + Dox(r) dr = nr, t) 


Cen) 
t 
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ist, so erhalten wir schlieBlich mit Riicksicht auf 8;(0) = #,(0) = 0 


r 


a Or, 2 
fees 
[ renin See i: dr dr . 


t 
oo 


2 2 Dy,,(r) 
B(r, t) = 7(r, t) Da [r? (2, 4 ®; aye | 
; 


Pan ha 


r. 
v5 


Diese Formel vereinfacht sich fiir den Fall hinreichend langsamer Anheizvorgange zu 
9,(0) In“ — (0) In 
AAT) ITSO) ee (160) 


Ing 


Sie ist nun ihrem Aufbau nach eine stationare Lésung, in der die Randwirkungen zeitlich verander- 
lich sind und mit dieser quasistationaren Lésung werden nun die in (141) bis (144) auftretenden 
Integrale ausgewertet. Es ergeben sich schlieBlich fiir den Fall diinner Rohre folgende Extremal- 
spannungen: 


a) Radialspannungen 
1 I 
Grexir(t) = — | 00 AV(t) z ; (161) 
b) Tangentialspannungen 


min p(t) = o9(t:.t) =—- 0 40(0) A a (a) | 


2 (ay | (162) 
eas i 


c) Axialspannungen im Falle des ebenen Verzerrungszustandes 


max Oa (¢) = GAT a, t) == = Oo AK(t) 


min ¢,(t) = 6.(r;, t) = — 0, |9;(t) — = (9,(t) + B.(t)) +» an) a | | 
(163) 
max 0.(t) = CALee t) =— 05 | A.(0 ry = (8,(t) Ja(t)) Ee oe =| | 


und im Falle freier Zylinderenden 


min o,(t) = o.(r;, t) = — +0 Ad(t) (2 =f = “al , 
(164) 
mare. (Yo N= 6 49() (I~ 5): 
z a\lao 2 0 a) Re 

3. Zusammenfassung. Die extremalen Tangential- und Axialdruckspannungen am Innen- 
rande sind also stets absolut gréfer als die entsprechenden (Zug-)Spannungen am Aufenrande. 
Im Falle freier Zylinderenden sind die extremalen Tangential- und Axial-Druckspannungen bzw. 
Zugspannungen gleichgroB, im Falle des ebenen Verzerrungszustandes hingegen sind die Axial- 
spannungen gréBer. Die Radialdruckspannungen sind in sémtlichen Fallen im Verhaltnis zu den 
Tangential- und Axialspannungen klein und daher fiir die Dimensionierung eines Rohres un- 
wichtig. ‘ 

Weiterhin erkennen wir an Hand der Naherungsformeln (162) und (164) im Zusammenhang mit 
Abb.8 und 9, daB grundsatzlich wahrend des Anheizvorganges héherer Spannungen auftreten als im 
stationdren Endzustand entsprechend der Tatsache, daB die Tempeaturdifferenzenkurve wahrend 
des Anheizvorganges einen Extremwert hesitzt (Abb. 7). Dieser Extremwert ist auch im allge- 
meinen Falle des Warmeiiberganges in ein Medium der Temperatur Null auf jeden Fall grofer als 
die stationare Temperaturdifferenz A/),, jedoch erkennt man aus Abb. 9, daf fiir den Fall schlechter 
Warmeleitfahigkeit (N= A/a h— 0) der stationére Temperaturdifferenzen-Endwert nicht viel 
kleiner ist als diejenige Temperaturdifferenz, die sich bei kurzfristigem Anheizen ergibt. Daraus 
kénnen wir vor allem fiir nicht zu schnelle Anheizvorginge folgern, daB die Dimensionierung von 
Rohren mit sehr schlechter Warmeleitfahigkeit (z. B. Betonrohre) sofern sie aufen nicht isoliert 
sind, mit geniigender Genauigkeit nach den stationaren Endtemperaturverhaltnissen erfolgen 
kann. Bei guten Warmeleitern (Stahlrohre) ist dieses vor allem bei kurzfristigem Anheizen un- 
gulissig, da hierbei die extremalen Temperaturdifferenzen und damit die gréBten Warme- 
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spannungen ein Vielfaches der stationiiren Endwerte ausmachen konnen. Im Falle eines auBen 
isolierten Rohres ist 49, sogar Null, so daB im Endzustand hier iberhaupt keine Warmespan- | 


nungen mehr herrschen. 
Wir betrachten abschlieBend ein Beispiel: Wie schnell darf ein 3 cm starkes, auBen isoliertes 


| 


Stahlrohr mit freien Enden vom Radius R= 10cm auf die Endtemperatur #;, = 200° gebracht 


werden, damit die zulassige Héchstspannung o,,, = 0,5 t/em? nicht tiberschritten wird? Gegeben 
sei E = 2,1 ° 10° kg/em?, v= 0,3; a, = 0,000 012. 


Die gréBten Spannungen sind die Tangential- bzw. Axialdruckspannungen am inneren Rande. | 


Sie betragen mit o, = Eo,/(1— v) fiir A0(t) = max Ad 


pa (1 |. : Ri) max Ay; 


On eas SS lo» mas — 2(1 =e v) 


und indem man [Grea == Walmae <6, setzt findet man 


ae a ena : 1 = max Ad = 0.1085 
max A? < ne ch ZAG neo ER oe: (PE )ja0, . 


“Man zieht nun in Abb. 8 eine Parallele zur Abszissenachse im Abstande 0,1085. Diese schneidet 
die Temperaturdifferenzenkurve f(f t.) ungefahr an der Stelle ft, = 30. Fir ein 3 em starkes 


Stahlrohr erhalt man B = (2/2 ah )? = 0,304/h? = 0,0338 und damit folgt schlieBlich aus f t, = 30 | 


30 


is 3 887,5 sek = 14,8 Minuten als kiirzeste Anheizzeit. 


VI. Achsensymmetrische Probleme bei diinnen Rohren (Kreiszylinderschalen). Nachdem man ° 
mittels der asymptotischen Entwicklungen der Zylinderfunktionen und der modifizierten Bessel- 
schen Funktionen die Lésung fiir das Temperaturfeld entsprechend vereinfacht hat, gelingt die | 


Spannungsermittlung in derselben Weise wie beim achsensymmetrischen Problem an dicken Rohren, 
Einem durch ein thermisch-elastisches Verschiebungspotential definierten Spannungszustand 
miissen wiederum Randspannungszustande zur Erreichung oberflichenspannungsfreier Rander 


iberlagert werden. Hierbei wird die Spannungsfreiheit an den Mantelflachen bereits mittels eines — 


zusatzlichen Membranzustandes erreicht, wahrend man das Problem der Oberflachenspannungs- 
freiheit an den Stirnflachen mit Riicksicht auf die Achsensymmetrie naherungsweise auf ein Rand- 
stérungsproblem eines nachgiebig gebetteten Tragers reduzieren kann. 


(Eingegangen am 28. Mai 1955.) 


Anschrift des Verfassers: Dr.-Ing. Rudolf Trostel, Technische Universitat Berlin, 
Lehrstuhl f. Mechanik (Prof. Dr.-Ing. Szabd) Berlin-Charlottenburg, HardenbergstraBe. 
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Bestimmung von Rahmenschwingungen mit Hilfe des Matrizenkalkiils 


Von H. Fuhrke 


1. Einleitung. In einer friiheren Arbeit: ,,Bestimmung von Balkenschwingungen mit Hilfe 

des Matrizenkalkiils“‘1 wurde die Matrizengleichung fiir die Querschwingungen des viel-feldrigen Bal- 
kens hergeleitet, wobei elastische oder starre Zwischenlager vorhanden sein durften. Die Methode 
soll nun ausgedehnt werden auf die Bestimmung der Eigenfrequenzen ebener Rahmen. Dort tritt 
eine Koppelung von Biege- und Lingsschwingungen oder von Biege- und Torsionsschwingungen 
auf, je nachdem ob der Rahmen in seiner Ebene oder quer zu seiner Ebene schwingt. 
___ Die Matrizen fiir die Biegeschwingungen kénnen aus der fritheren Arbeit iibernommen werden, 
die fiir Langs- und Torsionsschwingungen leiten wir hier ab, AnschlieBend formulieren wir die 
Ubergangsbedingungen an den Rahmenecken und erhalten schlieSlich die Eigenwertgleichungen 
fiir offene und geschlossene Rahmenwerke. 


2. Matrizengleichung fiir Langsschwingungen. Die Differentialgleichung fiir Stablangsschwin- 
gungen ist von der zweiten Ordnung. Damit wird auch die Matrix fiir Langsschwingungen sehr ein- 
fach: sie ist ebenfalls von zweiter Ordnung. Da es fiir die Herleitung der Differentialgleichung 
keinen wesentlichen Mehraufwand bedeutet, wollen wir zundchst auch innere und auBere Dampfung 
beriicksichtigen. Auf die numerische Anwendung verzichten wir allerdings. 

An einem Stabelement (Abb. 1) von der Lange dx wirke die AuBere, geschwindigkeitspropor- 
tionale Dampfungskraft k,udx; das Newtonsche Grund- 


: Feld i 
gesetz hei®t dann Nie MN, 
anes a) “LI 
. . . . . @ e: ur 
Die an einem Einheitsraumelement wirkende Spannung o sme 


ruft einmal eine Dehnung ¢ =u’ hervor; zum anderen 
mu sie aber auch die innere Dampfung tiberwinden, 
die der zeitlichen Dehnungsdnderung stets entgegen gerichtet ist, und deren GréBe wir pro- 
portional é annehmen wollen: g = Ee + ré, oder mit r/E = k, 

N , 5f 
Fir harmonische Schwingungen erhalten wir aus (la) und (1b) die Differentialgleichung fiir Langs- 
-schwingungen 


Abb. 1. Dehnstab. 


a” +(7) =0), (1e) 
wobei ohne Dampfung 
| 2 = Sep (2) 
und mit Dampfung 
oo (1 er et) 
2 ao BP (2a) 


™ EF (1 + io kp) 
‘ist. Die Dampfung kann also einfach durch eine ,,komplexe Masse“ bzw. einen ,,komplexen 
E-Modul“ erfaBt werden?, so daB sie in der weiteren formalen Behandlung der Differential- 
gleichung (lc) gar nicht mehr in Erscheinung tritt. 

Der Ansatz 


(= sin) +b cos A 
geniigt der Differentialgleichung (1c). Nach (1b) ist N = EFu’ (worin E bei Dampfung ,,komplex 
und frequenzabhangig“ ist). Fiir x = 0 (Abb. 1) ist u; = 6 und N; = EF (A/l) a. 
1 H. Fuhrke, Ing.-Arch. 23 (1955) S. 329. Gleichungen aus jener Arbeit werden hier in der Form (I,.. .) 


zitiert. 
2 K. W. Wagner, Schwingungen und Wellen, S. 262 und 370, Wiesbaden 1947. 
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Damit wird 


u(x) = = N; sind + u; cos 


und fiir x = 1 erhalten wir mit 


N = i Es (3) | 
l L EI 
und 
9, == (4) 
13 EF; 
die Matrizengleichung 
Raters 
ie [ cos 7 ee sin, u | | 
| =], Pe (5) 
N === sind  cosd | ENG 
I deed L Q eles eee 
oder 
ina] = [ci] = Leelee 


Die Komponente N des Vektors y ist auf Grund der Definition (3) dimensionsgleich mit u, wobei 
die vom jeweiligen Feld i unabhangigen BezugsgréBen 1 und EI aus Teil I, Ziff. 2 ubernommen 
wurden, da wir spater gekoppelte Langs- und Querschwingungen behandeln wollen. 


Sonderfalle sind: 


a) Furl; 0 oder EF;— co | b) Firy;—0 
mit 4; 1;—> m; | (oder m = 0, Statik) 
folety, == 05 (6,==0 | wird A; = 0 
und wir erhalten die Massen-Matrix und damit die Elastizitats-Matrix 
: : 
M, =| ane a 6) G&= lr a (7) | 
—o?m,B/EI 1 | | 


c) Kombinieren wir a) und b), so erhalten wir eine der Matrix J (I, 31) entsprechende Matrix 
3 =M- € fir ein elastisches, masseloses Feld mit Einzelmasse m am linken Ende: 


1 Oi A B 0 0 | 
= — = G; == (w ti F (7a) 
eens is EI 
——— 00 Se 1— a? UU Oe 1 0; 
EI EI 


Bi =—— 


zy 


Das in Teil T Gesagte ttber Abschnittsmatrix, Randbedingungen, Kigenwertgleichung und Eigen- 
funktionen gilt analog auch hier. Die Eigenwertgleichungen sind jetzt sehr einfach. Fiir jeden 


Rand gibt es nur zwei Méglichkeiten: Festgehalten: u = 0, oder frei: N=0. 


Fir den frei-festen, homogenen Stab z. B. erhalten wir aus (5): 0 = N, cosd. Die Matrix 9 fiir 
den homogenen Stab, oder allgemeiner die aus i Feldmatrizen nach (I, 29) entstandene Matrix Q fiir - 
den beliebig zusammengesetzten Stab enthalt bereits die Loésungen fiir alle méglichen Rand- | 
bedingungen (Abb. 2). 


3. Matrizengleichung fiir Torsionsschwingungen. Wie bei den Lingsschwingungen kann auch | 
bei Torsionsschwingungen die auBere und innere Dampfung durch eine komplexe ,,Masse“ oder durch 
einen komplexen ,,Gleitmodul” beriicksichtigt werden, wobei allerdings der Dampfungsfaktor nicht | 
nur material- sondern auch formabhangig ist!, Wir wollen die Matrizengleichung, in der die mes 
fung ohnehin nicht in Erscheinung tritt, fiir den Engcdam piles Torsionsstab herleiten. | 

Die Drehmasse des Stabelements dx (Abb. 3) ist dO= wi; dx. Die Bewegungsgleichung heibt 


dann «i, dx 7 = dMr, oder fiir harmonische Schwinguneen: | 


—wwiny = M7. 


' K. Klotier, Technische Schwingungslehre, 2. Aufl. Bd. 1, S. 237ff. Berlin, Gottingen, Heidelberg 1951. 
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Mit der Elastizita i hoes : : : : 
ign izitatsgleichung 4 = M7/GIy wird daraus die Differentialgleichung fiir Torsions- 


es ae = 0 mit! He =o hse I? 
] Ae Gift (8) 
Setzen wir noch entsprechend (I, 6) 
=i! (9a) 
Mr; = My, P/E (9b) 
i; EL 
und Sg poe = 
atGl ey) 
Ming X 
I ing Xivs Feld i Xi 
mis ( (EEN 
—— A y Stelle: i+7 tL Mp, 
27 22 x Y 
et 
Abb. 2. Abb. 3. Torsionsstab. 
so erhalten wir in Analogie zu (5) 
ers] if - iain ten. || Fig ee 
7 COSA 7 sin, | ¢ | 
= | Ic | a (11) 
Mr == Say} cosh | | Mr | 
Jjit+l ie be ji [ ji 
oder 
Sonderfalle sind: 
a) 1,0 oder GI7;— a ib) = ON 
mit Mi ini l; gs 0; ’ 
i 0 
a ip | i i 
HE Sota 1 (12) OR k He (13) 
EI 
i Vi j 0 0 
8:= 0 -G = 1 1 || =6.— a 0; —— 
3 a — wu? 0; Bs 1— ow? 6; 4; Bm Ee Ae EI . @ 
EI EI 1 4 


Fiir eine Lauferkette mit der Gesamtmatrix 2{ (Abb. 4) wiirde die Frequenzgleichung heifen? 
A 


4. Eckmatrix fiir Quer- und Langsschwingungen. Wenn Riegel und Stiele eines Rahmens 
gleichzeitig (iiber die Ecken gekoppelte) Quer- und Langsschwingungen in der Rahmenebene aus- 
fiihren kénnen, so laBt sich hierfiir die Frequenzgleichung in Form einer Matrizengleichung auf- 
stellen. Zwischen zwei an einer Ecke senkrecht aneinander anschlieBenden Schwingungsfeldern 


kann die Verbindung formal durch eine Matrix beschrieben werden. 


1 Der Parameter A ist natiirlich mit einem geeigneten Index zu versehen: /,, 1;, 4,, wenn sich mehrere 


Schwingungsarten mischen. 
2 Vel. C. B. Biezeno u. R. Grammel, Technische Dynamik 2. Aufl. Bd. 2, S. 347ff. Berlin, Gottingen, Hei- 


lelberg 1953. 


| 
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Diese Verbindung kénnen wir zunachst durch sechs Gleichungen herstellen: drei geometrische 
Beziehungen und drei Gleichgewichtsaussagen (Abb. 5): 


Noe nue Nigi= Qs 
LUE | ET Oy Qii1=—N;, (15) 
eS Wi» VA M; F 


Diese sechs Gleichungen lassen sich nun zu einer Matrizengleichung zusammenfassen mit der Eck- | 
matrix % : 


baa F 0-00) 100 a One 
| N (0 30) Oo Oy eaeeL N | 
ow Be we (16) 
y 09-0100 100 nO p 
M 10 0-0) > 0s.) Ae 
One NSO 30 OPO ia Osa 


oder 
[Pisa] = [ss] = [iets 


Die Vektoren }) setzen sich dabei aus den Vektoren y fiir Langsschwingungen und }) fiir Quer-| 
schwingungen zusammen. Dementsprechend miissen auch die anschlieBenden Feld- oder Ab- 
schnittsmatrizen zwischen den erweiterten Vektoren ) vermitteln; sie miissen also aufer einer Ma- 


} 
| 


| 
i 
| 
| 


Feld u+7 t 


Stele t+7 


Abb. 5. Rahmenecke (Biegung und Dehnung). Abb. 6. Rahmen. 


N 
N 
N 
\ 


Vdd 


trix 2 fiir Querschwingungen auch noch eine Matrix fiir Langsschwingungen 3{ enthalten. So mufi 


z. B. die Matrizengleichung zwischen ); und );_,; wenn %; , und Jt;_, die Matrizen fiir Quer- 
und Langsschwingungen im Abschnitt i — 1 sind, hei®en 


Uses ING NG 0 0 0 O° 7 ru 
N Nee ees t) 0 0 0 N 
w 0 0 Ay, A, is Ai ee | 
y “| 0 0 A,, Ay, Ay; Ag, | y j a | 
M 0 0 As, Asp Ass Az, | M 
L Q cee atl 0 Ay, Ay, Ay; ae, 0 Biren 


oder 


ip] = , (Bal [pia] 
Fiir den Rahmen der Abb. 6 erhalten wir z. B. die Matrizengleichung ' 
Pe = 85° Ba Bs Be Bip =H. - (18) 
Hierbei kénnen die Matrizen $8; wieder aus beliebig vielen Feldmatrizen durch Matrizenmultiplé- | 


kation entstanden sein. Jt und { kénnen zunichst getrennt aufgestellt werden. Sie brauchen nicht 
einmal die gleiche Anzahl von Feldern zu haben. | 


Wie in Teil I Ziff. 5 kénnen wir nun entsprechend den Randbedingungen aus der Matrix § (18) | 
die Kigenwertdeterminante gewinnen. Von den je sechs Komponenten der Rand-Vektoren sind 


XXIV. Band 1950 Fuhrke: Bestimmung von Rahmenschwingungen mit Hilfe des Matrizenkalkiils 3 


gewohnlich je drei gleich Null. Bei dem Rahmen nach Abb. 6 z. B. entsteht aus 


“as - EAL 
US ae 
0 | | 0 : Ay Hyy Hyg | 
a 9) es (19) © als Kigenwertgleichung 4A=/ H;, Hy H3,|=0, (20) 
M, | | 0 | Has Ay, F 4g | 
2 0, | | | | 0, ; also eine Determinante dritter Ordnung. 


5. AA-Matrix zur Matrix sechster Ordnung. Zu den in Ziff. 4 angegebenen Matrizen sechster 
Ordnung lassen sich die zugeordneten ,,4A-Matrizen“ aus Determinanten dritter Ordnung ent- 


sprechend I, Ziff. 6 bilden, wenn man das folgende Schema heniitzt: 
Zeilen bzw. Spalten-Kombinationen p bzw. q = 


en oer O oe let < 16 1gl = 20 
ied best Rd Tiss Wt Se 
DeDede ors Dee DDD). 
3h DrOLOe 33h OrosDe 
A. 4. 4.4, 4. 4.4. 4.4. A, Zeile bzw. Spalte. 
De Sy oa ar Sas 8h Os 
6. iy COO 6. 6.6. 6.6.6. 


(21) 


Damit erhalt man z. B. aus der Matrix ‘$ (17) die Matrix ‘$44, deren Elemente Kigenwert- 
Determinanten sind, und zwar fur die links neben der Matrix skizzierten linken und die oberhalb 


der Matrix skizzierten rechten Randbedingungen (Abb. 7). 
eee a eee eee 


uw 0 0 


Ny MX* 


Neg Neg" 


is 0 0 L «a 
4 | 
~*~ | eee eee 


Abb. 7. AA-Matrix mit den zugeordneten Randbedingungen. 


-nung) von Y. Zwischen 9(44 und der Kehrmatrix 2(-! besteht die Beziehung 


“ae a ae ee 

laa ae, wee 48, an | = Daal | ie 4a ae 

| AM Af4 Ai Als : | Eye Le Aa 

| eee 42 43 44 ie 41 31 ce 
| wobei fiir alle hier vorkommenden Matrizen die Hauptdeterminante |Q(-'| = 1 ist. 


)metrischen Balken ist ferner: A,, =(—1)?'4: Apq = AZ4 oder Y44 = Y. 


| Die Elemente AZ} der 4A-Matrix (44 in $844 sind die — in der Zeilenfolge p und Spalten- 
‘ folge q =1, 2, 5, 11 des Schemas (21) zusammengestellten — 4.4 Unterdeterminanten (dritter Ord- 


Fiir einen sym- 


in Langsrichtung starr (Abb. 8). Mit 
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Zur Matrix % (16) fiir den Eckiibergang gehoért 
Ee 0: 0 0-0) 00 1 10080 00s 20-0050 00200 
0000 000 °0°0°0) 0505.0 000s soe) | 
0.°0' 0.107 0-005 90505050, 0; 0c 0R0K0 RRO 0 ae a | 
0°.0.0°.0 » 0°00) 20) 0209 02012 -19050°07 05.00; Bae 
0:00 705 1900) 6705050 VOTO SO TOR 020 Sem OR 0 ae) 
00.0505 F0"T 0" 0) 070050 BO 0n0 v0 OKO 
—L000 000 000 00° - 000.0 0:0°0, 30 
0:00:04 6 010 050.007 OLOim0F0 0200 12050" er 
0:0 0°00. 0°00. 05050) “1209.0: OF OL 0 O20 nO Om] 
: 00:00: 0.0. 6.0 00 OL > 020:0,0 42070025 0s 
ey 3) 
0-0'0"0: 0: 020" f0> 1L0> 400.05 FO) 05050 OSOnO m0 
0':0:0°0° -0.0°0" 20.0 100.0" 20101050 ©0000 Os 
0 0.071 0.00) 0.070) O07 0.050505 0F0 Oe, 
0.0700 0/0°0, 9000) 0-05 801050205 05050 
0:0 0.0) 0.070" 10°00 P2020 0200 OSU sino am 
0.0 070. < <070.0' .070 0) 07055 0202051 0200 neo | 
00: 0°0" 00.0 17050" 00% 0010700 0 0 | 
ie 0: 10.0'0.0,0— (05050> 5010) 0505050 000m maG | 
0.0° 1,0. 0:0: 0~ 0020, ~ 0.0") 70 0050s 00st 0 | 
; 00 00 70 00 0 On0 00) 2 0)17020e OL0nO msm | 
Als Anwendungsbeispiel behandeln wir den unsymmetrischen Rahmen mit Stielen und Riegeln 


D = Bo BB. GB 


ist H24, = 0 (24) 
die Eigenwertgleichung. Dabei wird mit m = (m, + fu. 1, + mg) B/EL und 
10 } a | 0 | ed? 200) | 
| 3 0) | 0 0 
= 1 | —w?m 1 | Ome 
ot a | . ’ 8= | Bees | ? i= | | (25) 
L 0 C | Ie 0 B | U 0 2 at 
die 5. Zeile von die 16. Spalte von die 16. Spalte von 
(Bs 5) 44 (Bs 5 Bi) 44: = (5 Pie 
ee —By As; | <7 
0 ; | 0 
0 | 0 
0 | —Ad 
Cr Be “As, as BY AG, 0 
Ch | BY, an ale Jase Ass 0 
0 rx0 
ee | 0 
0 Ad, 
om | hd Ax, | (26) 
Cis | peel Ee Aneri Baga) Beaver i An 0 e || 
vis lee ene m (Bo. Ana 1 As) st Be Age Be Ad. | 0 | 
| 0 
0) Oo 
0) 
0 | AS 
y | ; Af, 
Ch, | Uae ae 0 
0 | 5 0 
0 0 
0 | 
| L Ai, 
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Folglich wird die Frequenzgleichung 
Fsto = Cty Bu Ag + Ch (BY, Af + BY, Ag) + CA(BA Ad + BA Ad) 
ra Cae ee a? m (Be, At, + Br, As) + Bis Age Be Ag.) 
= Cr — wm es As 7° Bs, Ax) “7 Bs, AG, rT By, Axe) a3 Ci Buy At =0. (27) 
6. Schiefwinkeliger Anschlu8. Entsprechend (16) 1a®t sich auch eine Eck-Matrix % fiir einen 


EckanschluB aufstellen, wenn zwei Felder unter dem Winkel y aneinander anschlieBen. Aus Abb. 9 
-liest man ab 


Ujii= UU, COSY + w; sing, 
W;41= —U; sin p + w; cos @ , 
Ye Wi» 


Nea N; COs — + Q; sin > 
Q341= —N; sin yp + Q; cos LY > 
Miyi= M;. 


Die Eck-Matrix in );.;—= % - p; wird also 


| cos p 0 sin p 0 0 (| 

0 COS P () 0 ) S=rn p | 

ee. jo Say 0 COS P 0 0 0 | 
. 0 0 0 li 0 0 a) 

0 0 0 Of eal ORs 

L 0 sin @ 0 0 0 cos 1 


Die AA-Matrix hierzu wird [siehe Gl. (28) S. 34]. 
Mit m = 0 wird hieraus die Einheits-Matrix, mit g = 90° die Matrix 444 (23). 


7. Eckmatrix fiir Quer- und Torsionsschwingungen. In der gleichen Weise, wie wir in Ziff. 4 
die Quer- und Langsschwingungen behandelt haben, wollen wir hier auch verfahren fiir den Fall, 
da8 zwei senkrecht aneinander anschlieBende Schwingungsfelder oder Feldketten nicht in ihrer 
Ebene, sondern senkrecht zu ihrer Ebene schwingen (Abb. 10). Wir kénnen wieder sechs Aussagen 
fiir das Eckfeld i machen: 


Xa1— Yi Miii= My; 
UF BAG a5 Ors Myivi1= — M;, (29) 
Ws Wj; Qi+ = OF 
Hieraus erhalten wir die folgende Matrizengleichung mit der Eckmatrix : 
fy 3 [ 0 0 0 1 0) 01 ey as 
Mr oa 0 0 0 d 0 | | Mr 
lw | | 0 0 1 0 ) 0. w | 
i ae 0 0 0 0 Mat ere: ay. 
M | 0 =i 0 0 0 0 -M 
OMe oP t 0 0 0 0 0 Lite *#1RG)? <I 


oder 
[Gini] = [S] GE 
Die anschlieBenden Feld- bzw. Abschnittsmatrizen miissen jetzt die Matrizen Q( fiir Querschwin- 
-gungen und & fiir Torsionsschwingungen enthalten. 
Die Matrizengleichung zwischen g; und qg;_; hei®t dann entsprechend (17) 
Re Ore ay a 31 
gi \ a Crt = ee « (31) 


Hiermit 1aBt sich z. B. ein Rahmen nach Abb. 6 erfassen, der quer zu seiner Ebene schwingt. 


8. Trigheitsfreie Anschlufglieder. Wir wollen die in I, Ziff. 8 gewonnene Matrix © fir die 
elastische Stiitze (I, 47) noch erweitern a) fiir zusatzliche Langsschwingungen (Abb. ]la) und b) 
fiir zusitzliche Torsionsschwingungen des Riegels (Abb. 11b). Auch hier seien die Anschluf- 
glieder (Stiele) masselos. 

3 
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% ee) 
o B oY 
f a 1 
iI S& 
6 Seco SOS SDSS SOON SS OOS e 
. o 
& 
S iI 
S q a ie 
ooye coooosm o8sgoeceo SES 
Os & g ° 
, Z ° 
Abb, 8. Rahmenbcispicl. o, 
S 
& =| 
~ q a e 
S = 
eco S522 re) poooee ety 
7) 3 o 
° | 
S; & 
cooo Soo foo scooscsse nooo 
bd } 
x ecoeoo esceosces SosoeoscoH SocSoSo 
Z osccoo eoceoseoso scoeoscHS coos 
Abb. 9, Ubergangsbedingungen an einem beliebigen Knick. & Se 
MOS | EON OSA ooogs 
wn 
. - o 
Feld t+7 z 
S & 
Stelle: v+7 cook oooece eogoos SOS° 
¥ 
LZ Yi 
Bay i S 
q & 
2 & S q 
OS Fo coco Oo% SHOSSS SOGS 
P 5 2 
3 n o ie 
} 
o 
Abb. 10. Rahmenecke (Biegung und Torsion). s = 
& S 
be S S ro] 
SoS ooooye woocecso — a) 
S a } & 
i} 
o 
: So 
= a 
Be s SS a 
eoge Sosceosog S4SSSS SSR 
Pa ) a 3 
3 =I 
eo 
& qe! 
a a 8S: é 
o3see SoS © -“zeccceo ogco 
“4: 5 a 3) 
Abb. 11. Mehrstielige Rahmer. 8 . 
Me Qa My, = 
- esceoo sceosco8feso scsoessos Hoos 
thy Z no DB 
YA {7-\KZZA|-) EZ : i 
(oj a 
a Qo x a [Se 
vA va seco oonococo oococo yy) 
b te 3 
Ls 
My coco omococoo coooocceco coco 
“— 
Ob 
ayy eooco nHooocooo —— KK — JK — = eoocooo 
b 
Ys <— S 
& 
L b oe oocoocoococo co fkooso —— 
2 i 
S 
= 
Abb. 12. Ubergangsbedingungen an einer Stiitze. & a B a 
a 
CW myat gH 5 oe 4 
4A b LLL ® 
Cy b Oy b rs) 
o & 
S ol oo <3 
Snoeoce coscs So Seooooo co foo 
| 5 Sues : 
oder . © Ss) 
| & As 
| OSES: a i 1) 
o 
| e& = at 
Cy | Cy 
al vy | 
Mp=1 


Abb. 13. Stiitzelastizitiiten, 
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a) Biege- und Langsschwingungen. Aus Abb. 12 lesen wir sechs geometrische Bezie- 
hungen und drei Gleichgewichtsaussagen ab: 


Ue) Wh » Ne. =N.z+Q, 
W, = Wy = — Uy, (32a) 0, =0,— Nz, (32b) 
Uy a Wb» | Me M, 42 Ma. | 


In den folgenden drei Elastizitatsaussagen fiir den Stiel bezeichnen wir die Federzahlen c 
nicht nach den Kraften oder Verformungen im Stiel, sondern nach den entgegengesetzt gleichen 
KraftgréBen im Riegel. Es, ist dann mit (32b) 

Ne a NG = Qs = Cn Wp — Cnm Wb 
(0,— 0.) = Ns = cg us (32c) 
M. — M, = M, = CMN Wp — Cm Yo - 


Nach dem Maxwellschen Vertauschungssatz ist: cyy = ¢nm- 

Zur Bestimmung der Federzahlen denken wir uns den Stiel jeweils so belastet, daB eine der 
FormanderungsgréBen am Stielende b (w,, u, oder y,) gerade = 1 wird, wahrend die anderen = 0 
sind. Die hierzu erforderlichen, am Stielende angreifenden Krafte sind dann die zugeordneten 
Federzahlen. So erhalten wir z. B. fiir den unten gelenkig gelagerten homogenen Stiel (Abb. 13) 
die Federzahlen cy und cywn aus 


vis 
WwW, = 1= ST und CMn = len . 


Bei den zwei wesentlichsten Lagerarten des FuSpunktes, Einspannung und Gelenk, ergeben 
sich so fiir den homogenen Stiel die folgenden Federzahlen: fiir die Einspannung am Stiitzenfub- 
punkt 


EF 12 0B 4 EI 6 EI ; 
CO Fie ° Ni 3 ° Cm >= ie ) CNM = ae ° (33) 


fiir die gelenkige Lagerung am Stiitzenfubpunkt 


EF Bie OW A 3 El SL 33 
C0 ae AN = EB? Cm = Te CNM = Sm ae ( a) 


Die Beziehungen (32a) und (32c) kénnen wir zu einer Matrizengleichung zusammenfassen, die 
zwischen den Vektoren ), an der Stelle c und jp, an der Stelle a vermittelt. Wir setzen noch 


is Z, is iS is = B 
OQ = UNS 9 Cyn = ¢n Fr ’ Cu = “MF ’ CnmM = lars (34) 
und erhalten 
Pig teat tcl 0 0 0 0; | u 
N mnie. 1 0 eet, eed) o| | N 
oe te ah). 0 0 1 0 0 Peels tae (35) 
7, 0 0 0 1 0 0; | py 
) |M | —tnm 0 0 Cm 1 0); | M | 
) ko. | pe Vic, 0 0 Tate On dee 
oder 
) p.] = [Sy] > [Pa] - 


aceon wir den Anteil von u und N weg, ‘streichen also von der Matrix ©y die zwei ersten Zeilen 
‘und Spalten, so fiihrt der Rest zur Gleichung (I, 47), natiirlich ohne m und @, die wir hier un- 
beriicksichtigt gelassen haben. 

b) Biege- und Torsionsschwingungen. Der Abb. 14 kénnen wir wieder sechs geome- 
trische Beziehungen und drei Gleichgewichtsaussagen entnehmen: 


| W,. = Wa = Wh » Q. = Q, ae Q, ? 
Wo = Wa = Xd 9 (36a) M. — M, ae Mr ? (36b) 
Mr. = Mra + M,. | 


Xe = Xa a Wo ? 
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Entsprechend (32c) schreiben wir mit (36b) die Elastizitatsaussagen in der Form 
.—2, = = Co Wp — COT Yb » | 
(M, M.,) = Mr = CM Xb > 


(36c ) 
Mr. Mr. = M; = CTQ Wp — CT Wb » | 


mit 

cra = Cor. 
Ist das AnschluBglied ein homogener Balken, so werden die Federzahlen, wenn der FuBpunkt 
cingespannt ist, 


GI, 12 EL 4 EI 6 EI | 
nye == ae ? co = ies y) Cr ih y) COT sat 2 “5 (37) 
und wenn der FuBpunkt gegen Biegung gelenkig gelagert ist, 
Gly 3 EI 3 EI 3 EI 
ti =- 7 ? CO Ss a? OTe Pe COs pee (37a) 


Aus den Beziehungen (36a) und (36c) kénnen wir wieder eine Matrizengleichung herleiten zwischen 
den Vektoren q, an der Stelle ¢ und gq, an der Stelle a. Wir setzen: 


Pyke to Spe one ee (38) 
Q Q cai” ie 1 EL’ M M EI’ TQ TQ ET 
und erhalten 

yas a te 0 0 0 0 07 ry 1 

| M cr J ero 0 0 0. Mr 

| w Es 0 1 0 0 0 | Ww (39) 

v 0 0 0 1 0 0 yp 

M 0 0 0 ¢m 1 0 M 
10 |, | —ere 0 20) 0 0 1} fO | 


oder 


Abb. 14. Ubungsbedingungen an einer Stiitze. Abb. 15. 


Die Matrizen Gy und Gr kénnen wieder in eine beliebige Feldkette nach Art der Gleichung (18) 


eingereiht werden. Alle Feldmatrizen eines Systems miissen dabei natiirlich auf den gleichen | 


Vektor, im Falle a) auf ), im Falle b) auf q bezogen sein. 


9. Geschlossene Feldkette. Mit den Matrizen § (16) und & (27) fiir den Eckibergang, ‘8 (17) 
fiir Biege- und Langsschwingungen und Gy (35) fiir das tragheitsfreie AnschluBglied haben wir 
die Méglichkeit, ein offenes ebenes Stabwerk zu erfassen, das aus beliebig vielen Feldern (homo- 
genen Balkenstiicken oder Einzelmassen) besteht, in seiner Ebene beliebig oft rechtwinkelig oder 
schiefwinkelig abgebogen ist und beliebig viele elastische, masselose Anschliisse hat. 

Kin Stabwerk kann aber auch so aussehen (Abb. 15), daB die Stelle i mit der Stelle 1 zusammen- 


trifft, d. h. es kann geschlossen sein. Zwischen dem Vektor p, und p; gilt dann die Matrizen- 
gleichung 


[pi] = 20) - [ps], (40) 
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und es miissen, wenn die beiden Felder in derselben Richtung liegen, alle Komponenten der Vek- 
toren ); und ), einander gleich sein: 


= 0s 4] 
Aus Gleichung (40) wird daher i (41) 
[i] = (20) - (pil 


(0] = (29 — (6) - tpl. (12) 
wobei [0] die Nullmatrix und [(] die Einheitsmatrix ist !. 
Soll der Vektor ), von Null verschieden sein, so mu die Determinante von ({QC] — [C]) ver- 
schwinden. Die Eigenwertgleichung heift damit 
4 =| [1] —[6]|=0. (43) 


a ee Stellen i und 1 nicht starr, sondern gelenkig verbunden, so miissen wir statt (40), (41) 
schreiben 


oder 


et eet 

N, | | | IN, 

w | 4 

SESE inert (44) 
ae YI 

On| | | 0) 

Boe eae ers 


Da wir auf der rechten Seite eine Unbekannte weniger haben, kénnen wir die Gleichung mit 7p; 
unterdriicken und erhalten als Eigenwertgleichung 


r 
et) Ay Aj; Ay Ai, | 
Ay, (Ay. — 1) Ap: Ay, Ay, 
A= As, Asp (A33 — 1) As, Ase = (45) 
As, A;. A:; As, As, 
Ay Ago Ags Aga (Ag6 — 1) 


Sind die Stellen i und | durch ein gelenkiges Lager verbunden, so ist die Verbindung nicht mehr 
echt. Die beiden Enden kénnten auch getrennte Lager haben. Zwischen den von Null ver- 
schiedenen Komponenten von )), und ); besteht keine Beziehung mehr, und wir erhalten die Lé- 
sung wie im Fall der Gleichungen (19) und (20). 

Wir sind nunmehr in der Lage, die Eigenfrequenzen eines 
Rahmens zu ermitteln, wenn er gekoppelte Biege- und Langs- 


Abb. 16. In der Ebene schwingendes Rahmenwerk. Abb. 17. Aus der Ebene herausschwingendes Rahmenwerk. 


schwingungen in seiner Ebene ausfiihrt. Abb. 16 zeigt einen solchen Rahmen. Die AnschluB- 
glieder sind als masselos vorausgesetzt. . 

Was wir hier fiir gekoppelte Biege-Langsschwingungen [mit p (17) als Vektor der mechanischen 
GréBen] gefunden haben, gilt in gleicher Weise fiir Biege-Torsionsschwingungen, wenn wir statt 
den Vektor g (30) einsetzen. Die Biegeschwingungen erfolgen dabei senkrecht zur Rahmen- 
ebene, wahrend die einzelnen Balken des Rahmens gleichzeitig Torsionsschwingungen um ihre 
Achse ausfiihren. Abb. 17 zeigt einen solchen Rahmen. Die Stiele seien auch hier masselos. 

Der Rahmen nach Abb. 17 kann aber auch in seiner Ebene schwingen. Dann gilt wieder der 
Vektor ), der sich in diesem Fall von q nicht nur dadurch unterscheidet, daB®B er statt der Tor- 


sions-Komponenten 7 und M, die Langsschwingungskomponenten u und N enthalt, sondern 
Wi 


1 R. Zurmiihl, Matrizen, S.6 und 11, Berlin, Gottingen, Heidelberg 1950. 
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auch dadurch, daB die Biegeanteile in beiden Vektoren ganz verschiedene Bedeutung haben: die 


Biegekomponenten des Vektors ) beschreiben die Biegung in der Rahmenebene, die von q be- — 


ziehen sich auf die Biegung senkrecht zur Rahmenebene. 
Im allgemeinsten Fall kénnen die Schwingungen in der Rahmenebene und senkrecht dazu 
auch gekoppelt sein, und zwar immer dann, wenn das Stabwerk (in seinen Massen oder Stielen 


z.B. Abb.17) nicht ebenen-symmetrisch ist. Um ein solches Schwingungssystem streng zu erfassen, — 
miissen wir die 2-6 Komponenten der Vektoren ) und q zu einem gemeinsamen Vektor und | 


die Matrizen § und 9, zu einer gemeinsamen Matrix zusammenstellen, wie wir dies bei der Koppe- 
lung von Biege- und Langs- oder Torsionsschwingungen gemacht hatten, (17) und (31). 


Nach Méglichkeit wird man durch vorherige Abschatzungen, Symmetriebetrachtungen, und | 


dergleichen versuchen, ein solches System so weit wie méglich zu vereinfachen. 


10. Dreifelder-Anschlu8. Wenn wir die Eigenwertgleichung fiir ein Stabwerk aufstellen wollen, 
das nicht nur aus einem durchlaufenden Balken besteht, sondern auch noch AnschluBglieder auf- 
weist, deren Achsen auf der Balkenachse senkrecht stehen, so ist das mit den bisher behandelten 
Mitteln nur méglich, wenn ein solches Glied an ein Ende des Balkens 
anschlieSt, so daB man etwa mit Hilfe der Eckmatrix 7% (16) von der 


dargestellten Rahmen miifbte also die Masse des mittleren Stieles ver- 


suchen, in die Schwingungsgleichung (Matrizengleichung) des vom 
he eee A es at Lager a zum Lager c durchlaufenden Grundsystems die fiir sich 
a Sticker. SChwingungsfahige Mittelstiitze aufzunehmen. 

a) Wir betrachten zunachst Schwingungen in der Ebene des Stab- 
werks, Als die Unbekannten fihren wir die sechs Vektoren (mit je sechs Komponenten) ),;, );, x 
und ),, )», p- ein. Zwischen ihnen bestehen 36 Gleichungen. Zunachst die 3-3 Randbedingungen, 
die je drei Komponenten der Vektoren ),, ),, ), festlegen, dann die drei Matrizengleichungen 
(= 18 Komponenten-Gleichungen) 


Dp; => pa pj = B° po, pe =C- be (46a) 
und schlieBlich die neun Ubergangsbedingungen am Knotenpunkt (Abb. 18) 


YILLLELL S42 


at 

N= —0,, 

Wy = Wj; > 

i ees (46b, c) 
Wh = YWi> Uy = Wj > 

M, = M; + M,, Wit US 

0, =Q; + N,, Ca Yj» 


die wir ebenfalls in Matrizenschreibweise zusammenfassen miissen. Auf der linken Seite der Glei- 
chungen (46b) steht der vollstandige Vektor p,, auf der rechten Seite der Vektor );, sowie die 


drei Kraftkomponenten des Vektors p;. Diese Beziehungen kénnen wir durch die Matrizen- 
gleichung 


a oe | ular, (47a) 


S89 So © 
[SS SS) SS) 
SQ 


SSS aS S|) 
SSS Ovo S 
i=) 


ms) 


0 45 


darstellen. Die Gleichungen (46c) enthalten die Verschiebungskomponenten der Vektoren p; und 
p;. Wir kénnen dafiir in Matrizenform schreiben 


lO 0 eee i O.- 0) hit Stpue ae me 
0 0 “1 00. top e125 eg 0 eee cman ee ATI 
(ate ee nos eure ire ur ¢ ene 


| 
| 


if 


einen Achse zur anderen ,,umsteigen“ kann, oder wenn die AnschluB- | 
glieder selbst masselos sind und als elastische Stiitzen mit den Ma- | 


trizen ©y (35) bzw. Gr (39) erfaBt werden kénnen. Bei dem in Abb. 18 


nachlassigbar klein sein. Ist dies nicht der Fall, so miissen wir ver- | 
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Wenn wir von den Stellen i, j,k zu den Stellen a, b, c ubergehen, indem wir die Gleichungen (46a) 
verwenden, erhalten wir aus (47a) und (47b) die Matrizengleichungen 


mt) 0 0 0 0 0 | 
| 0 0 0 0 0 —l 
ea. (ache oor 0G 
10), v (2¢ Da 45 0 0 0 0 0 0 * % : )s) o) (48) 
0 0 0 0 i! 0 
a) 1 0 0) 0 04 


1 0 0 0 0 0) 
Ole= 1.0) 0 il 0 0 0 
0 0 0 1 0 0 


0 0 —1 0 0 0 
“W-p,+)1 0 0 0 0 | “B-py. (49) 
0 0 0 —l 0 0 


Wir multiplizieren im Inneren der Gleichungen (48) und (49) die Matrizen [ ]-% und [ ]- 


aus, fassen die beiden Gleichungen zusammen und setzen den gemeinsamen Vektor i heraus. 
Jt, 


Dann wird aus (48) und (49) 


i 1 r 0 0 0 0 0 0 iT ae, 4 

| — 6 — By Bes Bos By; Bee Nee 

Ic | , vt 0 0 0 0 0 We 
ee 0 0 0 0 0 0 Wa 

Uae | Bs, Bs. Bsg Bsy Bss Brg Mz | 
Sek | at Be Bas Bag Bas, Bag 1) (| Oe 
Ay, Aj, Aygy Ayy Ajs Aig — Day = Bs, Be, — Baa —B,; —B3, \ eee 

Az, Asp Aggy Az4 Ag; Az, 11 By By By B,; By, No | 

J i Ay, Agy Aggy Agy Aas Ag 41 42 By By, By; By, see oan 
| Yo 
| M, 
Qs 


Formulieren wir nun die Randbedingungen, so werden rechts je drei Komponenten der Vek- 
toren ), und ), Null (wenn wir von elastischen Lagern absehen, deren Matrizen aber vor die Rand- 
vektoren gesetzt werden kénnen) und zusammen bleiben sechs unbekannte Komponenten wbrig. 
Links werden drei Komponenten des Vektors }), Null, so daB die linken Seiten von sechs Glei- 
chungen des Systems (50) verschwinden. Diese sechs Gleichungen fiir die sechs iibrig gebliebenen 
Komponenten von ), und }, stellen ein homogenes Gleichungssystem dar, dessen Determinante 
die Eigenwert-Determinante des Schwingungssystems ist. 

Sind in einem speziellen Fall die Randbedingungen gegeben, so kénnen wir in der Gleichung (50) 
alle die Zeilen bzw. Spalten der Matrizen Y, B und © streichen, die zur Aufstellung der Kigen- 
wert-Determinanten nicht gebraucht werden. 


Fur den in Abb. 18 dargestellten Rahmen erhalten wir so 


[ f Aye Ai; Aj 0 0 0 4 ] F : 
Cy Cy Cy Cy Cys Cig Ay, Ao; Age Bez —Be, —Beg | | | Na 
Cz, Cy Cup Cg Cs Cy6 As, As, Age | 0 0 0 | | M, 
4 a2 yg Cay Cys 46 Ag Ay, Ai, 0 0 04 | Qu 

| As, Ass Ase  Bsp Bog Bog | |-| Ne |- (51) 

k Ags Ags ; Ag Boo By, . Bog = | Yo 
Ay Ay; Aye Bs. —Bz_ —Ba¢ Qr 

Az Ags Ase By By 1d 

Ay, Ags Agg —Biyp—By —Bag 4 


Wir kénnen uns demnach schon bei der Aufstellung der Matrizen 2{, 8 und © darauf beschran- 
ken, nur die hier bendtigten Elemente zu bestimmen, also von jeder Matrix nur drei Spalten bzw. 
Zeilen. Multiplizieren wir nun noch die inneren zwei Matrizen von (51) aus, so erhalten wir drei 
Zeilen von je sechs Elementen, die zusammen mit den darunterstehenden drei Zeilen die Eigen- 
wert-Determinante liefern. 


. (50) 
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f 
b) Wir wollen jetzt einen Dreifelder-Rahmen betrachten, der senkrecht zu seiner Ebene | 
schwingt (Abb. 19). | 
Zwischen den drei Randern a, 6, c und den drei Knotenstellen 1, j, k bestehen jetzt die drei | 
Gleichungen 


Gp =) Gas q; =U- a, qe = - Gr. (52a) | 


wobei die Matrizen 2, 8, € und die Vektoren q die Biege- und Torsionsanteile nach (31) enthalten. | 
Als Ubergangsbedingungen am Knoten erhalten wir hier die neun Beziehungen 


Uk aa Ws ? 
Mr, = Mr; — M,, 1; =e 
a? “i die (52b, c) 
Pk = Yi» Ve eis 
M, = M; + Mz;; 
Abb. 19. Rahmen mit massebehafteten Stielen. 0, Pa 0; -+- Q; 


zwischen den Stellen i, j, k [vgl. (36a) und (36b)]. 


Wir gehen genau so vor wie unter a) und erhalten die allgemeine Matrizengleichung fiir Biege-_| 
und Torsionsschwingungen | 
ae | 
| Mra 
fe ven aie: 0 0 0 0 0 O° i aves 
Mr, | = By Be Big ee oe | Pa 
We 0 0 0 0 0 0 M, 
= Ge a as 
We 0 0 0 0 0 0 Qu 
Vs By By By By By, Bag | Xo 
Q. Be Bep Ben Bes Bes Beg |) | Mrv | 
we See Ay, Ay Ay Ay Ayy Aye: By Ba Be Bu Be Ba Wy | 
0 As; Ag, Agg Aga Ass Ase —Bz, —B3, —B33 —B3, —B;; —Bgg | Bb 
| Ag, Asg Aggy Ag Ass Age : —By —By —Big —By —By, —Bye | | M% 
LQ | 
Diese Gleichung vereinfacht sich wieder durch Einsetzen der Randbedingungen auf die Halfte 
des Umfangs, wenn wir dabei nur die Komponenten-Gleichungen mitnehmen, deren linke Seiten 


verschwinden. Die Determinante der ibrigbleibenden Matrix ist dann wieder die Eigenwert- 
Determinante. 

Wir wollen als Beispiel einen Rahmen (Abb. 20) betrachten, 
der zu einer Vertikalebene durch a, b,c symmetrisch ist. Die 
Langs- und Quertrager kénnen im iibrigen ganz beliebige Massen- 
und Steifigkeitsverteilung haben. Die Stiitzen seien masselos. 
Thre Langs- und Biegesteifigkeiten kénnen mit Hilfe der Matrix 
Gn (35) beriicksichtigt werden. Wir kénnen (beziiglich der Sym- 
metrieebene) symmetrische und antimetrische Schwingungen 
getrennt behandeln und den Ansatz jeweils auf eine Hialfte des 
Rahmens beschranken. Fiir eine solche Halfte sind die Stellen 
a, b, c Randstellen, und es gelten die folgenden Randbedingungen fiir jede der drei Stellen a, b, c: 


Abb, 20. Raumlicher Rahmen 
(Stiele masselos). 


1) Symmetrische 2) Antimetrische 
oe . . . . 
Eigenschwingung: Eigenschwingung : 
Mr Ea) ‘| ie NG | 
y= 8 (54b) 0 (54a) 
C= 077 Whales | 


Damit erhalten wir aus (53) die beiden Eigenwertgleichungen 
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1) fiir symmetrische Eigenschwingungen: 


PA ate 4c 0 0) 0) 
Cy Cy Cy Cy Cys Cog | Any A,, A,; —Bs, Bs; —Bss 
a Cy = Cyg Cy 55 Cag|| 4g, gg Ags 0 0 0 
61 62 Cos Cor Cos 66!| Ay, Agg Aas 0 0 0 
i (As, Ass Ass By, Bas Bo, || =0, (55) 
Ag, Ags Ags Bey Bes Bes -1| 
Ate Ara Big Bas By | 
As, Ass Ag5 Bs, Bss Bs; 
Ay, Ay; Aj; By, By B,; 
2) fiir antimetrische Eigenschwingungen: 
a | Ayp Ain AG 0 0) Qo 
Cy Cp Cis Cis Ci, Ci Ay, Any Aog Bs, Bs, Be 
31 Cp C53 34 C35 cn As, As, Age 0 0 0 
51 Oss 54 55 56 ao Aggy Ag 0 0 0 
7 52 As. As, By By Bag | =0. (56) 
Ao Aga Ave Bo Bes Beg J 
Pe AeeA eee. a. bee Dee 
13, Ag, Age By, Bsa Bygg 
Ay, Ay Ase: —Biy By Bye 


Wollen wir beide Schwingungsarten erfassen, so brauchen wir die vollbesetzten Matrizen YI, %, . 


ll. Zweifach geschlossene Feldkette. a) Biege- und Lingsschwingungen. Wenn wir die 
beiden Seitenarme des in Abb. 18 dargestellten Rahmens mit den Enden a und c bis an das Ende b 
des Mittelarms heranfiihren und die drei Arme a, b, c verbinden, so entsteht eine zweifach geschlos- 
sene Kette. Von irgendwelchen Stiitzen kénnen wir zunichst absehen, da fiir jede Stiitze ins- 
besondere jede elastische Stiitze eine entsprechende Matrix in die betreffende Abschnittsmatrix 
aufgenommen werden kann. Auch fiir diesen Rahmen gilt die G1. (50), und zwar mit den folgenden 
Ubergangshedingungen: 


U, = Ug 5 

Nee One, Uig= Oi 

a6 = Wa p me aa ie ae z (57a, b) 
We = Wa > Wo = Wa > 


M,=M,+ M,, 


0. = Q.— Ms, 


Setzen wir (57a, b) in (50) ein, so kénnen wir schreiben 
7 
| 
| 


Ug ] ee De Lae Daw | Ua 
IN = Q Dz, Na 
Wa ° Wa 
Pee ple | | Yo 
M,+M, ; Dp car MM, 
Q.—N _ . j-| @ |, (58) 
0 , | Wa | 
0 Nz 
L 0 a IL Dy, a —U, 
tn 
M, | 
i QO, 3 


also, in Matrizenform, neun Gleichungen fiir neun Unbekannte. Um die Eigenwertdeterminante 
zu erhalten, miissen wir die linksstehenden Komponenten von }, und ))p auf die rechte Seite bringen 
und alle Elemente zusammenfassen, die zu gleichen Komponenten der Vektoren }, und }), gehéren. 
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Soweit dabei in der Matrix D, die Elemente der Matrizen 9( und § erhalten bleiben (in den un-- 
tersten drei Zeilen), wollen wir sie auch mit A,,, bzw. B,,, bezeichnen, wenn wir D, ausschreiben: 


IP (D,, = D,,—1) D,, (Dis t D,,) (Dis 1 Dy 10) Ds D,, Dy Di Dy1> 71 Ir va 1 
| (Oye Dy) (Ds, — 1) (Dy; = D,,) (D,, <t Dy 49) Dy; D,, Dy, Ds43 (D5,12—1) Na 
(Day aa Dayo) Day (D3 site D,; —1) (Des =F Ds,30) D,. D,, D3 Ds,15 Ds, 19 wa 
| (Oe D;.) D,, (Py; at Dy) (Dys + Dy. —) D,; D,, D,s Dyas Diyo32 Ya | 
0= | (Di — Dyp) Ds» (Dss a9 D;,) (Ds. AP D510) (D5 xy 1) D;, Ds» (D3,11 an 1) Ds.» 2 Ma 0 (59) } 
(Dor — Doo) Ds, (Dos se D,;) (Dos i D, 10) D, (Doo aa 1) (Dos f 1) D, 11 D. 12 Qa 
(Ay, sts By;) Aj, (4,3 — B,,) (Aq, a B;,) A,; Ax, —B;, —B,,; = By, Np 
| (Ay, — By;) Ag, (Az3 a B,,) (Ag, Se B,,) A 55 As, By, B, Bi, Mp 
L (Ay ate B,;) Ay, (Ay, — B,) (A\,— B,,) A; Ay, —B,, ' &B,; —B,, = - Qp 4 


Mit Hilfe der Gl. (59) kénnen die Eigenfrequenzen von Rahmen etwa von der Art der Abb. 21a bis) 
21c fiir Schwingungen in der Rahmenebene bestimmt werden. Die stark gezeichneten Teile kénnen 
dabei ganz beliebige Massen- und Steifigkeitsverteilung haben, solange sie nur symmetrisch zur_ 
Rahmenebhene bleiben. Die einfach gezeichneten Stiitzen sind masselose AnschluBglieder. Bei dem 
Stabwerk der Abb. 21¢ mu vorher untersucht werden, ob die Koppelung mit Schwingungen senk- 
recht zur Rahmenebene und mit Torsionsschwingungen vernachlassigbar ist. 


SS 
LLLLLLLLLLLILLLILLLLLILILLL 


Nu 
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Abb. 21. 


b) Biege- und Torsionsschwingungen. In gleicher Weise wie wir fiir den zweifach ge- 
schlossenen Rahmen aus (50) die Gl. (59) fiir Schwingungen in der Rahmenebene gefunden haben, 
kénnen wir auch fiir Schwingungen senkrecht zur Rahmenebene die entsprechende Gleichung 
aufstellen. Wir gehen aus von der Gl. (53), denken uns die beiden freien Rander a und c (Abb. 19) 
mit dem Ende b des Mittelarms vereinigt und formulieren hierfiir die Ubergangsbedingungen, die 
wie unter a) in die Grundgleichung (53) eingesetzt und geordnet zu der Gleichung 


ff (Dy + D,,19—1) D,» (Dis + Dy») (D,,— D,;) D,; D,, D,, Diy D,,, aq T Xa 
(D,, aig Dy,10) (D..—1) (Dy5 ala Dy») (Do, — D,;) Dy; Dy, Days (Dp,11 —h) Dz, Mya 
(Ds, qe Ds 10) D3, (Dy; tS D;;—1) (Ds, ee D;;) Ds; D,, D3" Di aa D319 Wa 
(Dy ar D,,10) Dy (Dy; ap Dy») (Ds; — Dy—1) Dy; D,, D,s Dy Dy.15 Va 

0= (Di at Ds,10) D;» (Ds lz D5») (Ds = 1) (Ds; — 1) Ds, (Ds 1) Ds Ds. | Ma (60) 

(De: 39 Dg,10) Dy» (Des iP Dy») (Dos =—4 D;;) D, (Dgs a 1) Des Deyn (De,12—1) Qu 
(4, + By) Ay, (413 + Bys) (4.— Bu) 1, Ay, B,, B,; Bye Mr 
(43, — B,,) A,, (A,3— B;;) (Ag, ap B,,) A 3 A 6 ==! ee B;; —B;, My 

L (Ag — By) Ay, (Ays; — By) (Aj, + By) Ai; 14, tb a= lB, —B,, a L Qp el) 


fihren. Diese Gleichung gilt z. B. fiir einen Rahmen nach Abb. 20 ohne die dort vorgeschriebene 
Einschrankung der Symmetric, allerdings wieder nur dann, wenn keine nennenswerte Koppelung | 
mit Schwingungen in der Rahmenebene eintritt (siehe Ziff.9, Schluf). | 


12. Schiakwort. Mit Hilfe des handlichen Matrizenkalkiils lassen sich Balken von beliebiger — 
Felderzahl mit starren oder elastischen Stiitzen sowie offene oder geschlossene Rahmen mit Koppe- 
lung verschiedener Schwingungsarten exakt erfassen, ohne daB dabei der Uberblick verloren geht. 
Die Grenzen, die dem praktischen Rechner heute noch gesetzt sind, werden durch die standig fort- 
schreitende Technik des Maschinenrechnens immer weiter hinausgeriickt. Es ist zu erwarten, da} 
damit auch unser Kalkiil mit der immer wiederkehrenden einfachen Operation der Matrizenmulti- 
plikation immer miiheloser angewendet werden kann und daB die hier gezeigten Ansatze dann noch 
weiter ausgebaut werden kénnen. 


(Eingegangen am 7. Juni 1955.) 
Anschrift des Verfassers: Dr.-Ing. Hartmut Fuhrke (14a) Fellbach (Wiirttemberg), Heinrich-Heine-Str. 8. 
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Anwendung von komplexen Verinderlichen zur Synthese einer 
Kurbelschwinge mit vorgeschriebenen Grenzen 


der Abtriebs- Winkel-Geschwindigkeit 


Von N. Rosenauer 


1. Die Lésung. Das Erscheinen der interessanten Arbeit von K. H. Sieker: », Winkel-Langen- 
Zuordnung an der Schubkurbel, Darstellung mittels komplexer Zahlen“* im September-Heft 1954 
Jer ,, Konstruktion“: veranlaBte den Verfasser die Anwendung der komplexen Veranderlichen zur 
Lésung eines anderen Problems zu zeigen. 

Dieses Problem wurde erstmalig von 
S. Bloch! behandelt. Eine Lésung wurde 
auch vom Verfasser 1944 veréffentlicht 2. 
Da aber voraussichtlich kaum etwas von 
dieser Auflage erhalten geblieben ist, so 
ware es zweckmaBig, eine Lésung des Pro- 
blems hier zu zeigen, in der vom Ver- 
fasser wesentliche Vereinfachungen ein- 
vefiihrt sind. Lo 

Die komplexen Bezeichnungen der Glie- ~ Lee 7 
Re Celenkvierecks mach Abh ol ind Abb. 1. Gelenkviereck mit eingetragenen Bezeichnungen. 
OA=G=aei%, AB=b=bei®, DB=¢ =ce'®, DO =d=de'*, wo O,, O,, Og ver- 
anderliche Winkel darstellen, aber 6 ein konstanter Winkel ist, der spiter bestimmt werden soll. 

Das geschlossene Polygon des Gelenkvierecks ist durch die folgende Vektor-Gleichung aus- 
gedriickt: 


oo bee ld = 0. (1) 

Die erste und zweite Ableitung dieser Gleichung nach der Zeit liefert 
O,4+0,b6—w,e+0-d=0, (2) 
(@ 6,—w2) 4 + (t2,—@?) b— (62, wi) 6 +0-d=0, (3) 


wo W, = 0, » We = 0, , 3 = O, die Winkelgeschwindigkeiten der Glieder und e, = @, = 6,, 
a CO), = 0, , €3 = W; = O, die Winkelbeschleunigungen der Glieder bedeuten. 
Die Auflésung des Systems (1), (2), (3) ist aus der Matrix abgeleitet 


1 1 — 1 il | 
M =| Wy Dy — Ws 0 | , 
© . 6 . Dy 
Pie. i &— W3 —t &, + W3 0 


wobei der Proportionalitatsfaktor gleich 1 gesetzt ist: 


al 


= Wg Wy (3 — We) + (Wg & — My &) T, 

= () Wz (Wy — Ws) + (@ 3 — Wg €y) F (4) 
= Wy Wy (dy — Wy) + (Wy &y — Mg &1) t- | 

Bei gleichférmigem Antrieb des Gliedes OA ist ¢, = 0, und fiir die GréBt- und Kleinstwerte 


der Abtriebs-Winkel-Geschwindigkeit ist auch e, = 0. 
Dies vereinfacht die Ausdriicke (4) in die folgende Form: 


oS] 


= Wy Wy (W3 — We) + Ws Et » | 


= (, Ws (WH, — Ws) » (5) 
= Wy Wy (Wy — Wy) + 4 Egt- 


mS. Bloch, Zur Synthese von viergliedrigen Mechanismen. Nachrichten der Akad. der Wiss. U.R.S.S. 


i. N. sieshidaer Gelenkvierecke mit vorgeschriebenen GréSt- und Kleinstwerten der Abtriebs-Winkel- 


eschwindigkeiten, Mabau.-Getriebetechnik 12 (1944) S, 25—27. 


ol ol el 
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Die Ausdriicke (5) stellen die Glieder zweier verschiedener Gelenkvierecke dar; das eine ent+ 
spricht max j, und das andere min w,, die in verschiedenen Mafistaben erscheinen. | 

Bezeichnet man die Glieder des ersten Gelenkvierecks mit a,, b,, c;, d, und die Glieder des) 
zweiten mit dy, b,, ¢, dy, so hat man zwei Gleichungssysteme 


Gy = Wg1 Way (Wg; — Wy) + Ws Eo1 T | Gy = Wg Wop (C32 — Wop) + Wg E221 » 
by = Wy Wg; (1; — 1) » | by = @y Wg (1, — gp) » (6)| 


Cy == Wy Woy (Wy — Mo1) + y €q, T; Co, = Wy Wop (Wy — Wop) + Oy Eog ¥ 


Falls die Gleichungen (6) dasselbe Gelenkviereck in zwei verschiedenen Lagen darstellen sollen, 
miissen die MaBstabe gleichgesetzt werden. Dies wird erreicht, indem man b, = k b, setzt. Hier- 
aus folgt das Verhaltnis der MaBstabe k = b,/b, . 

Die Gleichsetzung der Quadrate der Glieder a und c unter Beriicksichtigung des Koeffizienten ke 
ergibt die folgenden Gleichungen: 

Ay ee Shee By aes 
CD51 C03, (Dg, — Wy)” + C3, E31 = Ke? C3» C32 (Wag — Map)” + hk? wp E50 » (7) 
CD31 (Cy — Wy) + £3, = KP Why (Oy — Wop)? + he? €> « (8) 
Die dritte Gleichung folgt aus (1) d = ¢ —a—b mit d,=kd, zu 
2 ye 
[204 (a4 (1 — 1) — (31 Woy (031 — Hg1) — 1 M1 (21 — gy) ]? + ( — 1)” €24 
Pp 

= K? [cy yp (01 — Wy) — Wgy (yp (C32 — Wop) — Cy Wap (021 — Wp) |? + Ke? (2, — Wp)? €2 - (9) 
Betrachtet man ¢,, und & als Unbekannte, so fiihrt das zum folgenden System von drei Glei-; 
chungen: | 

CD31 E51 — KP We E32 = kh? Why Wp (Wg2 — Hyp)” — 093, ", (31 — My)? 5 
£3, — h? &3, = FP oy (©, — Cdy9)? — 03, (Oy — Wei)? 5 
(101 — gy)? €3, — kh? (coy — Wyo)? 5 (10) 
= 2 
= FP [004 Wyo (Wy — Woy) — gp Woe (W232 — x2) — 4 Wap (1 — Ms2)] 
2 
— [@y C1 (2, — Gp) — 31 Woy (031 — We) — 4 Hg (M1 — M31) }? «J 


) 
| 
Die Elimination von ¢,, und ¢,, aus diesem System fiihrt zu einer kubischen Gleichung, die wy, | 
und W.. enthalt. 

Da man beweisen kann, daB 


“ = £3 = Ws (Ws = () (, +- Ws; — 3 (Wy) 


ist, so geniigt es, einen Wert fiir w,, so zu wahlen, daB die Ungleichung 
(gy (Wg, — 4) ( + Wg; — 3 Woy) <0 
fiir ein max «, befriedigt ist. Setzt man diesen Wert in die kubische Gleichung, so kann man all- 


gemein drei Wurzeln fiir w,, berechnen. 
Um die Giiltigkeit dieser Wurzeln festzustellen muB die Wurzel erstens die Ungleichung 
| 


(O32 (Wg — Gy) (Wy + W32 — 3 Wop) > 0 


fiir ein min ~, befriedigen, und zweitens miissen die Gleichungen (10) fiir ¢5, und ¢3, positive GréBen | 
ergeben. 

Falls diese Forderungen erfiillt sind, kann das gesuchte Gelenkviereck in beiden Lagen unter | 
Beriicksichtigung des Koeffizienten k gezeichnet werden. | 


Der soeben beschriebene Gedankengang der Lésung wird an dem folgenden Zahlenbeispiel | 
erlautert. 


2. Ein Beispiel. Gegeben sei 


Op =) 10 seen eee On 
One 6 sec , enue On 

pate Biscort =p 

Wa Seca £55 = 0k 


Das entsprechende Gelenkviereck soll in den zwei entsprechenden Lagen fiir die vorgeschriebenen 
Grenzen der Abtriebs-Winkel-Geschwindigkeit bestimmt werden. 
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Erstens berechnet man das Verhaltnis der Ma®stibe k aus den gegebenen GriGen: 


y hye by @y Mg, (@; — gy) _ | 6(10—6) 1 


CP di. (Oi — 4) 8 0 8) 


Setzt man den Wert von k und die gegebenen Werte von w,, und (go in die Gleichungen (10), 
so erhalt man 
81 €3, — 4 €5, = 4 wd, (8 + Wye)? — 81 w3, (6 — w,,)? , 
36 3; —€3 = 3 (10 — wy)? — 36 w3, (10 — w,,)?,, (11) 
16 €3, — 9 by = 9 (Wy — 2 Wy, — 80)? — 16 (w3, — 16 w,, + 60)? . 


Die Elimination von €), und ¢,, aus diesen Gleichungen liefert die pests die gleich Null 
gesetzt werden muB: 


81 4 4 Wp (8 + Woy)” ae 81 31 (6 Yay Wo)” 
36 1 (39 (10 — ao)? — 36 w , (10 — wo)? =O 
16 9 9 (W5. — 2 Woy — 80)? — 16 (wi, — 16 w,, + 60)? 


Die Auswertung der Determinante fiihrt nach einigen Vereinfachungen zur Gleichung 


T Wy + 20 w3, — 36 w., — 42 w3, + 372 w3, — 3840,, = 0. (12) 
Der Wert von w,, kann jetzt so gewahlt werden, da® die Ungleichung 
(gy (31 — 4) (My + 31 — 3 Woy) <0 (13) 


befriedigt ist. Diese gibt 

— 24 (16 —3,,) <0. 
Es ist vorteilhaft w,, = 0 zu wahlen, da in diesem Falle die kubische Gleichung (12) zu einer qua- 
dratischen reduziert wird. Eine Wurzel ist dann w,, = 0. und die zwei iibrigen erhalt man aus 


der Gleichung 
T Wo + 20 Wy — 36 = 0. 
Dies ergibt @,, = 1,25 und W.. = — 4,11. 
Die giiltige Wurzel muB erstens der Ungleichung 
| gq (32 — 4) (1 + Wg — 3 Wop) > 0 (14) 


geniigen und zweitens aus den Gleichungen (11) positive Werte fiir ¢3, und ¢3, liefern. Die Unter- 
suchung der Wurzel w,, = 0 zeigt, daB diese keine tauglichen Werte fiir ¢3, und ¢3, gibt. Die Wurzel 


Mg, = 1,25 geniigt der Ungleichung (14) nicht. Die Wurzel w,. = — 4,11 befriedigt die Unglei- 
chung (14) und gibt positive Werte fiir ¢5, und ¢3, aus den Gleichungen (11). Man erhalt dabei 
See 91,200 oo, eee 

eo 10 180,00" == E-61,1 see . 


“Mit Benutzung der oberen Zeichen erhalten die Glieder des gesuchten Gelenkvierecks in der ersten 
Lage aus (6) folgende komplexe Ausdriicke: 

a, — Ms; Yo4 (3; — Woy) + Ws, €5) 1 ae 0 4- 6. 147 = 84. 1 3 

b, = @; Mg; (, — 31) = 10- 6 (10 — 6) = 240, 

Cy = Dy Wo (Wy — Moy) + Hy Gt = 0410-147 = 1407, 

. d, = ¢, — a, — b, = 140i — 84 i — 240 = — 240 + 56i. 

Das entsprechende Gelenkviereck 0,4,B,D, in Abb. 2 zeigt die Lage, in der eine maximale Ab- 
-triebs- Winkel-Geschwindigkeit auftritt. 

Um die zweite Lage des Gelenkvierecks fiir die minimale Abtriebs-Winkel-Geschwindigkeit 


zu bestimmen, muB in die komplexen Ausdriicke der Glieder (6) der Koeffizient k = — 1/6 ecin- 
fe ubrt werden, was die folgenden Resultate ergibt: 


: 
: 
: 


G, = Ke Way Wyo (C3 — Wo9) + Ke Wyp 99 t = 21,3 — 81,51, 

b, = k ay Wag (M1 — Wg) = 240 , 

Cy = k Ody Woy (Wy — Gog) + ke OD Exg 1 = 96,6 + 101,87, 

d, = t, — 4, — b, = 96,6 + 101,8 i — 21,3 + 81,5 i —240 = — 164,7 + 183,37. 
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Das entsprechende Gelenkviereck fiir die minimale Abtriebs-Winkel-Geschwindigkeit 0,4) B’D, 

ist in Abb. 2 gezeigt. Falls dieses zweite Gelenkviereck um O, so gedreht wird, daB die Stege zu- 

sammenfallen, erhilt man dasselbe Gelenkviereck in seiner zweiten Lage O,A,B,D,. Die Kon- 

trolle der erhaltenen Lésung ist in Abb. 3 durchgefiihrt, wo beide Lagen des Gelenkvierecks ge- 
zeigt sind. 


At las 5, ase te In der ersten Lage 0,4,B,D, ist die Ab- 
Va ‘ . e 3 triebs- Winkel-Geschwindigkeit 
. a, 
/ 5 i Op = ee, = 10= Osee=, 
__| reelle Achse ft | ; ee |! 1; | 
-10, 50 Oy R50 | 00 0 0 30 wie vorgeschrieben. Auferdem sind die’ 
SNS j, HE tze. Geschwindigkeiten der Punkte A, und B,, 
KG Aes A DISS eres) einander gleich, also die relative Normal-- 
Ae ae he 8 beschleunigung der beiden Punkte gleich) 
4-700 x yo Null. Dies bedeutet, daB die Beschleuni-- 
Re x SS lke gung des Punktes B, mit dem Gliede B,D, 
ney, ‘SC eo & zusammenfallt und folglich die Lage einem: 
S ee Extremwert der Winkelgeschwindigkeit ent- - 
S Vy spricht. 
me In der zweiten Lage 0,A,B,D, ist die} 
Abb. 2. Zwei Lagen des gesuchten Coreae aero mit (0, — 10 pele Abtriebs- Winkel- Geschwindigkeit 
maxX (3 i eo eee ere ae8 epee Poe i | 
i Ws. = — rai Oy aaa 10 =— 8 sec 4 
aad 


wie vorgeschrieben. 

Die Beschleunigung des Punktes B, ist 
graphisch konstruiert! und ergibt sich als, 
B,N,, wiederum in der Richtung des Glie-| 
des B,D,. Dies bedeutet, da auch in 
dieser Lage ein Extremwert der Abtriebs- 
Winkel-Geschwindigkeit auftritt. 


Da die Wahl des Wertes von ~,, nur 
durch die Ungleichung (13) begrenzt war, 
ist die Anzahl] der Lésungen unendlich grof. 
Dies bedeutet, da auBer den vorgeschrie- 
benen Bedingungen noch andere erfiilli 


Abb. 3. KontrolJe der erhaltenen Lésung mittels graphischer werden konnen. 
Konstruktionen. 


Nimmt man fiir ¢,, und ¢,, die unteren 

Zeichen. so erhalt man ein Spiegelbild desselben Gelenkvierecks, also eigentlich dieselbe Lésung. , 
Einen anderen Beweis dafiir, daB die erhaltenen Getriebelagen den extremen Abtriebs- Winkel- 
Geschwindigkeiten entsprechen, liefert die von F. Freudenstein in seiner Arbeit? festgestellte | 
Higenschaft, daB in den genannten Lagen die Kollineationsachse zur Koppel senkrecht steht. , 
Diese Kontrolle ist in Abb. 3 nicht durchgefiihrt, sondern dem Leser iiberlassen. 


(Kingegangen am 16. Juni 1955.) 


Dr.-Ing. N. Rosenauer, A.M.I.E. (Aust.) N.S.W. University of Technology 
Broadway, Sydney, N.S. W. Australia. 


‘ N. Rosenauer und A. H. Willis, Kinematics of Mechanisms. S. 212—218, Sydney 1953. 
> I. Freudenstein. Trans. A.S.M.E. Paper No. 55—SA—20. 
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Uber die Verwendung von Bipolarkoordinaten zur Lésung 
einiger Probleme der Plattenbiegung 


Von S. Woinowsky-Krieger_ 


1. Einleitung. Die Halbkreisplatte wird gewéhnlich als ein Sonderfall der Kreissektorplatte 
betrachtet und das Problem ihrer Biegung dementsprechend in Polarkoordinaten behandelt. Ist 
nun die Halbkreisplatte an ihrem Rande eingespannt, so mu man sich dabei allerdings mit 
Naherungslésungen begniigen. Fiir die gleichférmig verteilte Belastung wurden solche Lésungen 
von R. Kuno1, von H. R. Hassé? und schlieBlich von H. D. Conway und M. K. Huang? ange- 
geben. Hassé benutzt in seiner Arbeit eine von A. Weinstein herriihrende Methode und beschrankt 
‘sich dabei auf die Berechnung der Durchbiegung der Platte. Conway und Huang gehen dagegen 
von der elastischen Flache einer gleichférmig belasteten eingespannten Kreisplatte aus und bringen 
die Durchbiegungen lings eines Durchmessers dieser Platte niherungsweise zum Verschwinden, 
indem sie daselbst entweder eine kosinusférmige oder cine aus passend gewahlten Einzelkraften 
bestehende reaktive Belastung anbringen. 

Die vorerwahnten Lésungen kénnen auch auf eingespannte Kreissektorplatten, zum mindesten 
auf solche mit bestimmten Zentriwinkeln ausgedehnt werden. Speziell fiir die eingespannte Halb- 
kreisplatte 1aBt sich aber auch eine Lésung in geschlossener Form angeben ®, wenn man an Stelle 
der Polarkoordinaten Bipolarkoordinaten verwendet. Fiihrt man diese letzteren ein, so erscheint 
es wiederum naheliegend, die Problemstellung zu verallgemeinern, und zwar in einer Weise, die 
der Eigenart dieses Koordinatennetzes entspricht. 


2. Die zwischen zwei Kreisbégen eingespannte Platte mit gleichférmiger Belastung. Die in 
Abb. 1 dargestellte Platte sei zugleich auf rechtwinklige Koordinaten x, y und auf die Bipolar- 
‘koordinaten €,7 mit den Polen in P, und P, be- 


J a Sin & ae asin) l 
~ Cof E+ cosy’ Y ~ Gop E+ cos 7 us 
Kin durch die Gleichung 7) = konst. bestimmter 

Kreisbogen P, P, hat, wie man sich leicht tiberzeugt, 

‘einen Halbmesser a/ | sin | und einen Zentriwinkel 

||274|. Die durch die Kreisbégen 7 =a und 7 ==) 

\(B>«) begrenzte diinne elastische Platte von kon- 

stanter Dicke mége langs ihrer Randlinie ein- 

‘gespannt sein und trage eine gleichférmige Belastung von der Intensitat q. Ist ferner D die 

\Biegesteifigkeit der Platte, so hat ihre elastische Flache der Differentialgleichung 


AAw = + (2) 


Abb. 1. Die durch zwei Kreisbégen begrenzte Platte im 
Koordinatensystem. 


“za geniigen. Der Ausdruck fiir die Durchbiegung w der Platte mége von der Form 


W = Wy + Wy (3) 


sein, worin 


q 2 22 
Wy = 6A D (a ae vy?) (4) 


sine partikulare Lésung der Gleichung (2) ist und die elastische Flache einer gleichformig mit 
| belasteten eingespannten Kreisplatte, die die Strecke P,P, zu ihrem Durchmesser hat, dar- 


1 R. Kuno, J. Soc. Mechanical Engineers of Japan 37 (1934), S. 101. 
2 H, R. Hassé, Quarterly J. Mechanics and Applied Mathematics 3 (1950), S. 271. 
3 H. D. Conway und M. K. Huang, J. Applied Mechanics 19 (1952), 5. 5. 

4 4. Weinstein und D. H. Rock, Quarterly Applied Mathematics 2 (1944), S. 262. 
5 §. Woinowsky-Krieger, J. Applied Mechanics 22 (1955), S. 129. 
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stellt. Indem man nun von den Beziehungen (1) Gebrauch macht, kann man die Gleichung (4) 
auch in der Form | 


qa cos? | 
(Co) & + cos 4) wy = =i D Gof é aa 1 0°) 


umschreiben. | 

Die zur Erfillung der auf 7 =a, / vorgeschriebenen Randbedingungen noch erforderliche/ 
Zusatzlésung w, driicken wir gleich in den Variablen £, 7 aus. Die Funktion w, hat dann dex| 
homogenen Differentialgleichung 


(Gof € +0087)" (ger + Gal] t= 9 (6) 


zu genigen. Geht man mit dem Ansatz 
(Cof &+ cos y) wy, = e'4* y(n) , (7}| 


in dem 4 eine beliebige reelle Zahl bezeichnet, in die Gleichung (6) ein, so lassen sich die Varia- 
blen € und 4 in der Tat trennen, und man erhalt fiir die Funktion (7) die Lésung 


pn) = C,Coj An cos 7 + C, Sindy cosy + CyCojAy sinn + C, Sindy sin , (8) 


wobei die Koeffizienten C,, C,, C,, C, zwar nicht von 7, wobl aber von dem Parameter/ abhangen. 
Um schlieBlich die Funktion w, den im unendlichen Bereich — co < &< + 00 vorgeschriebenen, 
Randbedingungen anzupassen, moge sie in der Form 


(Coj € + cos 7) wy, = ZO) cos) — d, (9) 
i | 

gewahlt werden. Die Funktion e'** des Ansatzes (7) wurde hier durch cos & ersetzt und zwa ! 
mit Riicksicht auf die durch die Art der Belastung bedingte Symmetrie der elastischen Flache; 
hinsichtlich der Achse & = 0. 
Entwickelt man nun den rechten Teil der Gleichung (5) als eine Funktion von & ebenfalls in 
ein unendliches Integral, indem man von der Fourierschen Integralformel Gebrauch macht, so 


lautet das Ergebnis 


(Co) & + cos y) wy = q.a* cos” 7 (oe cosdé dA 


8Dsinn | Gindn ~* (10) 
f) 

FaBt man jetzt die beiden Lésungsanteile (9) und (10) zusammen, so erhalt man fiir die Gesamt- 

durchbiegung w = wy + w, die Gleichung 


qa* Sindy cos? 7 


(Coj + cos 7) w = | is Sinan dee, (| cosA Ed), (11) 


0 


worin fiir @,(7) noch der Ausdruck (8) einzusetzen ware. Soll nun die Durchbiegung w auf den 
heiden Randern 7) =« und 7 =/ verschwinden, so miissen zu gleicher Zeit die beiden Gleichungen 
q at Sin A n cos? y 

8D Sindasinn 


mn) =0 (12}) 


H=a, 8 
bestehen. 

Um die beiden anderen Randbedingungen in einer passenden Weise zu formulieren, betrachten 
wir zunachst die zwei Identitaten 


i Pos ra) 
grad, w = = (Coj € + cos 7) a 5 | | 
ee : (13) | 
in [(Co} § + cosy) w] = (Co) & + cos 1) 7 w sin 7 . | 


Soll nun langs eines Randes 7 = konst. neben w auch grad, w verschwinden, so erweist sich offen- 
bar die Bedingung grad, w —0, d.h. Qw/dy = 0 als vollig aquivalent mit der Randbedingung 


ci 2 (oj —E + cosy) w] =0. (14) 
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Auf den Ausdruck (11) angewendet, fiihrt die vorstehende Gleichung zu den beiden Randbedin- 
gungen 

q a* AGojAn cos? y sin yn — Sindy cos 7 (1 + sin? n) , dprn) 

8D Sind a sin? vp ee: 8 Re ri 0, (15) 


die zusammen mit den Gleichungen (12) hinreichen, um die vier in den Ausdruck (8) eingehenden 
Integrationskonstanten zu bestimmen. Hiernach bestimmt sich auch die elastische Flache der 
Platte durch den Ausdruck (3) in Verbindung mit den Gleichungen (4) und (9). Die Biegemomente 
- der Platte ergeben sich schlieSlich mittels der bekannten Ausdriicke + 
D ia o ( : 
eae Coe + cosy)? (age Ly =) + (1—v) (Cof & -++ cos) (sin § + siny | , | 


on? 


D At Ge 2 - ie os 
m,——2. (oy £ + cosy)? (» = a 3) — (1—v) (Cj € + cosn) (sin iF + siny mi)| , | 


in denen y die Querdehnungszahl bezeichnet. 


3. Die eingespannte Kreissegmentplatte unter gleichférmiger Belastung. Die wirkliche Be- 
stimmung der Integrationskonstanten und die Berechnung des gréBten Biegemomentes soll nur 
fiir die Kreissegmentplatte als Beispiel durchgefiihrt werden; diese sei durch den Kreisbogen 
7 = mit dem Zentriwinkel 2 6 und durch die Gerade 7 = « = 0 begrenzt (Abb. 2). Das System 
der Gleichungen (12) und (15) liefert in unserem Fall die folgenden Konstanten: 


= gan 

o OD Sindx’ 
a qa* _ AGSinA B oj A B sinB — 2 Gin? 1B cosh + 7? sin? B cosp 

2 8D GinAz sinB Sin? A BP —2? sin? B ; (17) 
C, = —AC,, 
eee qa (Sind 6 cos Bf —A Coj AB sin f)? 

4 8D SinAmsin? B Sint2Ab—2 sin®Bo : 

Die elastische Flache der Platte ist also 
fod 2 2 ae, 1 ) 
Vv 64D i iis bees Goj € + cos n [ ven OES Oke a 


0 


) wobei g,(7) durch die Ausdriicke (8) und (17) 
bestimmt wird. Das gréBte Biegemoment der 
Platte entsteht in der Mitte des geraden 
Randes, also in =0,7=0. Nach Aus- 

_ fiihrung der vorgeschriebenen Differentiation 

an der Funktion (18) ergibt sich das be- 

treffende Moment zu 


q a? 
4 m, = —— 
ive) 


qa | (Sind B cos B —A CojAB sin f)? Add 


(Sin? A B —A? sin? 8) Sind a 
: (19) 


enw 

A 2 : F c=afsink f-atoz 

In Abb. 2 ist dieses Einspannmoment im Ver- fing f-ata7 

er, 5 eco og dere f ne ea. 

aac zurdens Wert gy, mit f—atg (fj), rs cent, carta cingepannten Kreisecementplatte in Ab- 
| und in Abhangigkeit von dem halben Zentri- hangigkeit vom halben Zentriwinkel 6. 


| winkel 6 des Kreisbogens dargestellt. 


2 sin? 7 


4. Die gleichmifig belastete eingespannte Halbkreisplatte. Die. Losung fiir die in Abb. 3 dar- 
gestellte Halbkreisplatte ergibt sich, indem man / = 7/2 setzt, als ein Sonderfall der soeben dis- 
kutierten Lésung. Der Ausdruck (18) nimmt sodann die Form 


ra 
= gp —!— 9! apg tram iH ae (20) 
ee (cine ry —1) Sinda 
coe 
1 Vel. z. B. W. Fliigge, Die strenge Berechnung von Kreisplatten unter Einzellasten, S. 15, Berlin 1928. 
A 
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an, wobei ; 
Oe Sm G e- nj (sin*y cos 1) + Ao] = sin n)—?* Coj A cos 4 (21) 


sein soll. 


Die gréRte Durchbiegung der Platte entsteht an der Stelle x =0, y = 0,486a (bzw. & = 0, | 


7, = 0,905) und ergibt sich zu 0,002021 q a*/D in guter Ubereinstimmung mit den von Hassé und 
auch von Conway und Huang berechneten Werten (0,002022 bzw. 0,0020202 q a’/D). Das groBte 
aller Biegemomente ist wieder das Einspannmoment in x = y = 0 (baw. £ =7 =0). Der von 


Wp allein herriihrende Anteil ist hier gleich (1 + 7) q a?/16. Das Glied » q a?/16 hebt sich durch | 
einen gleichwertigen Integralausdruck fort; ersetzt man nun das restliche Glied g a?/16 durch das _ 


Int egral a ue 
ra) ga 
qa? j' Ada —0,Y85L 
7 | SGindn ’ 


0 
so erhalt man im Endergebnis das EKinspannmoment 


a 3 
it ce - Varn Cees 0,0731 qa’. 
s (cine — i) ere 
(22) 


Auch die Biegemomente langs der Mittellinie der Platte 


-O020 xs 


0 


7-0 


7) tata’) 


Ag=+00o) 2 


—00731 


Abb. 3. Die Halbkreisplatte im Koordinatensystem. Abb. 4. Der Verlauf der tangentialen (m_ = me) al 
4, id 


der radialen (m, = m,) Biegemomente langs der Mittel- 
— a . “il . . . 
(x = 0 baw. € = 0) lassen sich nach den allgemeinen lini ciner eingespannten Halbkreisplatte bei gleich- 


Formeln (16) leicht ermitteln. Auf die Wiedergabe der crac Pelee ee 
betreffenden Ausdriicke sei hier verzichtet. Der Verlauf der unter der Annahme y = 0,3 be- 
rechneten Momente ist in Abb. 4 dargestellt. 


5. Die eingespannte Halbkreisplatte mit hydrostatischer Belastung. Die Halbkreisplatte (Abb. 3) 


trage nun eine nach dem Gesetz 


(pom ey ? (23) 


verteilte Belastung, worin q) eine Konstante ist. Die Lésung der Differentialgleichung A Aw = (qoy)/ 
(a D) moge jetzt, der Annahme (3) entsprechend, in der Form 
A Cresta . 
w= Wy Ww, = [ora Xs | bane toss p(y) cosa E d). (24) 
0 


angesetzt_ werden. Die Durchbiegungen wy verschwinden nunmehr lings der gesamten Rand- 


linie der Platte, ebenso wie der grad, wy auf dem Kreisbogen. Auf dem geraden Rande hat man 
dagegen 


Ow, (a? == x)? a’ 
Sy Woiyy. f= 0 _ do < q0 ; 
grad, Wy & ir . 192 aD 18 D (Goi £1)? (25) 
oder, wie dies aus den Bezichungen (13) folgt, 
a) 4 1 Pe 
Bat ((s¢ = qo @ qo @ (A? + 1)A cosa Edd 
on [(Coj é je cos 1) Wo| 48 D (Gof é — 1)2 19 D | ae S nia a ert) (26) 
F 


wenn man die Funktion (of € +1)? zum SchluB noch durch ein gleichwertiges Integral er- 


setzt. 
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Die Randbedingung grad, w = 0 oder, was dasselbe ist, 
a) < 
oy [(Co} & + cos 17) (wy + w,)] = 0 (27) 


auf 7 = 0 1aBt sich nun auf Grund der Gleichungen (24) und (26) in der Form 


qo at (4? + 1)A dein) 
72D Gindn dy ee y ce?) 
schreiben. Die drei weiteren Randbedingungen sind offenbar 
ac } dpi(n) 
eet) pcs EPI) 
[vAMn-0= 9, [pan], 2 =9, eee (29) 


Die Konstanten C,, C,, C3, Cy im Ausdruck (8) fiir die Funktion ~,(7/) lassen sich aus dem System 
der Gleichungen (28), (29) sofort bestimmen. Der Ausdruck (24) formt sich dann um in 


io.@) 


nee doy Gy)? WM at Q, (22 + l)Acosd & da 30) 
192 aD 72 D (Goj E + cos 7) Ah Dae ee eu, 
(ci a beg 72) Gin Aa 
. W 
mint 
celle Laas 
2Q5=GSin G ee i} Sin = sin 77 —/A GitAy cosy . (31) 


Das absolut gréBte Biegemoment der Platte tritt jedenfalls an einem Ende ihrer Mittellinie x = 0 
auf. Fir das Einspannmoment in der Mitte des geraden Randes (y = 0 baw. 7 = 0) findet 


man den Wert 


(oe) 


2 2 2 
Mh = a “ oe S 0,0276 qo a® (32) 
Zin? — 2 
und fiir das Kinspannmoment in der Mitte des Halbkreisbogens (y =a bzw. 7 = >) den 
Wert i 
ioe) 
quet) igyat (a2 4+ 1) a8 dd : 
ea a “od ! 3 Co. saline we ae Bae ios do. ae Be) 
Sin? > 7) Sin > 


0 


Fir die zugehérigen Momente m; gilt die einfache Formel m: = m,. 


6. Die zwischen zwei zueinander exzentrischen Kreislinien eingespannte Platte. In grund- 
satzlich ahnlicher Weise wie das Bogenzweieck aBt sich das in Abb. 5 dargestellte ringférmige 


' Gebiet, etwa bei gleichformiger Belastung, behandeln. Sind die Halbmesser der beiden Kreise 


L 


_rander in Bipolarkoordinaten nach den Formeln 


und die Exzentrizitat gegeben, so erhalt man den zugehérigen halben Polabstand a, die Lage 
der Pole zu den Kreisen und die Gleichungen der Platten- 


a 


pe are oe 


Ae? 
xq = ya? 4. d2, (34) a 
0) oe Sie 6 2 See earn 2 
Sy = Ww om—, §=0 an Sn re 
Im Ansatz (3) fiir die Durchbiegung der Platte moge wy 
wieder von der Form (4) bzw. (5) sein, und die allgemeine 
Lésung der homogenen Differentialgleichung (6) in der Form ee 
eimy y 
im “Coj E+ COS 7) pmlS) 8) Abb.5. Die Kreisringplatte im Koordinatensystem. 


gewahlt werden, worin unter m eine ganze Zahl zu verstehen ist. Setzt man diesen Ausdruck in 
die Gleichung (6) ein, so ergibt sich fiir m + 0 
@m(§) = Cy Gof (m +1) & + C, Gin (m + 1) E + Cy} (m— 1) + C, Sin (m—1)é& (36) 
4 * 


Ingenieur-Archiv | 
[ 


52 Woinowsky- Krieger: Uber die Verwendung von Bipolarkoordinaten 


und fiir m = 0 . | 
go(&) = C, Gof & + Cy Gin & + Cy ECof E+ C,E Sn E. (37) | 
Die Integrationskonstanten C,, C,, Cz, C, dieser Ausdriicke sind wiederum abhangig von der. 


jeweiligen Zahl m. 
Da die Funktion w,(y) nicht nur periodisch, sondern mit Riicksicht auf die angenommene 


Lastverteilung auch symmetrisch zu dem Durchmesser 7 = 0 sein muf, so bleibt sie auf die 
Kosinusglieder in 7 beschrankt, d.h. man hat 


(Coj & + cos 7) wy = 2 mle) cos m1 (38) 
anzusetzen. Entwickelt man die Pain ros (5) in der analogen Weise, so findet man 
¥ ; =—_ q at Co}? g 1 2 Wh pts 
(Gof € 4 cos 7) wy = 8D Ciné ie | ay 1)" e cosmy|. (39) 


Es bleibt jetzt nur noch iibrig, die zu jedem m gehérigen vier Konstanten C; in den Ausdriicken 
(36) und (37) aus den Randbedingungen 


(Cof € + cos 7) (wy + w,) = 0, (4 
2. [(Go} & + e081) (uy + 1)] =0 (41) 


auf € =y und  =6 zu bestimmen, indem man mit den Lésungen (38) und (39) in diese Bedin- | 
gungen eingeht. Es ist hierbei zu bemerken, da& wegen der Randbedingung (40) die Bedingung | 
(41) gleichwertig mit grad. w= 0 wird; dies ist aus den Gleichungen (13) sofort ersichtlich, | 
wenn man das Differentiieren dort nach £ anstatt nach ansetzt. Sind schlieBlich die Konstanten C; | 
ermittelt, so ist auch die elastische Flache der Platte durch den Ausdruck 


q 1 = | 
re ee tae re ee Dai cos m7] (42) | 


m= 


endgiltig bestimmt. Im Falle einer hinsichtlich der x-Achse antimetrisch, etwa nach dem Ge- 
setz p = (q ¥)/a, verteilten Belastung ware die Kosinus-Reihe in Ansatz (38) durch eine Sinus- 
Reihe zu ersetzen, im itbrigen wiirde sich aber an dem ganzen Vorgang nichts andern. 


(Eingegangen am 20. Juni 1955.) 


Anschrift des Verfassers: Professor Dr.-Ing. S. Woinowsky-Krieger, Quebec 6, Que, (Canada) 200 
Ave. Marguerite Bourgeoys Apt. 8. 
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Das Stabilitatsverhalten der nichtlinearen Biegeschwingungen 
des axial pulsierend belasteten Stabes. 


Von F. Weidenhammer. 


1. Einleitung. Der axial pulsierend belastete Stab ist das bevorzugte Modell fiir zahlreiche 
Untersuchungen ! des Stabilitatsverhaltens elastischer Kérper unter schwingenden auBeren Lasten. 
Dabei untersucht man mit der Methode der kleinen Schwingungen die Stabilitat der gestreckten 
Lage des Stabes, abhangig von der Amplitude der schwingenden Langslast und deren Frequenz. 
Bei technischen Fragestellungen liegt dabei der pulsicrende Lastanteil weit unterhalb der Knick- 
last und auch Langsschwingungen merklicher Amplitude werden im Stab nicht erregt. Unter diesen 

_Voraussetzungen sind diejenigen Frequenzen leicht anzugeben, bei welchen der Stab in seiner ge- 
streckten Lage instabil in bezug auf eine Querauslenkung wird. Die meisten Untersuchungen be- 
schranken sich dabei auf beidseitig momentenfrei und quer zur Stabachse unverschieblich gelagerte 
Stabe, und dieser Lagerungsfall soll auch in der vorliegenden Arbeit ausschlieBlich betrachtet wer- 


den. Ein Lager sei ferner in Langsrichtung ver- 42 
schieblich (Abb. 1) und die Kreisfrequenz « der pul- _|% 
(x,t) 
RE se) 
Abb. 1. Der Stab im Koordinaten- Abb. 2. Instabilitatsbereiche bis zur 4. Ordnung. Ordinate: 
system, Verhaltnis des pulsierenden Lingslastanteiles zur Knicklast. 


Verkleinerung der Bereiche durch Dampfung: — — — . 


sierenden Axiallast bleibe weit unterhalb der tiefsten Eigenfrequenz der Langsschwingungen. Dann 
findet man durch cine Stabilitatsbetrachtung, daf die gestreckte Lage instabil wird, wenn die 
Langserregungsfrequenzen in der Nahe der speziellen Werte 


2; 
Oy => (P=1.2,--,) (1) 


liegen?. In (1) bezeichnet w; (j =1, 2,...) jeweils eine Eigenfrequenz der Biegeschwingungen des 
Stabes, und es gibt zu jeder Eigenfrequenz gemaB den Ordnungszahlen p eine Folge von Werten «,, 
die jeder zu einem ganzen Bereich von Werten der Frequenz ~ fihren, welche Instabilitat der 
gestreckten Stablage zur Folge haben (Abb. 2). Auf Grund der groBeren Frequenzbreite der Be- 
reiche erster Ordnung (p = 1) und der hier am wenigsten stabilisierend wirkenden Dampfungs- 
krafte sind diese Bereiche am wichtigsten, wie man aus Abb. 2 entnehmen kann. Im weiteren soll 
daher nur das Schwingungsverhalten des Stabes in der Nahe der Frequenzen wm) = 2 w; genauer 
untersucht werden. Bei diesen Langserregungsfrequenzen ist die gestreckte Stablage nicht stabil, 
und der Stab wird nach einer anfanglichen, kleinen Querauslenkung in Querschwingungen end- 
licher Amplitude iibergehen. Die zugehérige Frequenz-Amplitudengleichung der dann eintretenden 
nichtlinearen Biegeschwingungen ist fiir den ungedampft schwingenden Stab bekannt* und ist in 
Abb. 3 fiir einen konstanten Wert der Langserregungsamplitude als Resonanzkurve aufgetragen. 
Als Amplitude ist dabei der Fourierkoeffizient tiberwiegender Grofe der nichtlinearen Schwingung 
bezeichnet und diese Amplitude wird in ihrer Abhangigkeit von der Frequenz betrachtet. Der 
gestrichelt eingezeichnete Ast der Resonanzkurve hat keine physikalische Bedeutung, da er einer 
instabilen Schwingung zugehért und also nicht beobachtet werden kann. 

Uber die soeben aufgezahlten bekannten Ergebnisse hinausgehend werden in der vorliegenden 
Arbeit folgende Probleme behandelt: 
‘ 0) Einschwingvorgange. Um einen tieferen Einblick in das Stabilitatsverhalten des Stabes 
zu erhalten, werden nicht nur die nach einer Anfangsstérung schlieBlich eintretenden stationaren 
Schwingungen untersucht, sondern es werden auch die Ubergangsschwingungen in diese End- 


1 Vel. das umfangreiche Literaturverzeichnis bei E. Mettler, Publ. sci. techn. Min. Air. Paris 281 (1953), 
Sh lee 

2 Vel. z. B. E. Mettler, Mitt. Forsch. Anst. GHH-Konzern 8 (1940), 5. 1. 

3 F. Weidenhammer, Ing.-Arch. 20 (1952), Se Sl. 


| 
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zustinde hinein berechnet. Dadurch kann man dic bisher offene Frage beantworten, ob ein Stab) 
bei einer Langserregungsfrequenz w > 2 @, in stationare Querschwingungen tibergeht oder in die! 
Ruhelage zuriickschwingt. Die bisherigen Untersuchungen haben lediglich fiir diese Langserregungs-~ : 
frequenz sowohl die Ruhelage (Amplitude A = 0 in Abb. 3) als auch Querschwingungen (A =: 0) 
in Abb. 3) als mégliche stabile Bewegungsformen erkannt, konnten jedoch die nun gefundene» 
Abhingigkeit von der Anfangsstérung nicht aufdecken. Die in diesen Fragenkreis neu ein-- 
gefiihrte Lésungsmethode der ,,langsam veranderlichen Phase und Amplitude‘ erweist sich da--|| 
bei als sehr zweckmabig!. V 

6) Einbeziehung dimpfender Krafte. Auch bei der Berechnung gedampfter Schwin-} 
gungen bewihrt sich diese Methode, denn der bei gedampften Schwingungen schwierige Stabilitats--) 
nachweis kann damit leicht iiber die H. Poincarésche Singularitatentheorie gefiihrt werden. So) 
findet man, daB dampfende Krafte eine ganz im Endlichen verlaufende Resonanzkurve mit einer 3 
endlichen Gré®tamplitude bedingen, wie dies von.) 
erzwungenen Schwingungen mit Dampfung her be-) 
kannt ist. 

y) Berechtigung der Naherungsrechnung., 
Da das Schwingungsverhalten des Stabes durch kom- - 
plizierte nichtlineare partielle Differentialgleichungen | 
beschrieben wird, ist man auf eine Naherungsrechnung ; 

i F angewiesen. In dieser fiihrt man iiblicherweise die Dis- | 

Tees So ™=SC(«é‘iéi@«‘ssssionn «es Schwingumngsverhaltens des Stabes an) 

J ° : : . “1s } 

ANR-S. Rosas cuce eierdoa Unreana pe cekeingenacneene noe Schwinger mit nur einem Freiheitsgrad durch. | 
Stab. Stabile Aste: = instabile Aste: --- —. Es ist daher wesentlich, nachzuweisen, da sich die | 
bisher gefundenen Schwingungseigenschaften auch in | 

einer sehr verbesserten Naherungsrechnung wiederfinden, in der man den Stab als System mit | 
mehreren, jedoch endlich vielen Freiheitsgraden behandelt. 

Die Untersuchung beschrankt sich auf kleine Nichtlinearitaten und kleine pulsierende Langs- | 
lasten mit Frequenzen weit unterhalb der Langsresonanzen. Auch wird stets mit einem linearen | 
Hookeschen Werkstoffgesetz gerechnet, so daB Nichtlinearitaten nur auf Grund der Kopplung von | 
Langs- und Querverschiebungen in die Rechnung eingehen. Die genannten Einschrankungen diirf- 
ten jedoch bei einem Stab als Bauglied in der Technik von selbst gegeben sein. 

Bei einem Stab mit komplizierterer Lagerung (z. B. eingespannten Enden) treten Instabili- 
taten der gestreckten Lage zusiatzlich bei den Frequenzen w,) = (w; + @,)/p (7 =1,2,...: 
k=1,2,...;j +k; p=1,2,...) auf. Derartige Instabilitaten zweiter Art kénnen jedoch im 
Falle momentenfreier Lagerung der Stabenden nicht auftreten, so daB die vorliegende Untersuchung 
zur Lésung des schwierigen Problems? der Instabilitaten zweiter Art keinen Beitrag geben kann. 


Amplitude 


2. Ausgangsgleichungen. Die Bewegungsgleichungen des Stabes lassen sich zweckmafig aus 
dem Hamiltonschen Prinzip herleiten: 


ty ty 
06 f(A: + V—T)dt+ ([bAdt=0, (2) 
t fy | 
worin A; die Formanderungsarbeit, V das Potential der auBeren Krafte, T die Bewegungsenergie 
und 0A die virtuelle Arbeit der geschwindigkeitsproportionalen Dampfungskrafte bezeichnen. Mit 
den Annahmen der technischen Biegelehre ist?: 


1/2 
EF Ne E. 
A= | G i i dx += | wi, dx, V = P(t)u(l2), | | 
—ij2 1/2 | | 
12 1/2 ) 
T= 5 [ (ui + wi) dx, 64 = | Bry dw de. | 
Sv —1/2 J 


Hierin bezeichnen x die Ortskoordinate, u(x, t), w(x, t) die gesuchten Verschiebungskomponenten 
in Stablangsrichtung und quer zur Stabachse, t die Zeit, J — konst. das fiir die Biegeschwingung 


1 Vel. fiir die entsprechende Fragestellung bei erzwungenen Schwi Ch l 1 
lations in non-linear Systems, S. 81, Osaka (lanan 1953. Tene ae 


2 Vel. E. Mettler, Ing.-Arch., 23 (1955), S. 354. 
° E. Mettler, Ing.-Arch. 13 (1942), S. 97, 


XXIV. Band 1956 Weidenhammer: Das Stabilitatsverhalten der nichtlinearen Biegeschwingungen 55 


maBgebende Flachentragheitsmoment, F = konst. die Querschnittfliche des Stabes, [) = konst. 
die Stabmasse je Langeneinheit, | die Stablinge, E den Elastizitatsmodul, P(t) = Py) + P, coswt 
die rein harmonisch pulsierende Axiallast, so daB P(t) > 0 Langsdruck auf den Stab bezeichnet, 
6 die Dampfungskonstante der Biegeschwingungen. 

Durch Ausfiihrung der in (2) verlangten Variationen erhalt man fiir die Verschiebungen u und 
w ein gekoppeltes System nichtlinearer partieller Differentialgleichungen mit zeitabhingigen 
Randbedingungen, deren strenge Lésung kaum méglich erscheint. Zum naherungsweisen Studium 
_des Zeitverhaltens werden daher die ,,gemischten Ritzansitze‘ 


rs = 1 
ut) =Y8 (+5) tnee Ul), (4) 
Ma ns l 
w(E,t) = Wax >) V2 sin j x (: ate 4] W;(z) (5) 
in 


mit vorgegebenen Ortsfunktionen und aus dem Variationsproblem zu bestimmenden Zeitfunktionen 
U(r), W(t) in (2) eingefiihrt. 

Dabei bezeichnen u,,,, und wy, die Betrige der gré8ten vorkommenden Verschiebungen, so 
daB im folgenden die GréBenordnungen der Nichtlinearitaten erkennbar werden. In den Ansatzen 
(4) und (5) ist die Langsverschiebung u als nicht wesentlich interessierend nur grob eingliedrig an- 
genahert, wahrend fiir die Querverschiebung w ein N-gliedriger Ansatz mit den alle Randbedin- 
gungen erfiillenden Eigenfunktionen des Stabes eingefiihrt wurde. Ferner ist die Zeitzahlung und 
die Ortskoordinate durch 


ta Ode —— > 
dimensionslos gemacht. Mit den Ansatzen (4) und (5) ist (2) nur noch eine Variationsaufgabe fiir 


die Zeitfunktionen U und W; allein und als zugehérige Eulersche Gleichungen findet man nach 
einiger Rechnung! 


9 Umax ce l Umax ) 2 Wmax ‘ 9 V3 12 
ce UF 4-3 2, OU sey ( sue) Doyen Wi = — PU) “py 5 (6) 
2 shin 1 /w a= 
4 W581 Wj + 0G Wi + V3 RO Uy Wy + (Pt) D>) tema Wi Wm Wn = 0 
¥ oy ‘ ais k,m,n=1 
(heed tN) (7) 
Hierin bezeichnen: 
o 
1 
VEJ/u 4 


die Langsanregungsfrequenz in einem dimensionslosen Mabf, 


Pe : 
4a pry et P, G2 ee 


mit der Eulerschen Knicklast Py, = E Jx?/I? die Biegeeigenfrequenzen des ungedampft schwingenden 
Stabes im gleichen Mab, 


ae ere 
‘2px /EBT[u 


die mit den ungedampften Eigenfrequenzen berechneten Dampfungsdekremente und 


ye 
Ur — FF 


den Flachentragheitsradius, wahrend Punkte nun fiir Ableitungen nach der dimensionslosen Zeit 
t stehen. Ferner kommen in (6) und (7) die Zahlenwerte 


Wea fu, wp =jt (f=, 2,...N) (8) 
und 
1/2 ! l\ .. . 
Vita ae f TT p eos px(é +5} 4 RU C= Uae eee kB) 


ib; p=j,k,m,n 


1 Vel. S. 317—320 der in FuBnote 3 auf S. 53 genannten Arbeit. 
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vor. Offenbar ist 
ine 
Nii — 4 mt ae und Xjkm r= 0 9 (9) 
wenn nicht mindestens zweimal zwei Zeigerwerte in beliebiger Reihenfolge iibereinstimmen. Ferner | 
folgt aus (8) und (9): 
4 
, pees ie 
Vivi Hl es Se 
was spater gebraucht werden wird. Da die Langserregungsfrequenz in (6) voraussetzungsgemaB 


weit unterhalb der ersten Langseigenfrequenz (x <x, ~ /3 l/ip) bleiben soll, kann man das Trag- 
heitsglied in der Gleichung fiir die Lingsschwingungen naherungsweise vernachlassigen. Dann 
kann aus (6) und (7) U(z) ganz eliminiert werden und man erhalt ein System von homogenen, nicht- 
linearen Schwingungsgleichungen fiir die Zeitfunktionen W,(7) allein: 


N 
Ty da Ip + Ey;; W; be 2 ya)—0 (PH 12.0427 es) 
+j 


worin die Abkiirziungen 


7 ; 
qj Ri Bj a7 a Oh Pa < 1 (11, 1—3) 
J 
7 Wij 
j= (11, 4) 
eee hin Vie 
VN sa ta Rae (ih) 
J $ 
Wax 234-51 — Ye Vuk il 
Vik =| ip ) ies 7 Fag oe 


verwendet werden. In (10) sind kleine Glieder durch den kleinen Parameter ¢ = P,/Py kenntlich 
gemacht, denn der pulsierende Lastanteil P, wird in der Technik stets sehr klein im Vergleich mit 
der Eulerschen Knicklast sein, so daB ¢e < list. Wahrend die Vorfaktoren der nichtlinearen Glie- 
der klein sind wegen der iiberhaupt nur zugelassenen kleinen Maximalamplituden (wWynax/ir< 1), 
ist die Dampfung meist nur durch die sehr kleine Werkstoffdampfung gegeben, so daB das zu- 
gehérige Dekrement ebenfalls klein ist: 

igre > 2 pra. 

Die Gleichungen (10) sind die grundlegenden Schwingungsgleichungen, aus denen heraus das 
Schwingungsverhalten des Stabes diskutiert werden soll. Die Besonderheit der zu untersuchenden 
Klasse von Schwingungen kommt in den Gleichungen (10) formal dadurch zum Ausdruck, da® 
die Erregung nur tiber periodisch veranderliche Koeffizienten in die Schwingungsgleichungen 
cingeht und da diese Gleichungen homogen sind und demnach stets auch die Lésungen 
W, =0(j =1,2,..., N) zulassen. Dieser Schwingungstyp wurde von K. Klotter! rheonichtlinear 
genannt, wahrend sich neuerdings die von N. Minorsky* gewahlte Bezeichnung dieser Schwingungen 
als parametererregt durchzusetzen scheint. Jedenfalls handelt es sich um einen speziellen Typ er- 
zwungener Schwingungen, der sich durch sein ganz andersartiges Resonanzverhalten von den viel 
behandelten erzwungenen Schwingungen des Duffingschen Problemes® unterscheidet. 


3. Einfiihrung der Methode der Jangsam veranderlichen Phase und Amplitude. Bei der Unter- 
suchung der stationdéren Schwingungen in der Nahe einer Langserregungsfrequenz Wy = 2 w; kann 
man sich auf die Diskussion der j-ten Schwingungsgleichung des Systems (10) allein beschranken, 
wie in Abschn. 6 gezeigt werden wird. Die Einschwingvorginge sollen daher ebenfalls aus dieser 
Kinzelgleichung allein berechnet werden. Daher wird in (10) fiir einen festen Zeigerwert j (etwa 
j=1 fiir die Grundschwingung) W;(r) == 0, W,=0(k=1,2,..., N, == J) gesetzt, und die dann 


K. Klotter, Techn. Schwingungslehre, Bd. 1, S. 350, 2. Aufl., Berlin 1951. 
N. Minorsky, Introduction to non-linear Mechanics, S. 355, Ann Arbor 1947. 
Vgl. z. B. S. 337 des in FuBnote 1 genannten Buches von K. Klotter. 
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verbleibende Schwingungsgleichung unter Fortlassung des Index j einfacher 
PW +ebBqw 4+ W(l—eacost +ey W2) =0 (12) 

geschrieben. Aus dieser nichtlinear erginzten Mathieuschen Differentialgleichung heraus sollen 
die in Abschn. 1 gestellten Fragen x) und £) beantwortet werden, wozu die Methode der langsam 
verdnderlichen Phase und Amplitude in der von B. v. d. Pol! angegebenen Form verwendet wird. 
Diese Methode gestattet fiir Kleine Parameter ¢ die durch (12) beschriebenen periodischen Schwin- 
gungen naherungsweise zu berechnen. Dariiberhinaus vermag sie auch einen Einblick in jene fast- 
periodischen Schwingungsbewegungen zu geben, die in unmittelbarer Nachbarschaft der statio- 
niren, periodischen Schwingungen verlaufen. Die Methode ist daher geeignet, zu entscheiden, 
welche der beiden stabilen Bewegungsformen der Stab annimmt, wenn fiir Frequenzen w > 2 1; 
sowohl die (periodischen) Bewegungen W = 0 als auch W == 0 méglich sind. 

Da fiir die zu untersuchenden Frequenzen q ~ 2 ist, wird zur naherungsweisen Lésung von (12) 
W(r) mit cos 7/2- und sin t/2-Funktionen von der halben Frequenz der Erregung in der Form 


W(t) = A(r) cos > + B(x) sin + (13) 


angesetzt, worin A(t) und B(t) noch aus (12) zu bestimmende Funktionen bezeichnen. Da diese 
Funktionen durch den Ansatz (13) und die eine Differentialgleichung (12) noch nicht vollstandig 
festgelegt sind, wird die Differentialgleichung 

Tt 


. T Cae 
A cos + B sin > 


= 0 (14) 


als Nebenbedingung zur Ausgangsdifferentialgleichung (12) hinzugefiigt. Dadurch erhalt W’ 
einfach die Form 


und wird damit frei von Summanden A’ und B’, die erst in der Beschleunigung W" vorkommen. 
Somit liefert die Nebenbedingung (14) wesentlich die Spezialisierung auf Schwingungen mit lang- 
sam veranderlicher Phase und Amplitude. Weiter wird die vorausgesetzte Kleinheit der Damp- 
fungskraft, der Langserregung und der Nichtlinearitat durch die Umschreibung der Differential- 
gleichung (12) in die Form 


g@w +W=«(—qpBW-+acost W—y W?) (15) 
mit dem Gréfenordnungsfaktor ¢ sichthar gemacht. Tragt man nun den Ansatz (13) in (15) ein, 


so kommen rechtzeitig auf Grund der Nebenbedingung (14) keine Zeitableitungen von A und B 
mehr vor und man erhilt fiir die gesuchten Funktionen A und B die Differentialgleichung 


Gh. in eee GE coe Ts. q° T ¢ at 
sat sin ~ 4 5 B cos; =A(G—leorg + B(G—1sing +e F(A, Baa), (16) 


worin F die Funktion 


3 
F(A, B,r) =—98(—4 sin Lae cos 7) cost (4 cos 5+ B sin) —(Aeosy +Bsin->} 


yezeichnet. Da das Frequenzverhaltnis q ~2 ist, darf fiir q?/4—1=0(e) vorausgesetzt werden, 


so daB alle rechtzeitig vorkommenden Summanden die GréBenordnung ¢ nicht iibersteigen. Durch 


Auflésung der Gleichungen (14) und (16) nach A’ und B’ ergibt sich nun 

Gay ts q re ey We eens te 

) 34 =—A(F—1Jeos Fein —B 1) sin 5 e F(A, B, t) sin, Re, 
fe ( 

FB = A(T — loos + BG —I)sin 5 cosy +e F(A, By) cov 5. | 


Dieses nichtlineare Differentialgleichungssystem (17) ist zufolge seiner expliziten Abhangigkeit von 
ler Zeit schwer zuganglich. Es ist daher wesentlich, da8 man stationare und auch Ubergangs- 
chwingungen aus dem zugeordneten zeitfreien System naherungsweise berechnen kann. Zum 
Jbergang auf das zeitfreie System werden die Gleichungen (17) rechtseitig iiber eine Periode 4 7 
-emittelt und dabei 4, B wie Konstante behandelt. Anschaulich ist durchaus verstandlich, da8 man 


1 Vel. z. B. A. A, Andronow u. C. E. Chaikin, Theory of Oscillations, S. 302, Princeton 1949. 
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4. 


auf diese Weise zu einer brauchbaren Naherung gelangen kann, da ja schwachveranderliche A(t) 
und B(r) berechnet werden sollen und diese Funktionen in den kleinen Summanden der rechten 
Seiten ohne wesentlichen Fehler fiir die Dauer einer Periode als konstant angesehen werden kénnen.. 
Fiir diese Mittelbildung hat man die in 


An 


2 2 2 : ee | 
i Ages 2 | E A (7 — 1) cos sin _ — BG — 1) sin? 7a e F(A, B, t) sin oe Wie 

A 7% 1 
ab Wf A AGE one Ch Sa ag Tw 
ee =i [4 (7 — 1) 00s x + B(G—1)sin 5 cos 3 + ¢ F(A, B,t) cos 5 dt | 


auftretenden Integrale auszuwerten, was nach kurzer Rechnung auf 


gd =B(I—T tes teqy [4 + Beh 5A, | 


gp B=—eptB—A (1 ee ee B]) | 
fiihrt. Dieses nunmehr zeitfreie System kann als eine nichtlineare Differentialgleichung erster Ord+ 
nung fiir A(B) in der Form | 
dA Biot oer a By) —ep 4 Z (19) 
(appeal es 3 ey [at a a) i | 


behandelt werden. An (19) kniipfen alle weiteren Rechnungen an. Die Untersuchung der Uber 
gangsschwingungen mit Hilfe des einfacheren Ersatzsystemes (18) oder (19) ist unter gewissen Be 
dingungen naherungsweise zulassig, wie J. Haag! nachgewiesen hat. Die Lésungen des urspriing-~ 
lichen zeitabhingigen Systems unterscheiden sich dann von denen des zugeordneten zeitfreie 
Systems bei gleichen Anfangsbedingungen nur um einen Betrag von der GréBenordnung e. _ 
Ersatz der Gleichungen (17) durch (18) und (19) mag daher fiir die vorliegende Aufgabenstellung al 
ausreichende Naherung gelten. 


Auch stationare periodische Schwingungen lassen sich einfach aus (19) berechnen. Denn fii 
stationare Schwingungen mu zufolge des Ansatzes (13) A = konst., B = konst. eine Lésung deni 
Gleichungen (18) sein, was offenbar nur méglich ist, wenn in (19) rechtseitig Zahler und Nenner zu 
gleich verschwinden. Da dann eine singulare Stelle der Differentialgleichung (19) vorliegt, kann ma» 
hier die H. Poincarésche Theorie der Singularitaten heranziehen, um iiber die Stabilitat der Schwin4: 
gungen zu entscheiden. | 


4, Stationire Schwingungen und deren Stabilitat. Obgleich das eigentliche Ziel der Unter- 
suchung die Berechnung der Einschwingvorgange ist, miissen zunachst die stationaren Endzustandey 
der Stabbewegung berechnet werden. Dies ist notwendig, weil zwar fiir den ungedampft schwingen-} 
den Stab diese stationaren Bewegungen bekannt sind, aber die entsprechenden Angaben fiir ge+} 
dampfte Stabschwingungen noch fehlen. 


a) Ungedampfte Schwingungen. Der Vollstandigkeit halber werden zunichst bekanntey 
Ergebnisse” fiir den ungedémpft schwingenden Stab mit Hilfe der in Abschn. 3 ee 
Methode nochmals kurz hergeleitet und vorangestellt. Zur Berechnung der ungedampften| 
Schwingungen hat man in (19) 6 = 0 zu setzen und diejenigen speziellen festen Werte von A und BI 
zu bestimmen, fiir die (19) singular wird. Dies fiihrt auf die Koeffizientengleichungen 


ie 3 
B, ees tet testy [Ai + B3}) =0, 


AS -e$ tebyta+ ay) =0 


J. Haag, Bull. Sci. Math. 72 (1948), S. 69. 


_® Ygl. 5.325 der in FuBnote 1 auf Seite 53 genannten Arbeit. Die dortigen Gleichungen (57), (59) stimmen 
mit (20,2), (20,3) tiberein, wenn man nur lineare Potenzen von « beibehilt. 
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mit den Lésungsméglichkeiten 


1) An == 0, Bo = 0. Wane (20,1) 
- qg° ee 8 x . a 
2) eae tere oy ey Ai = 0, W = + Ay cos, (20,2) 
2 3 : 2 
3) Ae Os 1 —t—e ¥ yey Be 0s W = + B, sin 5 ; (20,3) 


Wie bekannt gibt es also auSer der gestreckten Stablage ohne Querschwingungen (Fall 1), je 
einen cos- und einen sin-Typ stationarer Schwingungen mit den angegebenen Frequenz-Ampli- 
tuden-Abhangigkeiten (Fall 2 und 3). Da jedoch der sin-Typ nicht stabil ist!, kommt dieser Lo- 
sung keine physikalische Bedeutung zu. Hingegen kénnen der cos-Schwingungstyp und die Ruhe- 
lage unter gewissen Bedingungen sogar gleichzeitig stabil und damit beide méglich sein. Die Auf- 
findung dieser Stabilitatsbedingungen wiirde mindestens eine Betrachtung der Lésungen von (19) 
in der Nachbarschaft der ausgezeichneten Werte 4), B, erfordern. Hier kann man jedoch die L6- 
sungen in der ganzen (A, B)-Ebene ohne Einschrankungen geschlossen angeben und gewinnt durch 
eine derartige Betrachtung der Ubergangsschwingungen auch ein besscres Verstandnis der Stabili- 
tatseigenschaften der stationéren Schwingungen. Die Stabilititseigenschaften der Lisungsmég- 
lichkeiten (20,1—3) werden daher bei der 
Berechnung der Ubergangsschwingungen 
(Abschn. 5a) von selbst erkennbar und brau- 
chen daher hier nicht abgeleitet werden. 


VASB6 


b) Gedimpfte Schwingungen. Die 
in Abb. 3 fiir den ungedampft schwingenden 
Stab wiedergegebene Resonanzkurve der 
Querschwingungen kann insofern den 
physikalischen Sachverhalt nicht voll be- 
friedigend wiedergeben, als daraus keine 
endliche GréStamplitude fiir w > 2; zu 
entnehmen ist. Da beim Duffingschen Pro- 
blem der erzwungenen nichtlinearenSchwin- 
gungen fiir m>2.«j; eine derartige Be- 


grenzung der Amplitude besteht, wird man Abb. 4, Resonanzkurven fiir den gedimpft schwingenden Stab, 

i ‘ : : Stabile Aste: , instabile Aste: — — —. Dampfungskonstanten: 
auch ee vorliegenden Falle eine pone B,>++:> B,. Resonanzkurve fiir ungedimpfte Schwingungen (6, = 
Endlichen verlaufende Resonanzkurve ver- zum Vergleich: —-—- = 


muten. Wie sich zeigen wird, geniigt auch 

hier die Erfassung der dampfenden Krafte durch einen einfachen geschwindigkeitsproportionalen 
Ansatz. Die Berechnung der stationaéren Biegeschwingungen bei Anwesenheit dampfender Krafte 
kniipft zweckmaBig an die zeitfreie Differentialgleichung (19) und nicht unmittelbar an die Schwin- 
gungsgleichung (12) an, denn das Auftreten einer ersten Ableitung nach der Zeit infolge der Damp- 
fungskraft macht die Berechnung periodischer Schwingungen mit der iiblicherweise verwendeten 
Stérungsrechnung unméglich. Da jedoch stationére Lésungen aus dem zeitfreien System nahe- 
rungsweise berechnet werden kénnen, ist dieser Weg besonders auch im Hinblick auf den anschlie- 
Bend notwendigen Stabilitatsnachweis zweckmifig. Entsprechend dem Vorgehen in Abschn. 4a 
hat man auch hier die Singularitaten von (19) aufzusuchen, wenn man periodische Lisungen be- 
rechnen will. An Stelle von (20) ergibt sich zwar fiir die Ruhelage ohne Querschwingungen noch 
A, = B, = 0 aber fiir Querschwingungen selbst das gekoppelte nichtlineare Gleichungssystem 
fur A, und B,: 


2 3 f 
By (1 t beg tp elss + Bil)— eB y= 0, (21) 
Z L 3 : ‘ 
== GY0) 4 By Ay(1 i es Aceh [42 aie B}) —— 0 i (22) 


Fir A, + 0, By = 0 lassen sich die Gleichungen (21) und (22) leicht in bezug auf das »Ampli- 
tudenquadrat** [42 + B?] hin auswerten. So ergibt sich als Frequenz—Amplitudengleichung mit 


1 Vel. die in FuBnote 1 auf Seite 53 genannte Arbeit. 
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dem Frequenzparameter q =x/x; die bereits von N. W. McLachlan" angegebene Naherung 


aa | 


3 Re Ok 2 ped | 
ey [A? -- BY? — i fl ale ee eg? ibe ai ’ (23) | 


worin das Doppelvorzeichen den beiden zusammenlaufenden Asten der nunmehr geschlossenen und © 
ganz im Endlichen verlaufenden Resonanzkurve entspricht. Der untere Ast der Resonanzkurve 
gehért der negativen Wurzel von (23) zu. Er wird sich sogleich als zu einer nicht stabilen Schwin- | 
gung gehérend erweisen und ist daher auch in der Auftragung eines Zahlenbeispieles der Abb. 4 
nur gestrichelt eingezeichnet. Alle Kurven liegen ganz innerhalb der zum Vergleich strichpunktiert 
eingetragenen Resonanzkurve fiir den ungedampft schwingenden Stab. Ubrigens verlangt die 
Realitit der in (23) vorkommenden Wurzel, daB die dort durch ex eingehende Erregung einen ge- 
wissen als Schwellwert bekannten Wert tberschreitet. 

Die Beantwortung der Stabilitatsfrage wird ebenfalls an (19) angeschlossen, indem die auf Sta- 
bilitat zu priifenden Schwingungen durch die Funktionen A(r) und B(r) in der Nahe der aus- 
gezeichneten Werte A, und B, untersucht werden. Zu dem Zweck wird in (19) 

A(t) = Ay + 9(t), B(x) = By + C(t) | 
fiir kleine Werte von 7; und ¢ gesetzt und damit die Natur der Singularitaten der Punkte A,, By | 
untersucht. Eine solche Beschrankung auf differentiell kleine Nachbarwerte ist hier im Gegen- _ 
satz zum ungeddmpft schwingenden Stab notwendig, da die Differentialgleichung (19) fir 6 + 0 | 
nicht mehr geschlossen lésbar ist und daher A und B auch nicht allgemein formelmabig in der || 
ganzen (A, B)-Ebene diskutiert werden kénnen. Zur Untersuchung der Stabilitat der durch 7 | 
und € gegebenen Schwingungen kann man jedoch in bezug auf 7 und ¢ linearisieren und erhalt | 
dann die leicht auswertbare Ahnlichkeitsdifferentialgleichung 


dy _al+by (24) | 
dt cC+dn 


mit den Koeffizienten: 


2 3 5 ’ e Bae 
a=1—L+ep + yey [t+ 3B], b= Sey 4, B—ep st, 


if 2 
= ey 4B) ep i= (2 = eg + ger 34s + Bil). 


Je nach dem Verhalten der Lésungen von (24) kann man hier die Singularitaten A), B, als Strudel- 
punkte, Sattelpunkte oder Knotenpunkte klassifizieren und damit die Stabilitatsfrage beant- 
worten. Da alle in der Nahe eines Strudelpunktes verlaufenden Lésungskurven asymptotisch in 
diesen singularen Punkt einmiinden, gehéren diese Koordinaten Ay, By einer stabilen Schwingung zu. 
Im Gegensatz hierzu fiihrt ein Sattelpunkt auf eine instabile Schwingung. 

Zur Untersuchung der Singularitaéten hat man zundchst die Diskriminante 


D=(b—c)?+4ad (25) 
zu bilden *, Falls dieser Wert D positiv ausfallt, hat man noch 
A=ad—bce (26) 


zur Kntscheidung heranzuziehen. Wegen b +c <0 sind auBer den Strudelpunkten auch alle 
etwa vorhandenen Knotenpunkte stabil. Bildet man zunachst D, so erhalt man nach elementaren, _ 
jedoch langeren, hier unterdriickten Zwischenrechnungen aus (25): 


2 2 3 
D= 40 pet —3cy [At + Ba ( eT [43 + Bil), (27) 


wobei noch kein Gebrauch von der Frequenz—Amplituden-Beziehung (23) gemacht wurde. Daher 
gilt (27) noch fiir beide Aste der Resonanzkurve. 

Die weitere Auswertung soll nun zunichst fiir den unteren Ast der Resonanzkurve geschehen. 
Demgema nimmt D mit der jetzt geltenden Beziehung 


qi 3 gly. Gong 


* N. W. McLachlan, Ordinary non-linear Differential Equations in engineering and physical Sciences, 
S. 126, Oxford 1950. 
2 Vel. z. B. J. J. Stoker, Nonlinear Vibrations, S. 44, New York 1950. 
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einen stets positiven Wert an: 


2 2 2 2 
D,= 42 PPE + 3ey [43+ BY YE — apt so. (29) 


Es kann also nur ein Sattelpunkt oder ein Knotenpunkt vorliegen. Um zwischen diesen beiden 
Méglichkeiten zu entscheiden, muS (26) herangezogen werden. Nach kurzer Rechnung ergibt 
sich unter Verwendung von (28): 


2 py2 2 
A=+3ey [42 + B] ee o>o, (30) 


so daB ein Sattelpunkt vorliegt. Demnach ist der untere Ast der Resonanzkurve stets instabil 
und hat damit keine physikalische Bedeutung. Er ist in der Resonanzkurve daher nur der mathe- 
matischen Vollstandigkeit halber und gestrichelt eingetragen. 

Um die Stabilitét der dem oberen Ast der Resonanzkurve zugehérenden Schwingungen nach- 
| zuweisen, hat man nun die aus (23) folgende Beziehung mit dem positiven Wurzelzeichen in (25) 
einzutragen und erhalt dann an Stelle von (29): 


2 2 oy2 2 

D.=42° $4 —3ey [42 + B?] “+ — 67 B? v. 
Jedenfalls ist D, <0, wenn itiberhaupt keine Dampfung vorhanden ist (fb = 0) und fiir die in 
der Technik nur in Betracht kommenden kleinen Dampfungen diirfte dies auch gelten. Dann 
liegt also ein Strudelpunkt vor. Wird die Dampfungskonstante gréBer gewahlt, so kann D, Null 
oder sogar positiv werden. Man hat dann in beiden Fallen einen stabilen Knoten, da an Stelle 
von (30) jetzt A <0 wird. Damit ist die Stabilitat der dem oberen Ast der Resonanzkurve ent- 
| sprechenden Schwingung allgemein nachgewiesen. 
Anschaulich ist dieses Stabilitatsverhalten durchaus verstandlich, wenn man auf den un- 
'gedampft schwingenden Stab zuriickblickt. Denn der obere Teil der Resonanzkurve gehért zu 
}einem Schwingungstyp, der noch im wesentlichen ein cos-Typ ist, wenn er auch auf Grund der 
| Dampfung eine Phasenverschiebung aufweist. Doch ist in diesem Falle immer noch der Betrag 
) von A, wesentlich gréGer als der von B,. Entsprechendes gilt fiir den unteren Teil der Resonanz- 
/ kurve, die wegen des gréBeren Betrages von B, einem Schwingungstyp zugehért, der im wesent- 
lichen ein sin-Typ ist. 
| 5. Einschwingvorginge. Nach vorangegangener Untersuchung der stationaren Schwingungen 
/sollen nun die Einschwingvorgange in diese stationéren Endzustande hinein untersucht werden. 
| Dabei wird insbesondere die Abhangigkeit des Endzustandes von der Art der Anfangsstérung 
aufzudecken sein. 
| Seiner gréBeren Ubersichtlichkeit wegen wird wieder der ungedimpft schwingende Stab zu- 
jerst behandelt, obgleich erst die Einbezichung dimpfender Krafte ein eindeutiges Bild der Er- 
» scheinungen liefert. 


| a) Ungedimpfte Schwingungen. Zur naherungsweisen Berechnung der Ubergangsschwin- 
| gungen wird an das Differentialgleichungssystem (18) angekniipft und dieses in der Form der 
/ zeitfreien Einzelgleichung (19) ausgewertet. Dies ist fiir den ungedimpft schwingenden Stab leicht 
jméglich, da sich (19) fix 6 = 0 streng lésen 14Bt. Nach kurzer Rechnung ergibt sich mit der 


) Integrationskonstanten C die Lésung in impliziter Form zu: 


poy (ae + By +21 cals EBL o aa Bea CO, (31) 


4 2 


| Die zahlenmaBige Auswertung dieser Lésung fiir ein Beispiel mit vier verschiedenen Werten der 
} Langserregungsfrequenz ergibt die in den Abb. 5, 6, 7 und 8 gezeichneten (A, B)-Kurven?. Hierin 
i gehdren die verschiedenen (A, B)-Kurven’ einer Abbildung verschiedenen Werten der in (31) vor- 
{kommenden Integrationskonstanten C zu. Die Zeit ist allerdings in den (A, B)-Kurven nur als 
tunbekannter Parameter enthalten, so da® sich nur der Durchlaufungssinn fiir die Kurven an- 
lgeben laBt. Dieser ist auf den (A, B)-Kurven durch Pfeile in Richtung wachsenden Zeitpara- 
meters kenntlich gemacht. 
| Durch den Ansatz (13) 1a&t sich allgemein eine beliebige, kleine Anfangsstérung beschreiben. 
‘Denn A(t) und B(r) sind Losungen des Differentialgleichungssystems (18) und bringen daher je 


1 Bei den numerischen Auswertungen hat mich Herr H. Bése wirksam unterstiitzt, was durch eine Zu- 
wendung der Karlsruher Hochschulvereinigung e. V. in dankenswerter Weise erméglicht wurde. 


62 Weidenhammer: Das Stabilitatsverhalten der nichtlinearen Biegeschwingungen _Ingenieur-Archivy 


ie 


4A 


Abb. 5. (A, B)-Kurven fiir m = w,; ohne Dampfung. Abb. 6. (A, B)-Kurven fiir @ = 2 O73 ohne Dimpfung. 


c> 


A | 


4A 


= 


\bb. 7. (A, B)-Kurven fiir = ,; ohne Diimpfung. Abb. 8. (A, B)-Kurven fiir @ > 2 O55 ohne Diimpfung. 


& 


’ 
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eine Integrationskonstante in W(z) ein. Somit kénnen vorgegebene Werte der Anfangsauslenkung 
und der Anfangsgeschwindigkeit erfaBt werden und damit hat der Ansatz (13) zu einer hinreichend 
allgemeinen Lésung gefiihrt, die in der Form der (A, B)-Kurven ausgewertet werden kann. 
Abb. 5 gehért zu einer Anregungsfrequenz (w,, Abb. 9 oben), welche gerade auf der einen Grenze 
des Instabilitatsbereiches liegt, und als stationdre Bewegung nur die gestreckte Lage zulabt 
(Schwingungsamplitude Null, Abb. 9 unten). Diese Stabilitatsaussage liest man unmittelbar aus 
dem geometrischen Bilde ab, denn die gestreckte Stablage ohne Querschwingungen wird durch 
den einzig vorhandenen singularen Punkt 4, = B, = 0 wiedergegeben und dieser Punkt ist ein 
isolierter Punkt. Offenbar verlaufen demnach bei kleiner Anfangsstérung alle Schwingungen als 
Kleine Bewegungen um die Ruhelage, so da®B diese selbst stabil ist. 
In Abb. 6 ist eine Frequenz genau in Mitten des Instabilitatsbereiches (wo =2;, Abb. 9) 
angenommen, wozu eine instabile gestreckte Lage und eine stabile Querschwingung endlicher 
‘Amplitude gehéren (Abb. 9). Ersichtlich gibt sich die Instabilitat dadurch zu erkennen, da® der 
Punkt A, = B, = 0 gegeniiber Abb.5 seinen Charakter 


ie - . . 
gedndert hat, indem er nun ein Doppelpunkt geworden 
ist und daher alle in seiner Nahe verlaufenden Kurven 
ee eine sich entfernende Tendenz aufweisen. Hingegen sind 
a tAcost 
a= _@ Druck fe 
As + Be Lug a 
<—Se = - Se SS 
€=const Ww 
b ee 
Ww 
c 
es. o g 
Oe 20; W; D>20; 
| Abb. 9. (oben) Instabilititsbreiche. (unten) Resonanz- Abb. 10. a Zeitlicher Verlauf der pulsierenden Axiallast. 
/kurve, untersuchte Frequenzen durch © bezeichnet. b Instabile Querschwingung; sin 7/2-Typ. c Stabile Querschwingung: cos 1/2-Typ. 


|zwei weitere singulare Stellen als Einsiedlerpunkte in das Bild gekommen. Sie gehéren wegen 
A, += 0, B) = 0 zu den halbperiodischen Querschwingungen vom cos-Typ, deren Stabilitat aus 
der isolierten Lage dieser Punkte zu erkennen ist. Allerdings gibt Abb. 6 den ersten Hinweis 
idarauf, daB die Einbeziehung dimpfender Krifte wiinschenswert ist. Denn Schwingungsbewegun- 
igen mit derart schwingenden Amplituden und Phasen, wie sie die Abb. 6 fiir gréBere Betrage 
tvon A und B nicht ausschlieBt, sind wenig plausibel und tatsachlich auch nur wegen der bedenk- 


Hlichen Vernachlassigung dampfender Krafte aufgetreten. 


| In der Abb. 7 ist der Abb. 5 entsprechende Grenzfall aufgetragen. Die Frequenz liegt auf der 
jgegeniiberliegenden Stabilitatsgrenze (w,, Abb. 9) und aus dem Amplitudendiagramm liest man 
jab, daB hier die gestreckte Lage ohne Querschwingungen nicht stabil ist, jedoch stabile Quer- 
ischwingungen mit endlicher Amplitude méglich sind. Die Stabilitatseigenschaften lassen sich 
jwiederum aus der jeweiligen Natur der Singularitaten ablesen. 

_ In Abb. 8 ist der hier besonders interessierende Fall ausgewertet, der bei einer Anregungs- 
Hfrequenz w > 2, auf zwei stabile Bewegungsformen gleichzeitig fihrt. Die Stabilitatskarte 
\(Abb. 9 oben) sagt fiir diese Frequenz nur aus, daB die Stabilitat der gestreckten Lage ohne Quer- 
ischwingungen besteht, wabrend aus dem Amplitudendiagramm (Abb. 9 unten) stabile Schwin- 
igungen mit der Amplitude Null (Ruhelage), aber auch stabile Querschwingungen mit endlicher 
| Amplitude zu entnehmen sind. Entsprechend vielgestaltig ist der Verlauf der Kurven in Abb. 8, 
Jiiber den man zunichst schon sagen kann, dab auf Grund eines Satzes von H. Poincaré! die drei 
Wirbelpunkte notwendig zwei Sattelpunkte mit in das geometrische Bild der (A, B)-Kurven 
tbringen miissen. Wahrend offenbar jede Stabschwingung, die durch eine Stérung mit den ent- 
sprechenden Anfangswerten eingeleitet wird, ganz in der Nahe der Amplituden- und Phasen- 
Iwerte von A, = B, = 0 nach (20,1) oder By = 0, Ay + 0 nach (20,2) verlauft, entfernt sich jede 
iSchwingung in der (A, B)-Ebene aus der Nahe der durch (20,3) 4) = 0, By + 0 gegebenen in- 
istabilen Schwingung. Diese merkwiirdige Instabilitat der sin-Schwingung hatte sich bereits 


1 Vel. z. B. S. 48 des in der FuBnote auf S. 60 zitierten Buches von J. J. Stoker. 
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friiher ! mit der Methode der kleinen Schwingungen nachweisen lassen und ergibt sich hier formal 
daraus, da® der durch (20,3) gegebene Punkt ein Doppelpunkt der (A, B)-Kurven ist und sich 
daher alle Nachbarkurven von ihm 
entfernen. Die Instabilitatseigen-/ 
schaften der sin-Schwingung 1aSti 
sich insofern anschaulich verstehen’ 
(Abb. 10) als bei der in (3) gewahl- 
ten Vorzeichenfestsetzung fiir das: 
Potential V(r) der Langskraft 

&x cost > 0 | 


eimen Langsdruck im Stab bezeichnet 
und daher eine Querschwingung vom 
sin t/2-Typ immer gerade in die ge- 
streckte Lage zuriickgezogen wird 
wenn eine Querauslenkung vorhandem 
ist. Wie aus Abb. 10 fiir den eben- 
falls eingezeichneten cos t/2-Schwin- 
gungstyp zu ersehen ist, tritt dann 
eine andersartige Phasenbeziehungs 
ein, so daB der Stab im ausgelenktent 
Zustand gedriickt wird und daher eine} 
Querschwingung endlicher Amplitude 
aufrechterhalten werden kann. 


Auch in Abb. 8 sind fiir gréBere# 
Betrage von A und B in Kreise iiber 
gehende Kurven zu sehen. Diese ent 
sprechen ungedimpften, fastperi- 
odischen Schwingungen, die jedoch 
wenig plausibel sind und die auch bei 
Bericksichtigung dampfender Kraftet 
verschwinden. Allerdings wird (19) 
fir den gedimpft schwingenden Stal 
nicht mehr geschlossen zu lésen sein 
und dadurch die Ubersichtlichketif 
verringert, weshalb die durchsichtigr} 
und in sich geschlossene Theorie des} 
ungedampft schwingenden Stabess 
vorangestellt wurde. 


b) Gedampfte Schwingun+ 


gen. Die Betrachtung des un-+ 


gedampft schwingenden Stabes er+ 
gab insofern kein véllig befriedigen-; 
des Bild von den Einschwingvor-. 
gangen, als sich immer noch Be-!. 


wegungen ergeben, die nicht in eine; 


stationaére Schwingung oder in ie 


Ruhelage einmiinden. Die Einbezie-} 
hung dampfender Krafte erweist 
sich daher als notwendig und < | 
reits der einfache geschwindigkeits- 
proportionale Ansatz fiir die Damp- 
Abb. 12, (4, B)-Kurven fir @ > 2: mit Diimpfung. fungskraft ergibt eindeutige Ver-} 
haltnisse in dem Sinne, da8 danni 

zu jeder Anfangsstérung ein stationarer Endzustand anzugeben ist, inden der Stab einschwingt. 


' Vel. S. 326, Abschn. f) der in FuBnote 1 auf S. 53 genannten Arbeit. 
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Die Auswertung der Differentialgleichung (19) hat nun fiir 6 + 0 zu geschehen und ist offen- 
bar formelmafig nicht mehr méglich. Die Integralkurven kénnen jedoch fiir den vorliegenden 
Fall geniigend genau mit Hilfe des Isoklinenfeldes gezeichnet werden. Das Ergebnis dieser 
numerischen Auswertungen fiir ein Zahlenbeispiel mit den bereits beim ungedaimpft schwingenden 
Stab verwendeten Zahlenwerten ¢, y,ax und nun jedoch kleinem f + 0 ist in den Abb. 11 und 12 
aufgetragen. Dabei entspricht Abb. 11 der in Abb. 6 untersuchten Frequenz wm = 2; in Mitten 
des Instabilitatsbereiches, wihrend Abb. 12 dem in Abb. 8 aufgetragenen Falle fir w > 2 (Dj 
entspricht. Hingegen ist auf die Untersuchung der Frequenzen auf den Bereichsgrenzen (w, 
und w, aus Abb. 9) hier verzichtet worden. 

Fir m = 2 w; entnimmt man der Abb. 11 die Instabilitat der gestreckten Lage und die Stabili- 
tat der Querschwingungen endlicher Amplitude. Die Instabilitatseigenschaften der gestreckten 
Lage A) = 0, B, = 0 ergibt sich im geometrischen Bilde wiederum als Sattelpunkt zu erkennen, 
wahrend die stationare Querschwingung sich durch zwei Strudelpunkte mit A, + 0, By ++ 0 mit 
kleinem Betrag von B, anzeigt. Die Berechnung dieser Werte als singulire Stellen von (19) wurde 
in Abschn. 4b bereits erledigt und es geniigt hier darauf hinzuweisen, da®B man auf Grund der 
dampfungsbedingten Phasenverschiebung keine reine cos-Schwingung erwarten darf, sondern durch 
B, + 0 auch einen sin-Anteil im Ansatz (13) erhalt. 

Fir © > 2; ist in Abb. 12 der wichtige Fall aufgetragen, in dem auSer der Ruhelage zwei 
Querschwingungen stabil sind. Im geometrischen Bilde entspricht dies wiederum drei singularen 
Punkten, die als Strudelpunkte stabil sind. Die beiden weiteren Schwingungen, die dem sin-Typ 
im Falle des ungedimpft schwingenden Stabes entsprechen, sind auch hier als instabile Sattel- 
punkte zu erkennen. In den Abb. 11 und 12 ist nunmehr jeder beliebigen Anfangsstérung ein- 
deutig ein stationaérer Endzustand zuzuordnen, auf den sich der Stab im Laufe der Zeit einschwingt. 
Und zwar trennen diejenigen speziellen (A, B)-Kurven, welche durch die Sattelpunkte laufen, 
Gebiete der (A, B)-Ebene voneinander ab, deren Kurven simtlich auf den gleichen singuladren 
Punkt zulaufen. Jeder der genannten speziellen Kurven wird als Separatrix bezeichnet und diese 
sind in den Abb. 11 und 12 gestrichelt eingezeichnet. Offenbar geben diese trennenden (A, B)- 
Kurven in Abb. 12 die Entscheidung dariiber, ob eine spezielle Anfangsstérung wieder in die 
Ruhelage abklingt oder sich auf eine der beiden méglichen Querschwingungen gleicher endlicher 
Amplituden aufschaukelt. Demnach hangt es also von der Anfangsamplitude und Anfangsphase 
einer Stérung ab, ob sich eine Querschwingung einstellt oder nicht. Die in Abschn. 1 aufgeworfene 
Frage a) kann damit fiir jeden Einzelfall auch quantitativ beantwortet werden. Die (4, B)-Kurven 
enthalten die Zeitvariable nur als Parameter. Um den Zeitablauf kenntlich zu machen, sind auch 
in Abb. 11 und 12 Pfeile eingetragen, doch kann iiber den ZeitmaBstab auf den Kurven selbst 
kaum etwas ausgesagt werden. Wahrend die Strudelzentren mit unbegrenzt wachsender Zeit an- 
genahert werden, entsprechen die Sattelpunkte fiir jeweils zwei ihn durchlaufende Kurvenaste 
dem Wert t > + ©. wahrend die restlichen beiden nur fiir 7 > — © auf den Sattelpunkt fiihren'. 
Daraus folgt, daB der Zeitmafstab der in ihrer unmittelbaren Nahe verlaufenden Kurven sehr 
gedrangt ist und demgemafB die Kurveniste nur in langeren Zeitraumen durchlaufen werden. Da 
aber aus dem zeitfreien System tiberhaupt nur Aussagen zu gewinnen sind, deren Fehler mit 
wachsender Zeit zunimmt ”, sind hier keine scharfen Aussagen zu erwarten, was aber im vor- 
liegenden Falle unwesentlich ist, da die zugehérige Schwingung ohnehin instabil ist. 


6. Berechtigung der Naherungsrechnung. An Stelle des vollstandigen Systems der Schwin- 
gungsgleichungen (10) wurde bisher nur jeweils die eine verbleibende Einzelgleichung (12) mit 
der Eigenfrequenz ~; ausgewertet, wenn Schwingungen mit Erregerfrequenzen w ~ 2; unter- 
sucht werden sollten. Es mu®B daher noch gezeigt werden, da dieses Vorgehen in dem Sinne 
herechtigt ist, daB die bei der Einzelgleichung gefundenen Phanomene auch in den Schwingungen 
des gesamten Systems wiederzufinden sind. Fiir diesen Nachweis soll die von J. Haag * auch auf 
Schwingungssysteme mit mehreren Freiheitsgraden ausgedehnte Methode der langsam verander- 
lichen Phase und Amplitude verwendet werden. Zu dem Zweck wird das vollstandige System (10) 


in der Form 


F lv. & Ps 
Wi +a W;= a — Bq; Wj +0; cost Wj —y;; Wj) —W; Sr) si Nae 
i ie We (je 2 oN) 2 (10,1) 


1 J. Haag, Les Mouvements vibratoires, Bd. 1, S. 135, Paris 1952. 
2 Vgl. die in FuBnote 1 auf S. 58 genannte Arbeit von J. Haag. 
3 J, Haag, Ann. scient. de I’Ecole Norm. Sup. 3, 64 (1947), S, 285. 
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unter Beibehaltung der Abkiirzungen (11, 1—6) geschrieben. Mit den natiirlichen Zahlen m;, M 
und dem kleinen Parameter ¢ wird zunachst eine Naherung fiir 1/q; durch 


Lie oats ~ 
Pepe de (Fa, 2a ay) (32) 
eingefiihrt, worin die /; 2 0 reelle Konstanten bezeichnen. Dann kann man an Stelle von (10,1) 
schreiben 
‘ie m; 2 my; : . 
wy + () w= e(G) am +, 1 (33) 


wenn rechtseitig wie immer im folgenden nur erste Potenzen von é beibehalten werden und W 
und W' fiir die Gesamtheit der in F; vorkommenden Funktionen W; steht. Die gesuchten Zeit- 
funktionen W; werden nun mit den Amplituden y,(r) und den Phasen x,(z) formal etwas anders 


als in (13) in der Gestalt 
my; 


W; = y; cos (is + x) = y; cos D; 
angesetzt. Die Langsamveranderlichkeit von y; und x; soll durch die zusatzliche Forderung 


Mm; 
WV; =— Ta sin D; (34) 


zum Ausdruck kommen, so dah also 
— y; x; sin®; + y;cos®, = 0, (7 = lees ees) (35) 


sein muB. Bildet man durch Differentiation von (34) die Ableitungen WwW; und setzt diese in (33) 
ein, so erhalt man 


m, | 
— y; x; cos ®; — y; sin ®; = oes y; cos B; + Fy cos®, tae sin QD, | >» Aj asl 2 ake 


(36) 
was sich zusammen mit (35) nach x; und y; auflésen laht. Diese Auflésung ergibt das Differential- 
gleichungssystem erster Ordnung 


m; | m 5 
a D yj; cos O; — Fy cos @, —— ysin DM, r) 


yj Xj =n cos ®D; ‘ (37,1) 


sn@®, (j=1,2,...,N),. (87,2) 


r AUG F é ES @ m : 
Yj i M Aj Mi cos @; meee ay, x cos ae Mo sin ®, ia 


worin (37,1) im allgemeinen den AusschluB der Werte y; = 0 verlangt. Weil im vorliegenden Falle 

jedoch die Funktionen F; jeweils den Faktor y; abspalten, gilt hier (37,1) ohne diese Einschrin- 

kung. Da in (37) die Zeit explizit auftritt, geht man zweckmaBig durch Mittelung iiber eine 

Periode der rechten Seiten zum zugeordneten zeitfreien System tber, denn stationire Schwin- 

gungen sind niherungsweise auch aus diesem vereinfachten zeitfreien System noch zu berechnen. 
Das zeitfreie System 


sjaes, Vp esl a (eee (38) 
ergibt sich nach Ausfiihrung bestimmter Integrale mit den rechten Seiten 
21M 
A; = eae | a i y, cos O, — F; (y cos@, — ar y sin ®, ) cos @; dt , 
aM (37) 
Vie 5 a | ar A; y; 00s D; — F; (y cos @, — 7 y sin @, i) sin O; dr, 
0 


worin fiir die Mittelwertbildung unter dem Integralzeichen y; und x; wie Konstante zu behandeln 
sind. 
Da in den untersuchten Fallen die Erregerfrequenz annihernd gleich dem doppelten Wert 


einer Kigenfrequenz «; der Biegeschwingungen ist, hat man nach (32) die ganzen Zahlen m, und 


| 
| 
| 
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M so zu wahlen, daf fiir den festen Zeigerwert j gilt: 


oe) 


77 
Dann ist auf Grund der als verschieden. vorausgesetzten Kigenfrequenzen 


pe a eal 2 ne Ne) 


von selbst erfiillt. Nach Auswertung der Integrale (39) liefert (38) dann: 


eae iT AS 0; 3 Lig oo 
sea (2 4 008 2 ay tei IG + me rst) (40,1) 
+i 
1 r mM}, Ay 3 i a 9 
oe ar (3 a g vkR er a = Vkm 1) ? (40,2) 
Lara 
ee mj aj : Bj 
Pa \ ype 2 %;— y Yi}> (40,3) 
my, 
Y, =i | Ss 2) (RTE Bee (40,4) 
Fiir eine stationdre Schwingung miissen alle Amplituden y; und Phasen x;(j = 1,2,..., N) 


konstante Werte annehmen, was offenbar dann der Fall ist, wenn mit den verschwindenden linken 
Seiten von (38) auch alle rechten Seiten verschwinden, und dies sind die Gleichungen: 


l (eece 
a ee yoke 9, Yj=0, Ye=0 (k=1,2,...,N,4j). (40,14) 


Die hier interessierende Loésung der Gleichungen (40) ist diejenige, bei welcher Y;, = 0 (40,4) 
| durch y, = 0 fiir k = 1,2,...,.N, + j und Y; = 0 (40,3) durch y; + 0 befriedigt wird. Wenn 


so alle Amplituden mit Ausnahme yon y; verschwinden, bleibt offenbar in ke X; = 0 (40,1) 
Vi 


und Y; = 0 (40,3) nur y; und x; allein stehen, was genau der in Abschn. 4 und 5 geiibten Be- 
schrankung auf die Einzelgleichung mit dem festen Index] entspricht. Die Aussagen iiber die 


Phasen x; aus — X;, = 0 (40,2) sind hier ausnahmsweise gegenstandslos, da ja die zugehérigen 


Amplituden y, simtlich verschwinden. 

Allerdings mu die gefundene Lésung y;+ 0, y, = 0 fir k = 1,2,...,N, +] noch als 
stabil nachgewiesen werden. Die Methode der kleinen Schwingungen gestattet die Stabilitat der 
Lésungen y, = 0 einfach nachzuweisen, indem man die Nachbarbewegung 7, zu y, = 0 in (38) 
in (40,4) einsetzt. Man erhalt 


=B,, . 
1] = — est ne (k= Me roca ATS J) 
was sich unmittelbar durch 
my, Bh, 
nk = Cy e “™M 2 


mit den Integrationskonstanten C, lésen laBt. Demnach klingt jede etwa eingeleitete Stérung 
exponentiell ab und die Lésung y, = 0 fir k = 1,2,..., N, = jist daher stabil. Da die Phasen- 
-aussagen (40,2) nicht interessieren, bleibt nur noch die Priifung der Stabilitat der Losungen y; + 0 
aus (40,1) und (40,3). In diesen Gleichungen kommt jedoch nur noch y;, x; allein vor, und zwar 
+so, wie man es auch aus der Einzelgleichung allein herleiten wiirde. Die Stabilitatsdiskussion fiir 
‘die Einzelgleichung ist jedoch mit dem Ansatz (13) in den Abschn. 4 und 5 ausfiihrlich durch- 
gefiihrt worden, so da® der Nachweis bereits gefiihrt ist. Tatsdchlich lefert der Ansatz mit der 
‘Amplitude y; und der Phase x; die gleichen Stabilitatsaussagen. Fir die zeichnerische Dar- 
)stellung der Ergebnisse hatte jedoch die (A, B)-Ebene gewisse Vorziige. 

Anschaulich ist das Schwingungsverhalten des Systems (10) durchaus verstandlich. Wenn die 
| Erregerfrequenz w ~ 2 w, ist, so hat nur das Schwingungsglied mit der Eigenfrequenz w; die zur 
' Erregung erforderliche Abstimmung und damit cine erregende Kraft [vergleiche (40,1), (40,3)], 


*wihrend alle anderen Schwingungsglieder keine Erregerkrafte enthalten [vergleiche (40,2), (40,4)]. 
5x 
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Da aber in allen Fallen dampfende Krifte wirken, miissen notwendig die Amplituden der ver- | 
stimmten und damit nicht erregten Systemteile abklingen, so da im eingeschwungenen Zustand | 
das Schwingungsverhalten lediglich durch den Schwingerteil mit der Eigenfrequenz w,; allein be- 


stimmt wird. Die iibliche Beschrinkung auf die Untersuchung dieses eingliedrigen Schwingers 
diirfte daher durchaus berechtigt sein, wenn man Biegeschwingungen berechnen will, die mit der 
halben Frequenz der Langserregung verlaufen. 


7. Zusammenfassung. Der axial pulsierend belastete Stab kann in Querschwingungen endlicher 
Amplitude iibergehen, wenn die Lingserregungsfrequenz in der Nahe des doppelten Wertes einer 


seiner Kigenfrequenzen der Querschwingungen liegt (w ~2.,). Die dann eintretenden, nicht- | 


linearen Schwingungen wurden fiir den momentenfrei gelagerten Stab untersucht. Insbesondere 


wurden die Kinschwingvorginge mit der Methode der langsam veradnderlichen Phase und Ampli- | 
tude berechnet, um die Abhangigkeit der stationéren Schwingungen von einer Anfangsstérung in | 
Querrichtung aufzufinden. Dabei ergab sich u.a. fiir Langserregungsfrequenzen jw > 2 ~; die | 
Méglichkeit des Zuriickschwingens in die Ruhelage wie auch des Einschwingens auf Querschwin- | 


gungen endlicher Amplitude. Weiter wurde die Resonanzkurve fiir langserregte Querschwingungen 


unter Beriicksichtigung dampfender Krafte berechnet. Sie verlauft wie bei den erzwungenen | 


nichtlinearen Schwingungen ganz im Endlichen. Da die genannten Ergebnisse wie tblich nur 


aus einer Schwingungsgleichung hergeleitet wurden, war die Berechtigung dieser Naherungsrechnung | 


noch besonders nachzuweisen. Mit einer Theorie von J. Haag lieB sich dieser Nachweis fir ge- 
dimpfte Schwingungen leicht in dem Sinne fiihren, da die gefundenen Phanomene auch im voll- 
standigen System der Schwingungsgleichungen wiederzufinden waren. 


(Eingegangen am 6. Juli 1955.) 
Aus dem Institut fiir Mechanische Schwingungstechnik der Technischen Hochschule Karlsruhe. 
Anschrift des Verfassers: Priv.-Doz. Dr. F. Weidenhammer, Karlsruhe i. B., HertzstraBe 16, Bau 40. 
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